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1 Einleitung 


Blech ist ein weit verbreitetes Material in der Produktion von Gütern. Die 
Karosserien von Autos, die Oberflächen in gewerblichen Küchen, Compu- 
tergehäuse und Schiffsrümpfe sind beispielsweise alle aus Blech geformt. 
Bleche können aus unterschiedlichsten Materialien hergestellt werden. Am 
verbreitetsten sind Bau- und Edelstähle, aber auch Aluminium, Kupfer, 
Messing oder sogar Titan sind als Blech verfügbar [10]. 

Für die Bearbeitung von Blechen zu einem Endprodukt ist fast immer 
ein trennendes Fertigungsverfahren notwendig. Die ersten trennenden 
Verfahren waren spanabtragende Verfahren wie z. В. Sägen und Nibbeln. 
Die Bearbeitung mit mechanischen Verfahren ist aber relativ langsam 
und benötigt oft Nacharbeit, z.B. das Verrunden von scharfen Kanten am 
Bauteil. In der zweiten Hälfte des 20. Jahrhunderts wurden Verfahren wie 
das Wasserstrahl-, Plasma- und Laserschneiden entwickelt [39, 69]. 

Seit der Markteinführung von Laserschneidmaschinen für die trennende 
Bearbeitung in der flexiblen Blechfertigung in der Mitte der 1980er Jahre 
haben sie einen Großteil des Marktes erobert. Dies liegt an ihrer hohen Ge- 
schwindigkeit, oftmals geringen Nacharbeit und sehr geringen Rüstzeiten. 
Das Werkzeug muss nicht gewechselt werden und erleidet im Vergleich 
zu spanenden Verfahren sehr wenig Verschleiß. 

Beim Laserschneiden wird das Blech durch einen Laserstrahl von oben 
sehr schnell stark erhitzt. Das Material schmilzt und wird durch ein Pro- 
zessgas, і. d. К. Sauerstoff oder Stickstoff, nach unten aus dem Schnittspalt 
ausgetrieben [44, 75]. Der schematische Aufbau eines Schneidkopfs und 
seine Lage zum Blech sind in Abbildung 1.1 dargestellt. 

Das Prinzip des Laserschneidens ist in unterschiedlichen Maschinenkon- 
zepten in den Markt eingeführt worden. Diese unterscheiden sich darin, 
welche Bewegungen durch den Schneidkopf oder das Blech ausgeführt 
werden. Bei C-Rahmen-Maschinen (siehe Abbildung 1.2 links) wird das 
Blech in zwei Achsen bewegt und der Schneidkopf verfährt ausschließ- 
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Abbildung 1.1 Schematischer Aufbau des Schneidkopfs beim Laserschneiden. 


lich in der Höhe. Bei O-Rahmen-Maschinen (siehe Abbildung 1.2 mitte) 
übernimmt der Schneidkopf zusätzlich eine Achse der Bewegung in der 
Blechebene, sodass das Blech nur in einer Achse bewegt wird. Neben die- 
sen beiden Maschinenkonzepten haben sich vor allem Flachbettmaschinen 
(siehe Abbildung 1.2 rechts) durchgesetzt, bei denen das Blech zur Bear- 
beitung still liegt und sich der Schneidkopf in drei Dimensionen bewegt 
[41]. 

Damit die Schlacke nach unten aus dem Schnittspalt ausgeblasen wer- 
den kann, sollte der Bereich unter dem Blech frei sein. Dies widerspricht 
aber einer stabilen Lagerung des Blechs und der zu schneidenden Teile. 
Gelöst wird dies durch feste Auflagestege (siehe Kapitel 1.1). Nebenzeiten 
können bei der Laserflachbettmaschine (LFM) dadurch verringert werden, 
dass die Maschine zwei Paletten hat, von denen eine bearbeitet werden 


Abbildung 1.2 Schematische Darstellung der verschiedenen Maschinenkonzepte zum 2D- 
Laserschneiden. Links: C-Rahmen-Maschine. Mitte: O-Rahmen-Maschine. Rechts: Flachbett- 
maschine. 


1.1 Problemstellung 


kann, während die andere Palette außerhalb der Maschine für den nächsten 
Schneidvorgang gerüstet wird. Bei C-Rahmen- und O-Rahmen-Maschinen 
muss auch freier Raum unter der Schneidkopfpositionen vorhanden sein, 
um die Schlacke abführen zu können. Dies wird durch eine zylinderförmi- 
ge Öffnung sichergestellt. Um die Bewegung des Blechs zu ermöglichen, 
werden oftmals Bürstenteppiche aus Plastik genutzt. Bei beiden Maschi- 
nenkonzepten wird i.d.R. versucht, ein geschnittenes Teil sofort aus dem 
Restgitter zu entfernen, um Probleme bei der Bewegung des Blechs zu 
umgehen, da gekippte Teile sich am Bürstentisch verkeilen könnten. Die- 
se Unterbrechungen verlängern konzeptbedingt den Schneidvorgang im 
Vergleich zu einer LFM. 

Um auf einer heute marktgängigen LFM ein Teil ausschneiden zu kön- 
nen, sind drei Schritte nötig [10]. Erstens muss die Geometrie des Teils 
in digitaler Form vorliegen. Es folgt die Programmierung einer Tafel für 
die Bearbeitung auf der LFM. Dafür müssen alle zu schneidenden Teile 
auf dem Blech angeordnet werden. Dieser Prozess wird im Folgenden 
Schachtelung genannt. Nach der Schachtelung werden die Bearbeitungs- 
reihenfolge der Teile und die Verfahrwege des Schneidkopfs festgelegt. 
Damit ist die Logik des Programms festgelegt und kann in die Sprache der 
Maschinensteuerung übersetzt werden. Ein solches NC-Programm (engl. 
für numerical control) wird auf der Maschine ausgeführt. 
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Um das Blech während des Schneidens auf einer LFM zu lagern, werden 
von allen Herstellern wenige Millimeter dicke Auflagestege benutzt, die 
fast immer aus vielen dreieckigen Spitzen bestehen (siehe Abbildung 1.3). 
Diese Stege werden parallel zueinander in der Palette befestigt und sind 
oftmals aus Baustahl, können aber auch aus Edelstahl oder Kupfer gefertigt 
sein. 

Während des Schneidprozesses kann es zu Einbußen bei der Prozesssi- 
cherheit kommen, die nicht direkt mit dem Schneidprozess zu tun haben. 
Dies äußert sich vor allem durch folgende Punkte (siehe auch Abbildung 
1.4): 
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Abbildung 1.3 Beispiele von Auflagestegen an Laserflachbettmaschinen verschiedener Her- 
steller. Oben links: Amada. Oben rechts: Bystronic. Mitte links: LVD. Mitte rechts: Mazak. 
Unten links: Salvagnini. Unten rechts: TRUMPF. 


= Verkippen von Teilen durch die auf das Teil einwirkende Kraft des 
Prozessgases und unzureichende Lagerung durch Auflagespitzen. 


a Verschleiß der Stege vor allem durch vollständiges oder teilweises 
Abschmelzen von Spitzen, aber auch durch Anhaftung von Schlacke. 


= Verschweißen von Teilen mit einer Auflagespitze, wenn direkt über 
dieser Spitze geschnitten wird und die Schlacke nicht ausreichend 
durch das Gas entfernt wird. 


a Qualitätsprobleme, z. В. Schmauchspuren oder Schlackespritzer an 
der Unterseite eines Teils beim Schneiden über Auflagestegen. 
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Abbildung 1.4 Beispiele von unerwünschten Effekten der Stege. Links: Verkippte Teile. 
Mitte: Verschweißter Schneidbutzen auf einer Stegspitze. Rechts: Schlackespritzer am Teil 
nach Schnitt über einer Stegspitze. 


= Abriss des Schnittspalts beim Schneiden über einer Auflagespitze, 
wenn die Schlacke nicht ausreichend entfernt werden kann und die 
Laserleistung nicht ausreicht, um das Material trotz der im Schnitt- 
spalt verbleibenden Schlacke zu erwärmen. 


Die Auswirkungen dieser Punkte sind eine geringe Produktivität der 
Maschine, da es während der Bearbeitung eines NC-Programms zu unge- 
wollten Stillständen kommen kann. Um eine hohe Maschinenproduktivität 
zu ermöglichen, ist dann aber ein erheblicher Einsatz von Personal nötig, 
das diese Probleme erkennt und entweder verhindert oder schnellstmög- 
lich löst. Außerdem kann nur unzureichend auf Automatisierungslösungen 
zurückgegriffen werden, da Teile, die verkippt oder mit einer Stegspitze 
verschweißt sind, nicht automatisch entnommen werden können. 

Es ist auffällig, dass alle oben genannten Punkte mit den Stegauflagen 
in Verbindung gebracht werden und grundsätzlich deutlich abgemindert 
oder sogar vollständig vermieden werden können, wenn die Position der 
Auflagestege relativ zum Blech bzw. der zu schneidenden Teile bekannt 
ist. Der Ansatz dieser Arbeit ist, die Position und ggf. auch den Zustand 
der Auflagestege visuell zu erfassen. 

Außerdem kann die Lage des Rohblechs relativ zu den Stegen gemessen 
werden. Basierend auf den zu schneidenden Teilen, deren Position auf dem 
Blech, der Blechdicke und -dichte sowie dem Gasdruck kann dann das 
Verkippen der Teile vor dem Schneiden berechnet werden. Ebenso können 
die Stegbeschädigung, die Neigung zum Verschweißen oder Orte von 
möglichen Qualitätsproblemen vorhergesagt werden. Das NC-Programm 
kann dann entsprechend angepasst werden, um die Probleme je nach 
Nutzerwunsch bestmöglich zu vermeiden. Dazu zählt z. В. das Verschieben 
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von Teilen auf dem Blech, das Verschieben der Einstech- und Anfahrpunkte 
eines Teils oder die Veränderung des Verfahrwegs des Schneidkopfs bei der 
Positionierung zwischen zwei Teilen. Der Lösungsansatz ist schematisch 
in Abbildung 1.5 dargestellt. 


Kapitel 2 
u. 
Kapitel 2 i 
api Pen Kapitel 3 
l ! der Stege f \ 
Mögliche Berechnung der 
Sensoren Auswirkungen 


Erkennung 
der Blechlage 


Verbesserung des 
NC-Programs 


LI 
Kapitel 3 


Abbildung 1.5 Struktur der vorliegenden Arbeit anhand des Mess- und Berechnungsablaufs 
zur Generierung von prozesssicheren NC-Programmen bei LFMn. 


Natürlich sind auch andere Ansätze zur Milderung der geschilderten 
Probleme bekannt. Zum Beispiel ist eine andere Strategie, um ein Verkippen 
von Teilen zu vermeiden, die sogenannten MicroJoints [83]. Hierbei wird 
die Kontur eines Teils nicht vollständig geschnitten und es bleibt eine kleine 
Brücke (i. d. R. 0,5 mm bis 1 mm) zwischen dem Restgitter und dem Teil 
stehen. Diese Verbindung reicht aus, um ein Teil im Restgitter stabil zu 
halten. Um das Verkippen aller Teile eines Blechs garantiert zu verhindern, 
kann jedes Teil mit einem oder mehreren MicroJoints stabilisiert werden. 
Dies erfordert jedoch bedeutend mehr Nacharbeit, da beim Absortieren 
der Teile von der Palette zusätzlich jedes Mal mit deutlichem Kraftaufwand 
der MicroJoint gelöst werden und ggf. verbleibendes Material abgeschliffen 
werden muss. 

Da die Probleme spezifisch für die LFM sind, lassen sie sich durch den 
Einsatz von C-Rahmen- oder O-Rahmen-Maschinen zu großen Teilen ver- 
meiden. Insbesondere das Verkippen tritt bedingt durch die bessere Auf- 
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lage des geschnittenen Teils auf dem Bürstentisch deutlich seltener auf. 
Dafür haben diese Maschinen ihre eigenen Prozesssicherheitsprobleme 
und erreichen bedingt durch das Maschinenprinzip nicht die Bearbei- 
tungsgeschwindigkeit und Produktivität einer LFM. 


1.2 Eigener Beitrag 


Der Beitrag dieser Arbeit liegt vereinfachend zusammengefasst in zwei 
Gebieten: Zum einen in der visuellen Detektion und Vermessung der Auf- 
lagestege und ihres Zustandes und zum anderen in der Auswirkung auf 
die Schachtelung von Teilen auf einem Blech unter Berücksichtigung der 
Position der Auflagestege und weiterer Nebenbedingungen, die aus dem 
Laserschneidprozess resultieren. In Abbildung 1.5 entspricht dies den 
Schritten, die den Kapiteln 2 und 3 zugeordnet sind. An dieser Stelle sollen 
die Verbesserungen gegenüber dem Stand der Technik kurz zusammenge- 
fasst werden: 


= Laser-Triangulation zur Stegdetektion: Es wird ein neuer Aufbau 
für ein Laser-Iriangulationssystem erarbeitet, der sich für LFMn be- 
sonders eignet, weil es zu keiner Abschattung und wenig störenden 
Reflexionen kommt. Für letztere werden Kompensationsverfahren 
entwickelt. Außerdem wird eine Möglichkeit zur Kalibrierung des 
Messsystems direkt auf dem zu messenden Objekt konzipiert. 


Stegdetektion auf einzelnen Bildern im sichtbaren Spektrum: Die aus 
dem Stand der Technik bekannten Methoden zur Objektdetektion 
werden auf einzelnen Bildern der Palette einer LFM getestet. Dabei 
wird festgestellt, welche Features am besten zur Detektion geeignet 
sind und wie sich die Störeinflüsse durch die Produktionsumgebung 
auf die Detektionsverfahren auswirken. 


Stegdetektion mittels schmalbandiger Nahinfrarot-Beleuchtung: Es 
wird eine schmalbandige Beleuchtung im Nahinfrarot-Spektrum 
(NIR) eingesetzt, um die durch die Störeinflüsse der Produktions- 
umgebung ausgelösten Probleme bei der Detektion der Auflagestege 
auf einzelnen Bildern zu umgehen. Hierbei werden die Stege in ei- 
nem flachen Winkel von hinten beleuchtet und die Kamera mit einem 
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Schmalband-Filter versehen, wodurch sehr dunkle Bilder entstehen, 
bei denen fast ausschließlich die Auflagestege hell sind. Es werden 
Verfahren entwickelt, mit deren Hilfe diese Bilder um mögliche Re- 
flexionen bereinigt sowie die Stege und Spitzen effizient erkannt 
werden können. 


Schachtelung unter Berücksichtigung der Auflagestege: Die Schach- 
telung von zu schneidenden Teilen wird bisher ausschließlich unter 
Berücksichtigung des entstehenden Verschnitts durchgeführt. Dies 
wird um die Verkippwahrscheinlichkeit der Teile und eine mögliche 
Stegbeschädigung erweitert. Es wird ein formales Modell entwickelt 
und anschließend werden zwei mögliche Lösungsverfahren vorge- 
stellt. Zum einen handelt es sich dabei um eine Konstruktionsheu- 
ristik, die eine Menge von Teilen direkt verschachtelt und dabei alle 
Kriterien je nach einer gegebenen Gewichtung berücksichtigt. Die 
Anordnung der Teile erfolgt dabei nach der Bottom-Left-Heuristik 
und platzierte Teile werden abhängig von der Gewichtung in beiden 
Richtungen nachträglich verschoben. Das zweite Verfahren ist eine 
Verbesserungsheuristik, die die Schachtelung einer kommerziellen 
Software-Bibliothek als Ausgangslösung nutzt und dann nachträg- 
lich versucht, jedes Teil in eine bessere Position zu verschieben. Da 
das Blech in der vollen Breite belegt ist, werden die Teile bei diesem 
Verfahren nur entlang der Länge des Blechs verschoben. 


1.3 Gliederung der Arbeit 


Die Arbeit teilt sich in zwei große Kapitel: eines zur visuellen Detektion 
der Auflagestege und eines zur Schachtelung unter Berücksichtigung der 
Auflagestege. 

In Kapitel 2 werden zuerst die Grundlagen der Bildverarbeitung und 
wichtige Umgebungsbedingungen für die vorliegende Messaufgabe er- 
läutert. Dann werden drei Methoden zur Stegerkennung genauer betrach- 
tet, nämlich Lasertriangulation, die Stegdetektion auf einzelnen Bildern 
im sichtbaren Spektrum und die Stegdetektion mittels Nahinfrarot-Be- 
leuchtung. Für jede Methode wird zuerst der jeweilige Stand der Technik 
eingeführt, bevor die neuen Methoden und ihre Umsetzung in der Pra- 


1.3 Gliederung der Arbeit 


xis vorgestellt werden. Es folgt eine Validierung und Diskussion jeder 
Methode. 

In Kapitel 3 wird zuerst der Stand der Technik von Schachtelungsal- 
gorithmen vorgestellt, sowohl für die allgemeine Schachtelung als auch 
für den Spezialfall von Nebenbedingungen aus dem Produktionsprozess. 
Danach wird das Problem der Schachtelung unter Berücksichtigung der 
Auflagepunkte formal beschrieben und zwei Heuristiken vorgestellt, die 
zur Lösung eingesetzt werden können. Diese werden anschließend durch 
Berechnungen validiert und genauer charakterisiert. 

Abschließend fasst Kapitel 4 die Arbeit zusammen und gibt einen Aus- 
blick auf mögliche weitere Forschungsarbeiten. 


2 Vermessung der Auflagestege 
von Laserflachbettmaschinen 


Eine Vermessung der Auflagestege von LFMn ist auf viele Weisen um- 
setzbar. Als grundlegende Messprinzipien, die dafür eingesetzt werden 
können, sind 7. В. akustische, visuelle oder kapazitive Messungen denk- 
bar. Für eine visuelle Messung spricht vor allem, dass sich so eine große 
Fläche einfach in ihrer Gesamtheit oder zeilenweise erfassen lässt, wäh- 
rend für die anderen Prinzipien i.d.R. eine punktweise Messung nötig ist. 
Letzteres dauert oft deutlich länger, was zu unerwünschten zusätzlichen 
Nebenzeiten an der LFM führen würde. 

Im Folgenden werden zuerst die für dieses Kapitel notwendigen Grund- 
lagen der digitalen Bildverarbeitung eingeführt, die Umgebungsbedin- 
gungen an einer LFM diskutiert und mögliche visuelle Messprinzipien 
umrissen. Danach werden drei Messverfahren jeweils in einem Unterkapi- 
tel detailliert betrachtet. 


2.1 Grundlagen der digitalen Bildverarbeitung 


2.1.1 Grundlagen digitaler Bilder 
2.1.1.1 Aufnahme und Darstellung digitaler Bilder 


Die digitale Aufnahme eines Bildes geschieht mittels eines optischen Sen- 
sors, der aus einzelnen, rechteckig in einer Matrix angeordneten Zellen be- 
steht. Diese werden als Pixel (kurz für engl. picture element) bezeichnet. Eine 
Szene wird also diskretisiert wahrgenommen. In heutigen Digitalkameras 
werden dafür vor allem CMOS-Sensoren (kurz für engl. complementary 
metal-oxide-semiconductor) eingesetzt. Diese nutzen den photoelektrischen 
Effekt und mehrere Halbleitertransistoren, um das gemessene Signal zu 
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verstärken. Sie bieten im Vergleich zu Sensoren mit dem Charge Coupled De- 
vice Prinzip mehrere Vorteile, insbesondere hohe Geschwindigkeiten und 
geringe Kosten [7]. Die spektrale Empfindlichkeit eines CMOS-Sensors 
beginnt bei etwa 350nm Wellenlänge und endet bei ungefähr 1100 nm 
Wellenlänge. Damit decken sich CMOS-Sensoren weitestgehend mit dem 
wahrnehmbaren Spektralbereich eines Menschen, detektieren aber zusätz- 
lich auch Teile des NIR-Spektrums [46]. 

Je nach verwendetem Sensor in einer Kamera kann ein Grauwert- oder 
Farbbild aufgenommen werden. Zur Aufnahme eines Farbbildes wird 
der CMOS-Sensor mit spektralen Filtern vor jedem Pixel versehen. Der 
Wahrnehmung des menschlichen Auges folgend werden dabei drei Kanäle 
genutzt, nämlich rot, grün und blau. Jedes Pixel ist somit nur für einen 
Farbkanal nutzbar und die Information für die anderen beiden Kanäle 
muss für dieses Pixel interpoliert werden. Wenn der CMOS-Sensor ohne 
Modifikation betrieben wird, entstehen Grauwertbilder, die nur die Licht- 
intensität wiedergeben, aber keine Farbinformationen enthalten. Digitale 
Bilder werden oftmals in Datenformaten gespeichert, die das Bild als Ma- 
trix auffassen. Der Wertebereich für Grauwertbilder wird oftmals von 0 
bis 255 gewählt, was einer Darstellung mit 8 bit pro Pixel entspricht. 

Durch die Aufnahme mit einem Sensor, der aus diskreten Zellen besteht, 
lässt sich auch das Bild diskret darstellen. Ein Pixel kann über seine Koor- 
dinaten (u, v) identifiziert werden. Wie bei Matrizen üblich wird hier auch 
die erste Koordinate u für die vertikale Dimension genutzt und die zweite 
Koordinate v entsprechend für die horizontale Dimension. Der Grauwert 
eines Pixels ist durch die Funktion g(u, v) gegeben. 


2.1.1.2 Kantendetektion mittels erster Ableitung 


Eine Kante in einem Bild entspricht einer starken Grauwertänderung ent- 
lang einer Richtung, wobei die angrenzenden Gebiete jeweils möglichst 
homogen sein sollten. Die Grauwertfunktion weist also eine starke Än- 
derung auf, weshalb sich Ableitungsoperatoren im Grundsatz sehr gut 
nutzen lassen, um eine Kante zu detektieren. Die erste Ableitung weist ein 
Maximum auf, die zweite Ableitung einen Nulldurchgang. 

Dabei sollte ein guter Kantendetektor natürlich jede Kante in einem 
Bild erkennen und möglichst wenige falsch-positive Ergebnisse haben. 
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Außerdem sollte die Kante möglichst genau im Bild lokalisiert werden 
können [13]. 

Eine Ableitung im diskreten Raum wird statt Differenzieren durch ein- 
fache Differenzen gebildet. Aus den beiden Ableitungen in Richtung u und 
v lässt sich dann der Gradientenbetrag errechnen: 


en, Gog 


Wao) у (2520 


Dieser ist ein Maß für die Stärke der Grauwertänderung [7]. Mit einem 
Schwellwert können dann Kanten bestimmt werden. In verrauschten Bil- 
dern führt diese Art der Kantendetektion aber zu vielen Fehlerkennungen. 
Ein anderes gebräuchliches und deutlich weniger störungsanfälliges Ver- 
fahren ist der Einsatz eines differenzierten Gauß-Tiefpass-Filters [7]. Der 
Gauß-Tiefpass 


E | 


1 
heaug(U, v) = Sage exp (- 902 


mit о? als Varianz, wird differenziert und zur Approximation von Motu, v) 
eingesetzt: 
Уд(и, v) = Уйсаив(и, v) * gu, v). 


Diese Faltung entspricht einer Glättung mit anschließender Ableitung [7]. 
Die oben bereits erwähnten Kriterien aus [13], nämlich geringe Fehlerken- 
nungen und hohe Lokalisierungsgenauigkeit, lassen sich auch mathema- 
tisch formalisieren. Eine numerische Optimierung der Kriterien resultiert 
im sogenannten Canny-Filter, für das der differenzierte Gauß-Tiefpass eine 
gute diskrete Näherung darstellt [13]. 


2.1.1.3 Kantendetektion mittels zweiter Ableitung 


Die zweite Ableitung des Bildes lässt sich recht einfach durch ein Laplace- 
Filter approximieren. Dieses nimmt keine Glättung vor und ist entspre- 
chend anfällig für Rauschen [7]. Ein Weg, um Rauschen zu unterdrücken, 
ist die Vorverarbeitung des Bildes mit einem Gauß-Tiefpass und anschlie- 
ßender Anwendung des Laplace-Filters [60]. Dies wird als Laplacian-of- 
Gaussian-Operator bezeichnet. Zur diskreten Implementierung kann das 
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Difference-of-Gaussian (DoG) Verfahren genutzt werden. Dafür wird ausge- 
nutzt, dass sich der Laplacian-of-Gaussian-Operator als Differenz von zwei 
Gauß-Funktionen mit verschiedenen Standardabweichungen annähern 
lässt [7]. 

Ein großer Nachteil von Verfahren zur Kantendetektion mittels der zwei- 
ten Ableitung ist, dass ein Nulldurchgang der zweiten Ableitung eine not- 
wendige, aber keine hinreichende Bedingung für eine Kante darstellt und 
so Fehler in der Kantendetektion passieren können. Der wesentliche Vorteil 
im Vergleich zu Verfahren, die die erste Ableitung nutzen, ist die hohe Lo- 
kalisierungsgenauigkeit, die aus der eindeutigen Position der Nullstellen 
resultiert [7]. 


2.1.1.4 Eckendetektion durch den Harris Corner Detector 


Die Grundidee des Harris Corner Detectors ist es, ein Fenster über das Bild 
zu legen und ein kleines Stück zu bewegen. Ob eine Ecke vorliegt, wird 
daran entschieden, wie der Intensitätsunterschied zwischen dem origina- 
len Fenster und dem verschobenen Fenster genau aussieht. Gibt es keine 
großen Änderungen, ist weder eine Kante noch eine Ecke vorhanden, gibt 
es große Änderungen in einer Richtung, liegt eine Kante vor, nämlich recht- 
winklig zur Richtung der größten Änderung. Liegen große Änderungen 
in allen Richtungen vor, handelt es sich um eine Ecke [33]. 

Diese qualitative Darstellung lässt sich wie folgt formalisieren. Der In- 
tensitätsunterschied Е, lässt sich abhängig von der Verschiebung um (=, у) 
und dem betrachteten Fenster w(u, v) darstellen als: 


E,(z,y) A. w(u,v) (glu + то +y) — д(и,0))?. 


(uv) 


Diese Funktion ist allerdings aufwendig zu berechnen für alle Bildpunk- 
te, weshalb sie durch Taylor-Expansion angenähert wird. In Matrixschreib- 
weise ergibt sich dann: 


T mi 9, Iu 
Еау) = (e.y)M(e,y)" mit = 27 


wobei g,, und g,, die partiellen Ableitungen in die entsprechende Richtung 
sind. Die Eigenwerte von М seien Л, und А». 
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Im Sinne des oben vorgestellten Ansatzes müssen jetzt beide Eigenwerte 
groß sein, damit eine Ecke vorliegt. Wenn genau ein Eigenwert klein und 
der andere groß ist, liegt eine Kante vor. Wenn beide Eigenwerte klein sind, 
hat das betrachtete Fenster keine Ecke oder Kante und folglich eine eher 
homogene Grauwertverteilung. 

Ein Gütemaß R, dass diese Zusammenhänge ausdrückt, wird in [33] 
vorgeschlagen: 

В = АА — k(Ay + А,)?, 


wobei k eine Konstante ist. Das Gütemaß R ist dabei für Kanten negativ 
und für Ecken positiv. Wenn R unter einem bestimmten Schwellenwert 
liegt, handelt es sich um ein Fenster ohne Kanten und Ecken. 


2.1.1.5 Eckendetektion durch den Shi-Tomasi Detector 


Zwischen dem Harris Corner Detector and dem Shi-Tomasi Detector [73] 
gibt es große Parallelen [84]. Einzig das Gütemaß wird hier anders gewählt, 
nämlich: 

R=min{A,, А}. 


Hierbei muss das Gütemaß dann größer als ein vorher festgelegter Schwel- 
lenwert sein, damit der untersuchte Punkt als Ecke klassifiziert wird. Das 
Ausgangsproblem des Shi-Tomasi Detectors ist dabei das Auffinden von 
Ecken, die bei kleinen Bewegungen in einer Bilderabfolge benutzt werden 
können, um die Bewegung der Kamera durch die Szene festzustellen. Eine 
möglichst große Anzahl an Ecken zu finden ist hier also nicht das Ziel, son- 
dern nur solche, die mit hoher Wahrscheinlichkeit auch in einem späteren 
Bild noch als Ecke zu identifizieren sind [73]. 


2.1.2 Transformationen und Projektionen von digitalen 
Bildern 


2.1.2.1 Fourier-Transformation von Bildern 


Die Fourier-Transformation wird im Allgemeinen eingesetzt um ein Si- 
gnal g(x) in die komplexe Funktion G( f) zu transformieren, die von der 
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Frequenz f abhängig ist. Dabei wird das Signal in die beteiligten harmoni- 
schen Funktionen unterschiedlicher Frequenzen zerlegt [7]. Die Fourier- 
Transformation ist durch 


Fo) = Gs) = | "here 


gegeben. 

Eine sehr wichtige Eigenschaft der Fourier-Iransformation für die An- 
wendung in der Bildverarbeitung ist die Umkehrbarkeit. Die inverse Fou- 
rier-Iransformation ist als 


Fan) = | WEE 


definiert. 

Die Grauwertfunktion g(u, v) eines Bildes kann mittels der zweidimen- 
sionalen diskreten Fourier-Transformation (2D-DFT) in die beteiligten Fre- 
quenzen zerlegt werden. Die komplexe Funktion G(f,, fe) weist dann 
ebenfalls Frequenzen für die beiden Dimensionen auf. Bei vektorieller 
Notation mit w = (u,v)? und f = (fu, LIT ist die Definition der 2D-DFT 


F(g(w)) = G(f) = Г g(w) exp (-j: 27: fTw) dw. 


Die Transformation des Bildes aus dem Ortsraum in den Frequenzraum 
kann für die Bildverarbeitung hilfreich sein, da manche Operationen im 
Frequenzraum einfacher durchzuführen sind. Dies gilt 2. В. für die Unter- 
drückung von frequenzabhängigen Artefakten im Bild oder zur einfache- 
ren Auslegung von Filtern [28]. Mittels der inversen Fourier-Transformation 
kann das Bild danach in den Ortsraum zurück transformiert werden. 


2.1.2.2 Parallelprojektion 


Die Parallelprojektion ist eine Projektion, bei der das Projektionszentrum 
im Unendlichen liegt und Punkte im Raum entsprechend auf parallelen 
Projektionsstrahlen auf einen anderen Raum abgebildet werden [7]. Oft- 
mals werden Punkte aus einem dreidimensionalen Raum auf eine Ebene 
abgebildet. Die Projektionsstrahlen treffen dabei im Allgemeinen nicht 
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rechtwinkelig auf die Ebene. Falls die Projektionsstrahlen rechtwinklig auf 
die Projektionsebene fallen, wird auch von einer Orthogonalprojektion 
gesprochen. Bei der Verarbeitung zweidimensionaler Bilder entspricht 
eine Orthogonalprojektion entlang der Achsen des Bildes der Summe der 
einzelnen Pixel in einer Zeile bzw. Spalte [7]. 


2.2 Umgebungsbedingungen und 
Voruberlegungen 


2.2.1 Umgebungsbedingungen von 
Laserflachbettmaschinen 


Es gibt viele Publikationen zu visuellen Methoden zur Erfassung der Form 
eines Objektes. Welche Methode eingesetzt wird, hängt auch davon ab, 
welche für die Aufnahme relevanten Umgebungsbedingungen vorherr- 
schen. Für die Erkennung der Auflagestege bedeutet dies, dass das Mess- 
system mit sich ändernden Lichteinflüssen, Vibrationen der Maschine 
und unterschiedlichen optischen Eigenschaften der Auflagestege umgehen 
können muss. Letztere sind z. В. durch die unterschiedliche Reflektivität 
der Auflagestege verschiedener Materialien oder bei verschiedenen Ver- 
schlackungsgraden gegeben. Außerdem müssen physische Änderungen 
der Auflagestege, also verbogene Spitzen, Verschlackungen, verschweißte 
Butzen oder inkorrekt eingelegte Auflagestege, richtig erkannt werden. 

Die Lichtverhältnisse außerhalb des Maschinenkörpers können nicht 
beeinflusst werden, da eine Abdunkelung der gesamten Maschine in ei- 
ner Blechfertigung mit dem An- und Abtransportieren von Blechen nicht 
praktikabel ist. Dies ist ein Problem, da die möglichen Beleuchtungsszena- 
rien von einer kaum beleuchteten Halle während der Nachtschicht bis zu 
direktem Sonnenlicht auf der Palette (siehe Abbildung 2.1) reichen. 

Auch die möglichen Positionen für Sensoren sind eingeschränkt. Da 
Blech von den drei offenen Seiten des Palettenwechslers und von oben 
eingelegt werden können muss, bleiben außen nur Positionen direkt an der 
Maschineneinhausung übrig. Im Maschinenkörper muss berücksichtigt 
werden, dass der Schneidkopf ungestört verfahren kann, was ebenfalls die 
möglichen Positionen stark einschränkt. 
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Abbildung 2.1 Sonnenlicht auf der Palette einer LFM. 


Ein weiterer Faktor ist der Einfluss des Messsystems auf den Produk- 
tionsablauf. Vor allem bedeutet dies, dass die für die Messung benötigte 
Zeit den Maschinenzyklus als Ganzes möglichst wenig verlängert. Die Ma- 
schinenzeiten werden dabei gewöhnlicherweise in Haupt- und Nebenzeit 
eingeteilt, wobei die Hauptzeit das wirkliche Bearbeiten des Werkstücks 
bezeichnet. Die Messung der Stege wäre also ein Teil der Nebenzeit. 

Darüber hinaus lassen sich die Auflagestege unterschiedlich genau ver- 
messen. Da Stege nur in bestimmten, an der Palette bekannten Stegsteck- 
positionen eingelegt sein können, reicht es für eine erste Schätzung der 
Positionen der Auflagespitzen aus, das Vorhandensein eines Steges fest- 
zustellen, was im Folgenden Stegdetektion genannt wird. Ein deutlich 
genaueres Vorgehen ist selbstverständlich, jede einzelne Auflagespitze zu 
vermessen, was als Spitzendetektion bezeichnet wird. 
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2.2.2 Visuelle Detektionsverfahren 


Gegeben die Umgebungsbedingungen in der Blechfertigung, ergeben sich 
auf den ersten Blick mehrere visuelle Messverfahren als geeignete Kandi- 
daten. Hier sollen kurz vier Messverfahren dargestellt werden. 

Die am einfachsten umzusetzende Methode ist die Aufnahme eines 
einzelnen Bildes und anschließende Detektion von Objekten durch die 
Bildverarbeitungsalgorithmen. Dieser Ansatz wird insbesondere dann 
gewählt, wenn das zu erkennende Objekt vorher bekannt ist. Mögliche 
Algorithmen sind z.B. Texturanalyse, Ecken- und Kantendetektion oder 
Fourieranalyse. Die so generierten Features können mit klassifizierenden 
Verfahren, z. В. dem Matched Filter zur Stegerkennung genutzt werden [7, 
46]. Verfahren mit kontrollierter oder sogar aktiver Beleuchtung sollten 
besser funktionieren als Verfahren ohne kontrollierte Belichtungsumge- 
bung. Der große Vorteil an der Aufnahme einzelner Bilder ist die geringe 
zusätzliche Nebenzeit, die sich auf die Zeit zur Aufnahme eines Bildes 
beschränkt, also i. а. К. weniger als eine Sekunde. Dies ist nur möglich, 
wenn die Kamera außen an der Maschine befestigt ist, da es im Maschinen- 
körper keine Position gibt, von der aus alle Stege zu sehen sind. Auf dem 
Palettenwechsler kann nach dem Absortieren der geschnittenen Teile und 
des Restgitters die leere Palette mit allen Stegen aufgenommen werden, 
bevor ein neues Rohblech auf die Palette gelegt wird. 

Eine sehr weit verbreitete Methode zur dreidimensionalen Formerfas- 
sung ist die Laser-Triangulation. Dabei wird ein Laserpunkt oder eine 
Laserlinie auf das Objekt projiziert. Die Position des Lasers wird von einem 
Kamerasensor detektiert, der in einem festen Winkel und Abstand zum 
Laser steht. Von der Position, bei welcher der Laser auf dem Kamerasen- 
sor wahrgenommen wird, lässt sich auf die Form des Objekts schließen. 
Dabei sind Genauigkeiten von deutlich unter einem Zehntel Millimeter 
möglich [9, 18]. Größter Nachteil dieser Methode ist die relativ lange zu- 
sätzliche Nebenzeit, da eine Laser-Triangulation immer ein scannendes 
Verfahren ist. Für eine LFM bedeutet dies, das die Palette ungefähr 20 
Sekunden leer sein muss, falls der Antrieb des Palettenwechslers zum 
Scannen einer Laserlinie genutzt wird. Außerhalb der Maschine wäre auch 
eine Bewegung der Laserlinie über die stehende Palette denkbar, dort sind 
aber Fremdlichteinflüsse stärker vorhanden. 
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Ein weiterer Ansatz zur Erkennung der dreidimensionalen Form eines 
Objekts ist Stereo Vision. Der Ansatz folgt dem menschlichen Sehen, bei 
dem der perspektivische Unterschied zwischen den zwei Augen genutzt 
wird, um Informationen über die dreidimensionale Gestalt eines Objekts 
zu extrahieren [61]. Aus dem Winkel, unter dem derselbe Punkt eines be- 
trachteten Objekts von den beiden Kameras gesehen wird, kann die Distanz 
dieses Punktes berechnet werden. Das größte Problem hierbei ist, dass 
derselbe Punkt in den beiden leicht unterschiedlichen Bildern eindeutig 
identifiziert werden muss, was oft als Korrespondenzproblem bezeichnet 
wird [7, 65]. Daraus entsteht zur Vermessung der Auflagestege von LFMn 
auch der größte Nachteil: bei periodischen Strukturen oder wenn es zu 
Reflexionen kommt, ist das Korrespondenzproblem sehr schwer oder gar 
nicht zu lösen [42, 52]. Da Auflagestege ausschließlich aus einer periodi- 
schen Struktur, nämlich den Spitzen, und aus metallischem Material, das 
Reflexionen begünstigt, bestehen, erscheint der Stereo Vision Ansatz nur 
bedingt geeignet. Dies konnte in einem Vorversuch bestätigt werden. 

In den letzten Jahren hat das Time-of-flight Messprinzip deutlich an Ein- 
fluss gewonnen, insbesondere im Bereich des autonomen Fahrens [21, 53, 
56]. Wie schon der Name sagt, wird von einer Lichtquelle aus ein Licht- 
bündel auf das Objekt gestrahlt und in einem Sensor detektiert, wann das 
reflektierte Licht wieder beim Sensor eintrifft. Aus der zeitlichen Differenz 
und der Lichtgeschwindigkeit lässt sich die Distanz vom Sensor zum Ob- 
jekt berechnen [52, 65]. Der größte Nachteil dieses Ansatzes ist die geringe 
Messgenauigkeit. Für Messabstände von 7,5 m bis 10 m liegt die Genauig- 
keit im einstelligen Zentimeterbereich [31, 52]. Damit ist die Genauigkeit im 
Vergleich zur Größe der Auflagestege zu gering und um Größenordnungen 
schlechter als 2. В. bei einer Messung mittels Laser-Triangulation. 

Da die beiden zuletzt vorgestellten Ansätze keine ausreichende Ge- 
nauigkeit versprechen, werden die davor vorgestellten Ansätze vertieft. 
Zunächst wird die Umsetzung eines Messaufbaus zur Laser-Triangulation 
vorgestellt. Danach wird erst die Objektdetektion auf einfachen Grauwert- 
bildern dargestellt, um dann ein darauf basierendes Verfahren mit aktiver 
NIR-Beleuchtung zu betrachten. 
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2.3  Stegdetektion mittels Laser-Triangulation 


Das Prinzip, eine Distanz durch Berechnungen in Dreiecken zu messen, 
ist schon lange bekannt. Es beruht darauf, dass ein Dreieck vollständig 
definiert ist, wenn zwei Winkel und die Länge einer Seite bekannt sind. 
Beim Messen sind i.d.R. ein Winkel und die Länge einer Seite bekannt 
und der zweite Winkel wird gemessen. 

Die Ausführungen dieses Unterkapitels sind in ähnlicher Form schon in 
[89] veröffentlicht worden. 


2.3.1 Stand der Technik 


Bei einem Laser-Triangulationssystem (LTS) wird ein Laserstrahl als Punkt 
oder Linie auf die zu messende Oberfläche projiziert. Die Position des 
Lasers auf der Oberfläche wird mit einer Kamera gemessen. Bei einem 
üblichen LTS-Aufbau wird der Laser direkt über der zu messenden Ober- 
fläche aufgestellt, sodass der Laserstrahl rechtwinklig zur Oberfläche steht. 
Die Kamera ist so angebracht, dass ihre optische Achse in einem Winkel 
zum Laserstrahl steht (siehe Abbildung 2.2 links). Die Bilder, die bei der 
Vermessung einer Oberfläche entstehen, zeigen also immer das Höhen- 
profil an einer bestimmten x-Position. Die Auslenkung der Laserlinie auf 
dem Kamerasensor hängt direkt mit der Höhe des Messobjekts zusam- 
men. Die Oberfläche wird dann nach und nach gescannt, was mit einem 
punktförmigen Laser natürlich viel länger dauert als mit einem Linienlaser. 

Ein leicht veränderter Aufbau wurde in einer Studie benutzt, bei der es 
wichtig war, neben dem Höhenprofil auch ein Farbbild des Messobjekts 
zu erhalten. Deshalb wurde die Anordnung der Kamera und des Lasers 
vertauscht (siehe Abbildung 2.2 Mitte). So konnten zusätzliche Transfor- 
mationen bei der Zuordnung von Bereichen des Höhenbildes zu Bereichen 
des Farbbildes vermieden werden [49]. 

Ein Problem, das bei der Auswahl jedes LTS-Aufbaus beachtet werden 
muss, ist Abschattung. Da die meisten Objekte nicht aus einer glatten Ober- 
fläche bestehen, kann es sein, dass entweder Bereiche des Messobjekts nicht 
von der Position der Kamera aus zu sehen sind, da sie verdeckt werden, 
oder dass der Laser durch andere Bereiche des Messobjekts verdeckt wird. 
In abgeschatteten Bereichen kann also nicht gemessen werden. Der Aufbau 
eines LTS muss folglich ggf. auf das Messobjekt angepasst werden [7]. 
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Abbildung 2.2 Verschiedene Aufbauten für die Laser Triangulation. Links: Laser rechtwink- 
lig zum Objekt. Mitte: Kamera rechtwinklig zum Objekt. Rechts: Aufbau für diese Arbeit. 


2.3.1.1 Kalibrierung 


Die Kalibrierung eines LTSs kann in zwei Schritte unterteilt werden. Die 
intrinsische Kalibrierung sorgt dafür, dass mögliche Verzerrungen im op- 
tischen System der Kamera korrigiert werden. Durch die extrinsische Kali- 
brierung wird die Transformation vom Kamerakoordinatensystem zum 
Weltkoordinatensystem hergestellt [37, 72]. 

Eine weit verbreitete Methode zur intrinsischen Kalibrierung ist das 
Lochkamera-Modell. In diesem Modell wird eine perfekte Lochkamera 
angenommen und die auftretenden Verzerrungen werden anhand der Eck- 
punkte eines Schachbrettmusters berechnet. Dabei wird die sogenannte 
Kameramatrix, welche aus den Brennweiten des Objektivs und den Koor- 
dinaten des Pixels, der sich auf der Achse des optischen Systems befindet, 
besteht, berechnet. 

In der extrinsischen Kalibrierung wird die korrekte Transformation zwi- 
schen Kamerakoordinatensystem und Weltkoordinatensystem berechnet. 
Dies wird i. d. R. dadurch erreicht, dass ein Kalibrierobjekt mit sehr genau 
bekannten Abmessungen im LTS vermessen wird [78]. Wenn man aus- 
schließlich an der Höhe des Objekts interessiert ist, ist dies gleichbedeutend 
mit einer passenden Skalierung der Höhenkoordinate (z. B. in [49]). Die 
extrinsische Kalibrierung kann eine Fehlerquelle für die Messung sein, 
wenn das Kalibrierobjekt und das zu messende Objekt nicht sehr ähnliche 
Oberflächeneigenschaften haben [64]. 
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Die Kalibrierung kann auch in explizite und implizite Methoden aufge- 
teilt werden. Explizite Methoden versuchen, physikalische Parameter zu 
bestimmen, wie z.B. die Brennweiten des Objektivs oder die Koordinaten 
des optischen Zentrums auf dem Kamerasensor. Implizite Methoden hin- 
gegen haben ausreichend Parameter, um die möglichen Fehler auszuglei- 
chen, diese lassen aber keinen direkten Rückschluss auf die physikalischen 
Größen des Messsystems zu [37]. 


2.3.1.2 Detektion der Laserlinie 


Die Messungenauigkeit eines LTS kann verschiedene Gründe haben, die 
folgenden Kategorien zuzuordnen sind [7]: 


e Projektionsqualität der Kameraoptik, 
= Projektionsqualität des Lasers, 

= Auflösung des Kamerasensors, 

= Speckle Noise. 


Der letzte dieser Punkte ist im Detail erforscht worden, da es sich dabei 
um den physikalisch limitierenden Faktor des Laser-Iriangulationsprin- 
zips handelt [18, 34]. Die Projektionsqualität der Kameraoptik und die 
Auflösung des Kamerasensors sind durch die Auswahl der Hardware 
gegeben. Eine geringe Projektionsqualität des Lasers, also eine relativ breite 
Linie auf dem Messobjekt, kann bis zu einem gewissen Grad abgemildert 
werden. Wenn die Laserlinie mehrere Pixel auf dem Sensor belichtet, lässt 
sich die Position der Laserlinie mit Subpixelgenauigkeit berechnen [66]. 
Dafür sind in der Literatur mehrere Verfahren bekannt, von denen hier 
drei vorgestellt werden sollen. 

Das Massenschwerpunktsverfahren (CoMP) errechnet für jede Spalte 
des Bildes den Massenschwerpunkt der Grauwerte [45], wobei f(r) den 
Grauwert in der Zeile r bezeichnet. Wenn nur ein klares Maximum in den 
Daten vorhanden ist, kann das Verfahren auf die gesamte Spalte ange- 
wendet werden. Wenn es aber unerwünschte weitere Maxima gibt, kann 
die Berechnung auch auf einen Bereich zwischen zwei Zeilen [k,l], mit 
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k < l eingeschränkt werden. Die gesuchte subpixelgenaue Position der 
Laserlinie Das, kann dann wie folgt berechnet werden: 


. see 
У 


Da eine Laserlinie als erste Annäherung eine gaußverteilte Intensität hat, 
ist es ein weiteres, recht intuitives Verfahren, das Maximum der am besten 
zu den Grauwerten der Spalte passenden Gauß-Verteilung als Position der 
Laserlinie zu nehmen. Die Gauß-Verteilung h(a) ist dabei gegeben durch: 


h(x) = exp ES . 


20? 


Das Maximum pmax wird bestimmt, indem die Parameter der Gauß- 
Verteilung variiert und mit den Grauwerten gefaltet werden: 


Pmax = argmax{ f(r) ж h(r)}. 


Das Maximum der am besten passenden Gauß-Verteilung wird dann als 
Position der Laserlinie betrachtet. Dadurch ist die Bestimmung von Pmax 
mit dem Verfahren der Annäherung einer Gauß-Verteilung (GaussP) im 
Vergleich zu den vorher vorgestellten Verfahren sehr rechenintensiv. 

Ein recht bekanntes weiteres Verfahren ist der nach seinen Erfindern 
benannte Blais-Rioux-Filter [8]. Dazu wird die Ableitung des Signals gebil- 
det und Rauschen durch einen Finite Impulse Response Filter unterdrückt. 
Der große Vorteil dieses Verfahrens ist die schnelle Berechenbarkeit, was 
für den ursprünglichen Anwendungsfall von großer Wichtigkeit ist. In 
der Literatur existieren natürlich noch viele weitere Verfahren, die aber 
eine untergeordnete Rolle spielen und hier nicht weiter betrachtet werden 
sollen. 


2.3.2 Aufbau und Kalibrierung 


Es gibt mehrere mögliche Anordnungen zum Aufbau eines LTS zur Ver- 
messung der Auflagestege einer LFM. Um allen Einflüssen einer realen 
Produktion ausgesetzt zu sein, wurde entschieden, den Aufbau in einer 
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LFM zu realisieren und nicht ein kleines Stück der Palette in einem Labor 
für die Tests zu benutzen. 

Der gewöhnliche Aufbau eines LTS ist in Abbildung 2.2 links dargestellt. 
Dieser Aufbau führt wie schon angesprochen zu Abschattungsproblemen 
bei den äußeren Spitzen. Außerdem enthalten nur sehr wenige Bilder einer 
Messung tatsächlich nützliche Informationen, schließlich ist zwischen zwei 
Stegen, die 3 mm breit sind, ungefähr 30 mm Abstand. Dadurch wird bei ei- 
ner gegebenen Geschwindigkeit des Palettenwechslers die Anforderung an 
die Bildwiederholfrequenz sehr hoch. Ein besserer Aufbau für diese Mess- 
aufgabe ist also so ausgerichtet, dass die Seitenfläche der Stege vermessen 
wird. Dadurch kann sowohl das Abschattungsproblem umgangen als auch 
mit jedem Bild nützliche Information gewonnen werden. Ein Nebenef- 
fekt dieses Aufbaus ist, dass die Schlacke an den Stegen besser vermessen 
werden kann. Der Aufbau ist in Abbildung 2.2 rechts dargestellt. 

In einer LFM kommt hinzu, dass der Schneidkopf die Maschine praktisch 
bis unters Dach ausfüllt, weil zwischen den Kabeln und dem Dach nur 
wenige Zentimeter Platz sind. Die Kamera und der Linienlaser müssen 
also oben an der Seite der Maschine angebracht werden (siehe Abbildung 
2.3). Eine Messung von unten ist zwar zur Detektion von Stegen möglich, 
die Spitzen werden dort aber nicht zu erkennen sein. Außerdem könnten 
die Kamera und der Laser durch Schlackespritzer beschädigt werden. 

Die Kamera wurde mit dem in Kapitel 2.3.1.1 vorgestellten Schach- 
brettmuster-Verfahren intrinsisch kalibriert. Die extrinsische Kalibrierung 
des Messsystems stellt eine Besonderheit dar, da sie nicht mit einem zu- 
sätzlichen Kalibrierobjekt, sondern mit neuen Stegen durchgeführt wird. 
Dies ist deshalb möglich, weil die neuen Stege mit einer Toleranz von 
+ 0,05mm auf der Maschine selbst geschnitten werden. Dieses Vorge- 
hen bringt zwei Vorteile für die industrielle Anwendung. Erstens ist es 
so möglich falls notwendig laufend neu zu kalibrieren und zweitens ist 
kein spezielles Kalibrierobjekt nötig, das ggf. verloren gehen könnte. Eine 
neue Kalibrierung wird notwendig, wenn der Maschinenkörper verzogen 
ist, z.B. wenn das Fundament der Halle sich setzt. Außerdem können die 
Stege aus unterschiedlichen Materialien und Materialqualitäten sein, was 
deutlichen Einfluss auf ihre optischen Oberflächeneigenschaften hat. Die 
Oberflächeneigenschaften des Messobjekts können als Fehlerquelle aus- 
geschlossen werden, wenn die Stege als Kalibrierobjekt genutzt werden 
[64]. 
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Abbildung 2.3 Der Laser-Triangulationsaufbau im Innenraum der LFM. 


Zur extrinsischen Kalibrierung soll aus den Aufnahmen die Position der 
Senken und Spitzen bestimmt werden. Von diesen ist durch die bekann- 
te Geometrie der Stege auch die tatsächliche Position bekannt. Aus der 
Abweichung können dann die Parameter der extrinsischen Kalibrierung 
errechnet werden. 

Um die Position der Senken und Spitzen zu bestimmen, werden die 
Bilder zunächst mit einem sehr geringen Schwellenwert binarisiert. Dies 
unterdrückt auch in den Aufnahmen vorhandenes Dunkelrauschen des 
Kamerasensors. Die Bilder, die zu einem Steg gehören, können eindeutig 
über einen Schwellenwert identifiziert werden, weil auf ihnen viele Pixel 
hohe Grauwerte haben. Zwischen zwei Stegen sind immer Aufnahmen, 
auf denen alle Pixel sehr geringe Grauwerte haben. 

Die Bilder eines Steges zeigen zunächst eine durchgehende Linie, die 
an den Senken aufbricht und gestrichelt ist, bis an den Spitzen nur noch 
weiße Punkte in einer Linie vorhanden sind (siehe Abbildung 2.4). Die 
Addition dieser binarisierten Bilder resultiert in einem Bild des gesam- 
ten Steges, wie in Abbildung 2.5 oben dargestellt. Anhand dieser Bilder 
kann die Kontur der Stege mittels Border Following [76] festgestellt werden. 
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Abbildung 2.4 Drei Bilder der Aufnahmesequenz eines Steges. Oben: Der untere, vollstän- 
dige Teil des Steges. Mitte: Unterbrochenes Linienmuster in der Mitte der Auflagespitzen. 
Unten: Stark ausgedünntes Muster in der Nähe der Spitzen. 


Border Following bedeutet dabei lediglich, dass die Pixelkoordinaten der 
Grenze zwischen dem schwarzen und weißen Gebiet ermittelt werden. 
Die Position der Senken und Stege kann so noch nicht einfach extrahiert 
werden, da die Kontur in dieser Form zu viel Rauschen enthält. Dieses 
kann mit dem Ramer-Douglas-Peucker Algorithmus entfernt werden [19, 
68]. Dieser Algorithmus prüft für jeden Punkt eines Polygons, ob der Ab- 
stand zwischen dem betrachteten Punkt und der direkten Verbindung des 
Punktes davor und danach einen gegebenen Schwellenwert unterschreitet. 
Wenn dies der Fall ist, wird der Punkt entfernt und das Polygon verliert 
eine Ecke. Sonst wird der Punkt beibehalten und der nächste Punkt des 
Polygons wird geprüft. Dadurch kann das Rauschen effektiv unterdrückt 
werden. 

Als Verfahren zur Detektion von Ecken kommen die in Kapitel 2.1.1 
vorgestellten Harris Corner Detector und der Shi-Tomasi Detector in Frage. 
Letzterer zeigte dabei deutlich bessere Ergebnisse für diese Aufgabe, weil 
er immer alle wahren Senken und Spitzen finden konnte und gleichzei- 
tig sehr wenige Ecken fälschlicherweise als solche klassifizierte. Deshalb 
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sind alle weiteren Berechnungen mit dem Shi-Iomasi Detector durchge- 
führt worden. Die Ergebnisse werden im folgenden Schritt auf Plausibilität 
geprüft, wobei vor allem die Distanzen zur nächsten Senke bzw. Spitze 
und die Höhe der erkannten Ecke genutzt werden. Das Ergebnis ist in 
Abbildung 2.5 unten dargestellt. 


Abbildung 2.5 Oben: Die aus mehreren Bildern zusammengesetzte Form eines Steges. 
Unten: Schätzung der Kontur mit dem Shi-Tomasi-Detektor. 


Da die Koordinaten der Senken und Spitzen im Bild (u, v) jetzt bekannt 
sind, kann die eigentliche Kalibrierung mit zwei Gleichungen angenähert 
werden. Die erste ist eine perspektivische Transformation, die dazu führt, 
dass alle Spitzen so dargestellt werden, als hätte man sie direkt von vorne 
gesehen (siehe Formel 2.1). Das Pixel (и, v) des Bildes wird dabei auf das 
Pixel (из, vp) abgebildet, wobei n ein Skalierungsfaktor ist. 


N Une Ou 412 443] |U u 
пи | = |421 аә @а23| |v| =А |v (2.1) 
n WEN @32 @33 1 1 


Diese Punkte (из, 0, „) werden dann durch die Matrix В in Maschi- 
nenkoordinaten transformiert (siehe Formel 2.2). w ist dabei die Nummer 
des Bildes in der Aufnahmesequenz. 


u u 
Tact bn bus big бд Sen А 
Pact = |Yacı | = |621 652 0653 bau 2 =B 2 (2.2) 
Zact bs, b32 653 bau 1 1 
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Für jede Spitze und Senke ist durch die bekannte Geometrie der Stege ei- 
ne als wahr angenommene Position in Maschinenkoordinaten p,, bekannt. 
Diese kann von der tatsächlichen Position leicht abweichen, da die Stege 
beim Einlegen in die Palette einer leichten Torsion ausgesetzt sind. Die 
Abweichungen bewegen sich im Bereich von ungefähr einem Millimeter. 

Die Matrizen A und B werden nun durch Minimierung der euklidischen 
Distanzen zwischen allen Punkten ру, und Pact approximiert: 


e= У || Pyt — Pact ||- 
Pact 
Dabei wird ein kombiniertes Verfahren aus Basin Hopping [79] für die 
globale Optimierung und Sequential Least Squares [50] für die lokale Opti- 
mierung eingesetzt. Die Implementierung stammt aus dem Software-Paket 
SciPy. Das Verfahren wurde für 15 Iterationen ausgeführt, wobei das Stop- 
Kriterium für die lokale Optimierung durch = < 1076 gegeben war. 


2.3.3 Bildverarbeitung und Transformationen 


Sobald die Kalibrierung berechnet worden ist, können die Stege vermes- 
sen werden, wozu allerdings noch ein paar Bildverarbeitungsschritte und 
räumliche Transformationen notwendig sind. 

Als erstes werden die Bilder auf den Bereich begrenzt, der relevante 
Informationen enthält, in diesem Fall von 90 x 2440 Pixel auf 45 x 1660 
Pixel, was einer Messbreite von ungefähr 755 mm entspricht. 

Da die Laserlinie auf vielen Bildern durch die Abstände zwischen zwei 
Spitzen unterbrochen ist, wird ein Schwellenwert für den maximalen Grau- 
wert einer Spalte gesetzt, der überschritten werden muss, damit die Spalte 
überhaupt weiter verarbeitet wird. Dann werden zwei von den in Kapitel 
2.3.1.2 vorgestellten Liniendetektions-Verfahren eingesetzt, um die Posi- 
tion der Laserlinie subpixelgenau festzustellen. Zum einen ist dies die 
CoMP Methode, die einmal nur auf die 11 Pixel rund um den maximalen 
Grauwert beschränkt ist (CoMP11) und einmal auf der gesamten Spalte 
(CoMP45). Zum anderen wird das GaussP Verfahren getestet, um anschlie- 
ßend diese drei Möglichkeiten vergleichen zu können. Die CoMP Methode 
wurde dabei gewählt, weil sie in einer anderen Studie besonders gut mit 
metallischen Reflexionen umgehen konnte [49]. Die GaussP Methode wird 
als der intuitive Ansatz ebenfalls getestet. Andere Ansätze sind vor allem 
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auf hohe Berechnungsgeschwindigkeit ausgelegt, was in diesem Fall aber 
nur eine untergeordnete Rolle im Vergleich zur Genauigkeit spielt, und 
sollen deshalb im Folgenden nicht weiter betrachtet werden. 

Das Ergebnis der Liniendetektionsverfahren ist eine Liste von Messpunk- 
ten р = |z, у, 2]. Dabei ist die x-Koordinate die Nummer des betrachteten 
Bildes, die y-Koordinate die Nummer der Bildspalte und die z-Koordinate 
das Ergebnis des Liniendetektionsverfahrens für die Bildspalte. Die Va- 
riablen x und у sind also ganzzahlig, während z eine Gleitkommazahl 
ist. 

Nun werden die während der Kalibrierung bestimmten Transformatio- 
nen zur Entzerrung auf den Messpunkten ausgeführt. Da die Kamera und 
der Laser nicht rechtwinklig zur gemessenen Oberfläche stehen und der 
Aufnahmezeitpunkt eines Bildes in der Kalibrierung nicht berücksichtigt 
worden ist, sind die Stege nicht aufrecht in Richtung der z-Achse, son- 
dern wirken nach hinten gelehnt. Durch den ungewöhnlichen Aufbau des 
hier implementierten LTS müssen noch zwei räumliche Transformationen 
durchgeführt werden, die eine geometrische Scherung der Messpunkte 
bewirken. Zuerst sollen die Stege in Richtung der z-Achse ausgerichtet 
werden. Der dafür nötige Winkel y lässt sich leicht aus der Geschwindigkeit 
des Palettenwechslers vp, der Höhe einer Spitze h, dem Abstand zwischen 
der Senke und der Spitze eines Steges in Millimetern d, der Anzahl an 
Bildern 1 und der Bildrate der Kamera f errechnen: 


ERTL (eg Pd 
y = arctan а ‚mi У ?. 


Die Rotation der Punkte р zu p,.,, ist dann durch: 


Lot 1 0 0 x 
Yrot | = 0 cos E sin ү H 
Zrot 0 siny cosy 2 


gegeben. Die einfache Rotation führt aber dazu, dass die Stege entlang der 
x-Achse nicht mehr den richtigen Abstand haben, später aufgenommene 
Stege höhere Werte in der z-Dimension haben und in der z-Dimension 
gestaucht werden. Diese Probleme können durch eine weitere Transfor- 
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Abbildung 2.6 Farblich codierte Höhendarstellung der gesamten Palette, bestehend aus 
mehreren Millionen Messpunkten. 


mation behoben werden, deren Ergebnis die korrekten Messpunkte ри, 
sind: 


"Frot 
Im sin y 
Pm = |Ym | = Yrot cos у ` Zrot 
Zm Žrot 
cos 7 


Das Ergebnis dieser Bearbeitungsschritte ist eine Punktewolke der ein- 
zelnen Messpunkte. Ein Beispiel mit ungefähr zwei Millionen Messpunkten 
ist in Abbildung 2.6 dargestellt. 


2.3.4 Methode zur Unterdrückung von Reflexionen 


Wie schon in Abbildung 2.4 zu sehen ist, sind die Aufnahmen mit Ausnah- 
me der Laserlinie fast vollständig dunkel. Allerdings gibt es an zwei Stellen 
entlang der Palette immer wieder helle Punkte über der eigentlichen La- 
serlinie. Diese treten dort auf, wo die Palette zwei Längsträger hat, die die 
Auflagestege jeweils nach einem Drittel der Palette unterstützen (vgl. auch 
Abbildung 2.15). Es wird also Licht von den Längsträgern zur Kamera 
reflektiert, obwohl diese nicht direkt vom Laser angestrahlt werden. 
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Abbildung 2.7 Ein Beispiel für ein durch den Längsträger ausgelöstes Artefakt. 


Wenn die Aufnahmen wie in den beiden vorhergehenden Unterkapiteln 
beschrieben bearbeitet werden, führt das zu Artefakten an einer bestimm- 
ten y-Position der Palette. Schließlich wird jeder detektierte Punkt so inter- 
pretiert, dass er auf der Ebene des Lasers liegen muss. Dies ist für das Licht 
vom Längsträger aber nicht der Fall. Die Artefakte sind zylinderförmig 
und in einem leichten Winkel zum Steg angeordnet (siehe Abbildung 2.7). 
Ungefähr 75 % der Stege auf den Aufnahmen zeigen diese Artefakte. Die 
Längsträger sind in der Punktwolke sonst nicht zu sehen, da sie zu tief 
liegen, um noch detektiert zu werden. 

Eine einfache Methode um dieses Problem zu beheben, ist, die Schwel- 
lenwerte, ab der eine Spalte überhaupt verarbeitet wird, anzuheben. Die 
Artefakte sind schließlich recht dunkel im Vergleich zur wirklichen La- 
serlinie. Dies führt aber zu unerwünschten Effekten an den Flanken und 
Spitzen der Auflagespitzen, die auch dunkler sein können als die mittle- 
ren Teile der Auflagespitzen. Da sonst kein Umgebungslicht im Maschi- 
nenkörper nachgewiesen werden konnte, erschien auch der Einsatz eines 
Schmalbandfilters mit Wellenlänge des Lasers keine geeignete Lösung. 
Das Licht auf dem Längsträger wird mehrfach reflektiert, wahrscheinlich 
von den Flanken der Auflagespitzen, die durch die Position des Lasers an 
der Maschinenseitenwand in einem dafür günstigen Winkel beleuchtet 
werden. 

Da die Artefakte strikt hinter den Auflagestegen auftreten und nicht 
darüber oder sich überkreuzen, kann die Anzahl an Messpunkten in einem 
Quader betrachtet werden, der durch die Punktewolke geschoben wird. 
Dieser erstreckt sich über die gesamte z-Achse, 30 Pixel in y-Richtung um 
die Position des Längsträgers und 5 Pixel in x-Richtung. Die Schrittweite 
bei der Verschiebung ist 1,5 Pixel in x-Richtung. Die Filterantwort ist die 
Anzahl an Messpunkten im Quader. Falls diese unter einem Schwellenwert, 
in diesem Fall 50, liegt, werden alle Punkte im Quader gelöscht. 
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Abbildung 2.8 Ein Ausschnitt der aufgenommenen Palette. Oben: Vor der Filterung der 
Artefakte. Unten: Nach Filterung der Artefakte. 


Auf den Auflagestegen aufgeschweißte Schneidbutzen haben zwar ähnli- 
che Charakteristika, dafür aber deutlich mehr Messpunkte als die maximal 
80 Messpunkte eines der zu filternden Artefakte. Bei der Anwendung auf 
den Aufnahmen wurden keine Messpunkte entfernt, die nicht zu den Ar- 
tefakten gehören. Gleichzeitig wurden 99 % der Messpunkte, die zu einem 
Artefakt gehören, entfernt. Dies ist auch eindeutig in Abbildung 2.8 zu 
erkennen. 


2.3.5 Umsetzung 


Ein entsprechender Aufbau wurde in einer LFM vom Typ TRUMPF Tru- 
Laser realisiert (siehe Abbildung 2.3). Während der Aufnahmen wurde 
die Arbeitsbeleuchtung ausgeschaltet, sodass der Maschineninnenraum 
dunkel war. Die verwendete Kamera war eine FLIR Blackfly S mit einem 
monochromen Sony IMX264 CMOS-Sensor. Die Bildhöhe wurde auf 90 
Pixelreihen begrenzt, um die Aufnahmegeschwindigkeit auf 150 Bilder pro 
Sekunde anzuheben. Jedes einzelne Bild hat deswegen eine Abmessung 
von 90 x 2448 Pixel. Die Kamera wurde mit einem Computar M0824-MPW2 
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Objektiv mit 8mm Brennweite kombiniert. Als Linienprojektor kam ein 
Z-LASER ZM18RF686 mit 640 nm Wellenlänge zum Einsatz. 

Der Messaufbau hatte einen Messabstand von ungefähr 880 mm und 
eine Messbreite von 755 mm, was ziemlich genau die Hälfte der Palette 
ist. Um die gesamte Palette zu scannen, wären also zwei Kameras oder 
ein passendes Weitwinkelobjektiv nötig. Die Geschwindigkeit der Palette 
während des Messvorgangs war ungefähr 70 mm/s. 

Jede Konfiguration an Stegen z.B. verschiedene Materialien, stark ver- 
schlackte Stege oder verbogene Spitzen wurde zwei Mal aufgenommen: 
Das erste Mal als tatsächliche Messung ohne die Arbeitsbeleuchtung der 
Maschine. Beim zweiten Mal wurde die Arbeitsbeleuchtung eingeschaltet, 
um bei unerwarteten Effekten in den Aufnahmen besser deren Ursprung 
feststellen zu können. 


2.3.6 Messgenauigkeit des Laser-Triangulationssystems 


Bei dem hier vorgestellten LTS haben die Kameraauflösung und die Auf- 
nahmegeschwindigkeit den größten Einfluss auf die Messgenauigkeit. Die 
Auflösung ist in allen drei Dimensionen ungefähr 0,45 mm/ Pixel. 

Da es sehr schwer ist, eine Grundwahrheit mit einer ausreichend ho- 
hen Auflösung aufzunehmen, wird die Validierung der Messgenauigkeit 
mit den Stegspitzen selbst durchgeführt. Deren Position ist, solange die 
Stege neu sind, aus der Steggeometrie bekannt. Dafür sind zwei Ansätze 
denkbar. Erstens kann die absolute Position der Spitzen über den gesamten 
Steg betrachtet werden. Zweitens kann der Abstand zwischen jeweils zwei 
Spitzen betrachtet werden. Letzteres ist deutlich weniger anfällig für eine 
mögliche Torsion des Steges in den Stegsteckplätzen, die zu einer Abwei- 
chung entlang der x-Achse führen kann, obwohl sich ein Steg natürlich an 
einer festen x-Position befindet. 

Um von den Daten der Punktewolke zu einer zweidimensionalen Pro- 
jektion der Stegspitzen zu kommen, wird der Bereich der z-Achse, in der 
sich die Stegspitzen befinden, durch Parallelprojektion auf die x-y-Ebene 
reduziert. Das Ergebnis dieser Projektion ist ein schwarzes Bild mit weißen 
Punkten an der Stelle des jeweiligen Messpunktes, bzw. Punktewolken 
für jede Stegspitze (siehe Abbildung 2.9). Der Schwerpunkt der Punkte- 
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wolke einer Stegspitze wird für die Validierung als Position dieser Spitze 
angesehen. 


Abbildung 2.9 Ein Beispiel der zur Validierung genutzten Daten mit Punktewolken von 
vier Auflagespitzen. 


Tabelle 2.1 Abweichung zwischen der gemessenen und der wahren Position in y-Richtung 
in Millimetern. 


CoMP11 CoMP45 GaussP 


Durchschnitt -0,047 0,001 -0,063 
Median -0,090 -0,068 -0,097 
Standardabweichung 0,313 0,324 0,317 
Minimum -0,845 -0,713 -0,825 
Maximum 0,688 0,789 0,657 


Die Validierung der Messgenauigkeit wurde mit zehn neuen Stegen 
mit jeweils 51 aufgenommenen Spitzen durchgeführt. Wie oben erwähnt, 
wird in einem ersten Ansatz die Position der Spitze mit der theoretischen 
Position dieser Spitze auf der Palette verglichen. Die Resultate für die 
drei betrachteten Liniendetektionsalgorithmen CoMP11, CoMP45 und 
GaussP sind in Tabelle 2.1 dargestellt. Der CoMP45 Algorithmus hat dabei 
die geringste maximale Abweichung mit + 0,8mm. Die Ergebnisse der 
verschiedenen Algorithmen weichen aber nicht stark voneinander ab, mit 
Differenzen der maximalen Abweichung von 0,056 mm und 0,036 mm. 

Der zweite Ansatz errechnet den Abstand zwischen zwei benachbarten 
Spitzen als euklidische Distanz. Die Verfahren erreichen mit diesem Ansatz 
bessere Ergebnisse (siehe Tabelle 2.2). Dies scheint zu bestätigen, dass die 
Stege eine leichte Torsion aufweisen. Der GaussP Algorithmus hat hier die 
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kleinste maximale Abweichung, aber die anderen beiden Verfahren sind, 
wie beim ersten Validierungsansatz auch, nicht signifikant schlechter. 


Tabelle 2.2 Abstand zwischen zwei Auflagespitzen in Millimetern. 


СоМР11 CoMP45 GaussP 


Durchschnitt 14,798 14,799 14,797 
Median 14,809 14,808 14,797 
Standardabweichung 0,167 0,158 0,157 
Minimum 14,106 14,191 14,258 
Maximum 15,273 15,250 15,171 
Differenz (Maximum - Referenzwert) 0,694 0,609 0,542 
Differenz (Referenzwert - Minimum) 0,473 0,450 0,371 


Abbildung 2.10 Die vier manuell verbogenen Spitzen. Oben: Aufnahme des Steges von 
oben. Unten: Der entsprechende Bereich der Punktewolke. 


Um die Messgenauigkeit bei Spitzen, die nur wenig verbogen wurden, 
sicher zu stellen, wurden vier Spitzen manuell absichtlich verbogen (sie- 
he Abbildung 2.10). Die Auslenkung dieser Spitzen konnte nur manuell 
gemessen werden, was mit einer Genauigkeit von ganzen Millimetern 
möglich war. Um die Auslenkung einer Spitze in der aufgenommenen 
Punktewolke zu berechnen, wurden dieselben oben beschriebenen Verfah- 
ren der Validierungsansätze verwendet. Die manuellen Messergebnisse 
und die Messergebnisse des LTS sind in Tabelle 2.3 gegenübergestellt. Die 
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Tabelle 2.3 Der Abstand zwischen manuell verbogenen und geraden Spitzen. 


SpitzeNr. Manuelle Messung LTS Messung 


1 -3mm -3,33mm 
2 -3mm -3,05 mm 
3 2mm 2,59 mm 
4 3mm 3,59 mm 


Abweichungen entsprechen der zuvor festgestellten Messunsicherheit mit 
Werten um + 0,6mm. Auch wenn eine Stichprobengröße von vier Spit- 
zen natürlich statistisch nicht aussagekräftig ist, zeigt sie, dass es keine 
Probleme gibt, die die Messgenauigkeit an gebogenen Spitzen dramatisch 
verringern. Als Validierungsansatz ist dieses Vorgehen jedoch nicht sinn- 
voll, da das Biegen und manuelle Messen von Hunderten von Spitzen nicht 
durchführbar ist. 


2.3.7 Robustheit des Laser-Triangulationssystems 


Die Robustheit des vorgestellten Messsystems wird anhand der vorher be- 
schriebenen Kategorien - unterschiedliche Stegmaterialien, Verschlackung, 
verbogene Spitzen, aufgeschweißte Butzen, falsch eingesetzte Auflagestege 
und Vibrationen - qualitativ bewertet. Die Probleme mit Artefakten, die auf 
Mehrfachreflexionen zurückzuführen sind, und Methoden zur Milderung 
dieser Probleme wurden bereits in Abschnitt 2.3.4 diskutiert. Aufgrund 
des Messaufbaus im Maschinenkörper verursacht das Umgebungslicht 
keine weiteren Probleme. 

Auflagestege können aus verschiedenen Materialien hergestellt werden, 
meist Kupfer, Edelstahl oder Baustahl. Diese drei Werkstoffe wurden getes- 
tet und alle verhielten sich im Grundsatz gleich. Die stärkere Reflexion an 
der Vorderseite des Stegs, insbesondere bei Edelstahl, führt zu einer allge- 
mein dunkleren Laserlinie, da mehr Licht von der Kamera weg reflektiert 
wird. Dieser Effekt fällt jedoch eher gering aus und deshalb kann leicht 
ein geeigneter Schwellenwert gefunden werden, der für alle Materialarten 
gute Ergebnisse liefert. 
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Einige Auflagestege waren stark verschlackt oder zeigten andere Formen 
von Verschleiß, wie z. В. geschmolzene Spitzen. Ein Beispiel für einen Steg 
mit Verschlackung ist in Abbildung 2.11 dargestellt. Natürlich wird die 
Form der Stege nur von vorne vermessen, da sowohl die Kamera als auch 
die Laserlinie an der Hinterseite abgeschattet sind. Durch die große und 
unregelmäßige Oberfläche war die Erwartung, dass Reflexionen hier ein 
Problem darstellen könnten. Dies konnte nicht bestätigt werden. 


N 


Abbildung 2.11 Ein stark verschlackter Steg zwischen zwei neuen Stegen. 


Dass Auflagespitzen nach vorne oder hinten verbogen sind, kommt beim 
produktiven Einsatz der Maschine oft vor. Dabei ist es für die Betrach- 
tungen dieser Arbeit ausschlaggebend, wie weit eine Spitze verbogen ist. 
Leicht verbogene Spitzen können durch das LTS gemessen werden (vgl. 
Tabelle 2.3 und Abbildung 2.10). Manchmal gibt es auch stark verboge- 
ne Spitzen, wie z.B. in Abbildung 2.12 gezeigt wird. Auch hier hat das 
LIS keine grundsätzlichen Probleme mit der Vermessung. Je nach Aus- 
lenkungen des verbogenen Auflagestegs kann es vorkommen, dass die 
verbogenen Bereiche des Steges durch den Steg davor verdeckt werden. Für 
das Messergebnis ist dies aber nicht entscheidend, da so stark verbogene 
Spitzen keine Auflagepunkte für ein Blech darstellen. 

Eine weitere Quelle möglicher Probleme sind während des Schneid- 
prozesses auf den Auflageleisten verschweißte Schneidbutzen. Bei auto- 
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Abbildung 2.12 Ein stark verbogener Steg zwischen zwei normalen Stegen. 


matisierten Anlagen führen sie dazu, dass das nächste Rohblech nicht 
korrekt auf die Palette gelegt werden kann und durch den Butzen eine 
Wölbung nach oben bekommt. Dies kann zu vermehrten Vibrationen und 
Maßungenauigkeiten führen. In den Aufnahmen ist ein solcher Butzen ent- 
halten (siehe Abbildung 2.13). Die Detektion von aufgeschweißten Butzen 
ist in der Punktewolke einfach durchzuführen. Dafür muss nur die An- 
zahl der Punkte im relevanten Bereich der z-Achse gezählt und mit einem 
Schwellenwert, der nicht überschritten werden darf, verglichen werden. 


u 


Abbildung 2.13 Ein sechseckiger Butzen, der an einer Stegspitze verschweißt ist. 


Da es ein wenig Kraft erfordert, einen Steg in die Stegsteckposition 
einzulegen und dies bei allen LFMn manuell durchgeführt wird, kann es 
vorkommen, dass Stege nicht korrekt eingelegt sind und stattdessen zu 
hoch stehen. Dies führt wie bei Butzen zu einer möglichen Wölbung des 
Blechs. Genau wie bei auf den Auflagestegen verschweifsten Butzen kann 
dies leicht durch einen Schwellenwert im relevanten Bereich der z-Achse 
erkannt werden. In Abbildung 2.14 werden zwei solche Stege gezeigt. 

Die Bilder wurden auch auf weitere mögliche Anzeichen von abneh- 
mender Robustheit wie z. B. Vibrationen untersucht. Dazu wurden meh- 
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Abbildung 2.14 Zwei nicht korrekt eingelegte Stege, deren Spitzen höher stehen als die 
umgebenden Spitzen. 


rere Tests durchgeführt, wie 7. В. ob es eine periodische Schwankung der 
Ausrichtung der Bilder einer Sequenz gibt. Es konnten keine negativen 
Auswirkungen auf die Messung gefunden werden. 


2.4 Stegdetektion auf einzelnen Bildern im 
sichtbaren Spektrum 


Im Vergleich zur Laser-Iriangulation bietet die Detektion der Stege auf 
einem einzelnen Bild der gesamten Palette der Flachbettmaschine den 
Vorteil, dass die Aufnahme weniger Zeit benötigt und sich reibungslos in 
den Produktionsprozess integrieren lässt, da die Palette nach dem Entfer- 
nen der geschnittenen Teile und des Restgitters auf dem Palettenwechsler 
leer ist. Dieser Moment kann für eine Aufnahme genutzt werden (siehe 
Abbildung 2.15). 
In [86] wurden bereits Teile der Inhalte dieses Unterkapitels publiziert. 


2.4.1 Detaillierte Problembeschreibung 


Bei den Bildern der leeren Paletten sind einige Störeinflüsse zu berücksich- 
tigen, die diese Detektionsaufgabe deutlich erschweren. Da die Kamerapo- 
sition auf die Rückseite der Maschine beschränkt ist, können die Bilder nie 
aus der Vogelperspektive aufgenommen werden. Entsprechend sind die 
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Abbildung 2.15 Oben: Ein Beispiel eines Bildes der Palette. Die Längsträger sind besonders 
im linken Drittel zu erkennen. Unten links: Detailausschnitt aus der linken oberen Ecke. 
Unten rechts: Detailausschnitt aus der rechten oberen Ecke. 


vorderen Stege der Palette fast direkt unter der Kamera und werden von 
oben gesehen (siehe Abbildung 2.15 unten links), während die hinteren 
Stege deutlich von der Seite gesehen werden (siehe Abbildung 2.15 unten 
rechts). Dadurch werden auch die verschiedenen Verschlackungszustände 
bei den hinteren Stegen verstärkt aufgenommen. 

Außerdem können auf dem Palettenwechsler die beiden Paletten der 
Maschine gleichzeitig übereinander stehen. Dann sind durch die obere 
Palette auch die Stege der unteren Palette sowie Bleche, die möglicherweise 
auf der unteren Palette liegen, sichtbar. Ein ähnlicher Störeinfluss ergibt sich 
durch Bleche oder Dreck auf dem Fußboden unter dem Palettenwechsler, 
der insbesondere durch die vorderen Stege zu sehen ist (siehe Abbildung 
2.16). 

Ein weiterer störender Faktor sind die unkontrollierten Lichtbedingun- 
gen. Da die LFM nicht samt Palettenwechsler eingehaust werden kann, ist 


41 


2 Vermessung der Auflagestege von Laserflachbettmaschinen 


ill | 
|| 


et 


Abbildung 2.16 Beispiele für die Störeinflüsse auf den Grauwertbildern. Links: Verschiedene 
Verschlackungszustände. Mitte: Zwei Paletten übereinander, dabei liegt links ein Blech auf 
dem Boden und rechts ein Blech auf der unteren Palette. Rechts: Die Spitzen und Senken 
benachbarter Stege sind im rechten Teil des Bildes sehr nah beieinander. 


die Qualität der Aufnahmen bis zu einem gewissen Grad von der Beleuch- 
tungssituation in der Produktionshalle abhängig. Hier reicht das Spektrum 
von spärlichster Beleuchtung während der Nachtschicht bis zu direktem 
Sonnenlicht auf der Palette zur Mittagszeit. 


2.4.2 Stand der Technik 


Das Template Matching ist eine der verbreitetsten Methoden zur Detektion 
von bekannten Objekten in Bildern. Hierbei wird zuerst eine Schablone 
erstellt, die entweder aus dem Objekt selbst besteht, wenn dieses nur ei- 
ne einzige Erscheinungsform hat. Bei Stegen ist dies z. B. auf Grund von 
Verschlackung und unterschiedlicher Beleuchtung nicht der Fall. Hier 
kann auch eine gemittelte Erscheinungsform des Objekts herangezogen 
werden. Die Schablone wird mit dem Bild gefaltet und Ergebnisse über 
einem gewissen Schwellenwert, bzw. das maximale Ergebnis, werden als 
Ort des gesuchten Objekts gezählt [7]. Statt einer Faltung, die als Kreuzkor- 
relation formuliert werden, kann auch die normalisierte Kreuzkorrelation 
berechnet werden. Letztere ist weniger anfällig für allgemeine Helligkeits- 
oder Kontrastunterschiede zwischen dem zu untersuchenden Bild und der 
Schablone und erleichtert die Suche nach einem passenden Schwellwert 
[14]. 

Eine weitere Analysemethode, die sich auf die Textur eines Bildbereichs 
fokussiert, sind Laws’ Energiemaße [54]. Solange das gesuchte Objekt eine 
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deutlich unterscheidbare Textur im Vergleich zu allen anderen Bildberei- 
chen hat, kann dadurch das Objekt erkannt werden. 

Laws’ Energiemaße bestehen dabei aus einer Menge von quadratischen 
Matrizen variabler Größe. Am verbreitetsten ist der Einsatz von Matrizen 
der Größe 5 x 5, da diese in numerischen Studien die besten Ergebnisse 
geliefert haben [54]. Die Matrizen werden aus dem dyadischen Produkt 
von vier Vektoren gebildet, die Geraden (15), Hochpunkte (35), Wellen (rs) 
und Kanten (e,) repräsentieren: 


15 Ke (1,4,6,4, 15, S5 = (—1, 0,2,0, —1)Т, 
г5 = (1, —4, 6, —4, LP: еѕ = (—1, —2,0, 2, EE 


Die 16 Matrizen werden jeweils mit dem Bild gefaltet und erstellen 
dadurch 16 Features. Ein Objekt mit bekannter Textur kann nun in diesen 
Features gesucht werden. Eine weitere Möglichkeit besteht darin, den 
Durchschnitt oder einen gewichteten Durchschnitt der Matrizen zu bilden 
und nur dort nach der bekannten Textur zu suchen. 

Viele andere bekannte Ansätze aus den letzten Jahren sind für das vorlie- 
gende Problem nicht anwendbar. Zum Beispiel hat der SIFT Algorithmus 
die Objektdetektion bei skalierten und gedrehten Objekten durch eine 
geschickte Darstellung der Features entscheidend vorangebracht [59]. Die 
Stege sind in den Bildern aber nie gedreht oder skaliert vorhanden. 

Einer der neueren und in der Bildverarbeitung außerordentlich erfolg- 
reichen Ansätze ist eine Klasse von speziellen künstlichen neuronalen 
Netzen, die Convolutional Neural Networks (СММ). Die ersten Ideen dieses 
Ansatzes gehen einige Jahrzehnte zurück, aber die ersten Beschreibungen 
des heute angewandten Konzepts finden sich in den späten 1990er Jahren 
(siehe z. B. [55]). Aktuelle Beiträge zur Objektdetektion und -lokalisierung 
mit CNNs wurden in [24] und [35] geleistet. Dabei ist zu beachten, dass 
für das Training eines CNNs 1. 4. К. mehrere tausend Trainingsbeispiele 
notwendig sind, da sie oft eine große Anzahl an Parametern haben. Da nur 
ungefähr 220 Bilder vorliegen, wurde dieser Ansatz nicht weiter verfolgt. 
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2.4.3 Methoden zur Stegdetektion auf einzelnen 
Grauwert-Aufnahmen 


2.4.3.1 Kanten- und Eckendetektoren und Laws’ Energiemaße 


Die Stege können als eine bestimmte Anordnung von Kanten und Ecken im 
Bild aufgefasst werden. Deshalb erscheint es sinnvoll, Kanten- und Ecken- 
detektoren zu betrachten. Als Kantendetektoren werden die in Kapitel 2.1 
vorgestellten Verfahren mittels differenziertem Gauß-Tiefpass oder DoG 
als Eckendetektoren, die respektive die erste oder zweite Ableitung der 
Grauwertfunktion nutzen, verwendet. Als oft gebrauchter Eckendetektor 
wird der Harris Corner Detector benutzt, der ebenfalls in Kapitel 2.1 vor- 
gestellt wurde, und als Texturanalyseverfahren die im vorhergehenden 
Abschnitt eingeführten Laws’ Energiemaße. 


Abbildung 2.17 Beispiele für drei Features. Links: Differenzierter Gauß-Tiefpass. Mitte: 
Laws’ Energiemaße. Rechts: Binarisierter Harris Corner Detector. 


Laws’ Energiemaße werden mit Matrizengröße 5 x 5 gebildet. Die Ma- 
trix, die aus der Multiplikation von 1, mit sich selbst entsteht, wird dabei 
vernachlässigt, da sie nur einen gewichteten Durchschnitt bildet. Nach 
der Faltung der 15 Matrizen mit dem Bild wird hochfrequentes Rauschen 
durch Faltung mit einem Gauß-Tiefpassfilter der Größe 5 x 5 gedämpft. Da 
in Vorversuchen gezeigt wurde, dass die einzelnen Matrizen nur wenig 
Information enthalten, wird anschließend wird der Durchschnitt über die 
15 Matrizen gebildet. 


2.4 Stegdetektion auf einzelnen Bildern im sichtbaren Spektrum 


2.4.3.2 Fregenzbasierte Features 


Da die Stegspitzen einen regelmäßigen und vorab bekannten Abstand 
haben, kann das Bild mittels diskreter Fouriertransformation in seine Fre- 
quenzen zerlegt und anschließend gesuchte Frequenzen gefiltert werden. 
Nach einer Rücktransformation sind dann nur noch die Bildanteile enthal- 
ten, die der gefilterten Frequenz entsprechen. 

Da die Stege gepunktete, vertikale Linien im Bild sind und die Spitzen 
einen bestimmten Abstand zueinander haben, sind im Fourier-Betragsqua- 
dratspektrum horizontale Linien mit schwankender Intensität zu sehen. 
Die Distanz d der ersten parallelen Linien zur Achse kann aus dem verti- 
kalen Abstand zweier Spitzen d. der Höhe des Bildes umax und der Größe 
eines Pixels Au berechnet werden: 


d = = Au: Umax: 
5 

Die prognostizierten horizontalen Linien sind in Abbildung 2.18 ein- 
deutig zu erkennen und treten im erwarteten Abstand auf. Die auffälligen 
Strukturen im Betragsquadratspektrum werden also durch die Stegspitzen 
ausgelöst. Da es noch weitere Oberschwingungen gibt, die durch weitere 
parallele Linien erkannt werden können, werden drei Bandpass-Filter (BP1 
bis BP3) definiert, um das Signal zu filtern: 


1, wenn 30 < |f,| < es und |f,| <5 
BPI Sen fy) =] 0 sonst IF | А 


1, wenn |, < е, und а-5< < 4+5 
BP2( fp, fy) =] А Е If, | 1 ЈА 1 


1, wenn |/,| < e, und 
0, sonst. 


Die Ausdehnung der Filter ist dabei durch eine Grenze е, gegeben, da 
die Energie des Signals für höhere Frequenzen abnimmt. Aus dem selben 
Grund werden auch nur drei Filter betrachtet. Die Bandpass-Filter sind in 
Abbildung 2.18 als weiße Overlays eingezeichnet. 

Nach der Filterung werden die Bilder in den Ortsbereich zurück transfor- 
miert, wodurch ein Bild entsteht, das nur noch die gefilterten Frequenzen 
enthält. Ein Beispiel ist in Abbildung 2.19 zu sehen. 
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Abbildung 2.18 Links: Das Betragsquadratspektrum eines Bildes mit Stegen. Rechts: Das- 
selbe Betragsquadratspektrum mit den drei Bandpassfiltern als weiße Overlays. 


Abbildung 2.19 Resultat der Filterung im Frequenzbereich. Von oben nach unten: Grau- 
wertbild, Filterung mit BP1, Filterung mit BP2, Filterung mit BP3. 


2.4.3.3 Template Matching Klassifikator 


Template Matching basiert auf der Faltung einer Schablone mit dem Bild. 
Die Schablone muss im vorliegenden Anwendungsfall also geeignet sein, 
um alle Stege zu erkennen, und darf trotzdem gleichzeitig nicht zu anfällig 
für Störeinflüsse sein. 

Die Schablonen für das Template Matching wurden alle von Hand erstellt 
und auf den jeweiligen Merkmalsraum angepasst. Bei den einfachen Grau- 
wertbildern kommt das Bild eines neuen Steges aus der Mitte der Palette 
zum Einsatz. Bei den Fourier-Features kommen sowohl ein Bandpass- 
gefilterter Steg (BP2 Steg) zum Einsatz als auch zwei binäre Schablonen 
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(BM, und ВМ.), die anhand der geeigneten Breite der Stege w, folgender- 
maßen definiert sind: 


BM, (u,v) = 1, wenn sw, <u < 3, 
5 0, wenn 0 < v < zw, und zw, S v < ws 
1, wenn v < tw 
BM,(u, 0) = 0, wenn tw ua 
р 9 5 > Eh 


Sie sollen die Erscheinung von Stegen im gefilterten Bild nachahmen, näm- 
lich zum einen einen hellen vertikalen Streifen zwischen zwei dunklen 
Bereichen und zum anderen einen hellen vertikalen Streifen, bei dem nur 
rechts ein dunkler Bereich anschließt. Letzteres ist der Tatsache geschuldet, 
dass bei Stegen im rechten Bereich des Bildes nicht nur die Spitzen, son- 
dern auch die Senken sichtbar sind und aufgrund der Distanz zwischen 
Senken und Spitzen auch ein sehr schmaler dunkler Streifen zwischen den 
Spitzen und Senken eines Stegs vorhanden sein kann. Ein symmetrisches 
Template wie BM, kann dann dazu führen, dass falschlicherweise zwei 
Stege erkannt werden. Entsprechend ist der Ansatz bei ВМ. nur dann po- 
sitiv zu klassifizieren, wenn neben einem hellen Streifen auch ein dunkler 
Streifen mit mindestens halber Stegbreite folgt. Entsprechend lässt sich 
auch das Template eines bandpass-gefilterten Steges anpassen. Hier wird 
dem bestehenden Template einfach ein schwarzer Streifen rechts ange- 
fügt, was im Folgenden auf Grund der Asymmetrie als „Asym. BP2 Steg” 
abgekürzt wird. 

Um die Ähnlichkeit zwischen der Schablone und dem Bild zu bewerten, 
wird die normalisierte Kreuzkorrelation benutzt. Das Resultat ist ein ein- 
dimensionales Signal. Da Stege nur alle 33,5 mm vorhanden sein können, 
wird gefordert, dass ein Steg ein lokales Maximum im Bereich + 15mm 
um die detektierte Position ist. Um bei unbesetzten Stegpositionen keine 
Fehldetektionen zu erhalten, muss außerdem ein gewisser Schwellenwert 
überschritten werden. 


2.4.3.4 Parallelprojektions-Klassifikator 


Ein weiterer möglicher Klassifikator entsteht durch die Parallelprojekti- 
on. Da die Stege in erster Näherung als vertikale Linien in den Bildern 
dargestellt sind, ist es möglich, eine Parallelprojektion entlang der Stege 
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durchzuführen und aus dem entstandenen Signal die Spitzen zu klassifi- 
zieren. 

Letzteres ist allerdings kein einfaches Problem, da das Signal, das aus der 
Parallelprojektion resultiert, z.B. durch Beleuchtungsunterschiede oder 
durch hellere Schlacke auf den Stegen beeinflusst wird. Hier wird die in 
[20] vorgestellte Wavelet-Iransformation basierte Pattern-Matching Me- 
thode verwendet, um die Maxima zu detektieren. Dazu wird zuerst die 
Faltung von diskreten Ricker Wavelets mit 20 verschiedenen Breiten mit 
dem Signal berechnet. Die Ergebnisse dieser Faltung werden in einer Matrix 
festgehalten, wobei jede weitere Zeile das nächst größere Wavelet enthält 
und jede Spalte zu der entsprechenden Stelle des diskreten Signals gehört. 
Dann werden in der Matrix die Maxima in jeder Reihe mit denen in den 
benachbarten Reihen verglichen und Maxima über mehrere Reihen hinweg 
zu sogenannten Ridge Lines zusammengefasst. Diese Ridge Lines werden 
dann einem Pattern-Matching unterzogen, um Stege aus den möglichen 
Maxima zu detektieren. Dazu wird die Form und Länge der jeweiligen 
Ridge Line herangezogen. 

Der Vorteil dieser Methode liegt darin, dass Maxima verschiedener Höhe 
leicht erkannt werden können und die Form der Maxima berücksichtigt 
wird, was zu einer höheren Informationseffizienz führt [20]. 


2.4.3.5 Transformation in Stegpositionen 


Ein Steg kann nur an ganz bestimmten Positionen auf der Palette sein, da 
es regelmäßige Stegsteckplätze gibt. Daher ist eine Transformation von 
der Bildspalte eines erkannten Steges auf einen Stegsteckplatz notwendig. 
Bei der später präsentierten Umsetzung handelt es sich um 3100 Bildspal- 
ten und 93 Stegsteckplätze. Da die Palette in jedem Bild unterschiedlich 
ausgerichtet ist, muss zunächst die Position des Blechanschlags Tstop er- 
kannt werden. Die Distanz zwischen dem Blechanschlag und dem ersten 
Stegsteckplatz ist ungefähr gleich der Distanz zwischen zwei Stegsteck- 
реп d, Die Position des Blechanschlags ist als Referenzpunkt also 
ausreichend und wird mittels Template Matching in einem definierten 
Bereich in der entsprechenden Ecke des Bildes festgestellt. Die Schablone 
dafür wird für jede Maschine neu von Hand erstellt. 
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Der Stegsteckplatz п € {1,2,...,92,93} einer bestimmten Bildspalte v,, 


kann durch: 
А (= = “stop ) 
n = nint | ————— 


st 


berechnet werden. 


2.4.4 Umsetzung, Ergebnisse und Diskussion 


Die Aufnahme der Bilder wurde an LFMn der Baureihen TRUMPF TruLa- 
ser 3000 und 5000 gemacht, deren Palettenwechsler fast identisch aufgebaut 
sind. Insgesamt wurden 215 Bilder verwendet, wovon 27 Bilder von einer 
Testmaschine bei TRUMPF Werkzeugmaschinen stammen und 188 Bilder 
von Blechfertigern. 

Der Kameraarm war dabei an der Maschinenrückseite angebracht, so- 
dass die gesamte Palette mit einer Kamera betrachtet werden kann (siehe 
Abbildung 2.20). 


Abbildung 2.20 Die Position des Kameraarms an der Maschine. 


Die Ergebnisse der Methoden bei Benutzung der Parallelprojektion bzw. 
des Template Matchings als Grundlage der Klassifizierung sind in den 
Tabellen 2.5 und 2.4 zu sehen. Dabei ist die Accuracy wie in Formel 2.3 
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Tabelle 2.4 Ergebnisse der verschiedenen Features bei Klassifizierung auf Grundlage des 
Template Matching. 


Feature Template КРК КМК Accuracy 
Grauwerte Neuer Steg 0,204 0,642 0,434 
BP2 BP2 Steg 0,948 0,948 0,948 
BP2 Asym. BP2 Steg 0,966 0,958 0,961 
BP1 BM, 0,764 0,715 0,738 
BP2 BM, 0,921 0,857 0,887 
BP3 BM, 0,930 0,870 0,898 
BP1 BM, 0,416 0,649 0,538 
BP2 BM, 0,926 0,947 0,937 
BP3 BM, 0,922 0,930 0,926 


definiert. Die Klassifikationsergebnisse sind mit ihren Anfangsbuchstaben 
abgekürzt, sodass z.B. RP für richtig positiv und FN für falsch negativ 
steht. Die Richtig Positiv Rate (RPR) ist dabei der Anteil der positiven 
Instanzen, die korrekt als solche erkannt wurden. Die Richtig Negativ Rate 
(КМК) ist analog definiert. 


RP +RN 
Accuracy = RPAFP+FNIRN (2.3) 

Grundsätzlich schneiden die Fourier-Features deutlich besser ab als 
jede andere Methode, mit einer Accuracy von i. а. R. rund um 90%. Die 
nächst beste Methode ist DoG mit 72,5 % Accuracy. Grauwerte ohne weitere 
Features zu benutzen, funktioniert mit Parallelprojektion besser als mit 
Template Matching. Letzteres ist mit einer Accuracy von 66,7 % eine der 
schlechteren Methoden. Dies zeigt, dass kein gutes Template gefunden 
werden kann, was dafür spricht, dass die Störeinflüsse tatsächlich eine 
große Rolle spielen. 

Ähnlich verhält es sich mit den Ecken- und Kantendetektoren, die zu oft 
durch Rauschen, entweder durch Strukturen unter der Palette oder durch 
verschlackte Stege, zu fehlerhaften Detektionen gelangen. Sie erreichen nur 
eine Accuracy zwischen 49,7 % und 72,5 %. Hier zeigt sich, wie sensibel 
Verfahren sind, die den Spitzenabstand nicht berücksichtigen können. 
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Tabelle 2.5 Ergebnisse der verschiedenen Features bei Klassifizierung auf Grundlage der 
Parallelprojektion. 


Feature КРК КМК Accuracy 
Grauwerte 0,860 0,486 0,667 
BP1 0,855 0,412 0,625 
BP2 0,916 0,655 0,781 
BP3 0,897 0,805 0,850 
Laws Energy Measures 0,903 0,453 0,670 
Harris Corner Detector 0,605 0,396 0,497 
Рос 0,880 0,581 0,725 
Differenzierter Gauß-Tiefpass 0,807 0,336 0,562 


Die meisten Fehler werden bei jedem Verfahren im rechten Drittel des 
Bildes gemacht. Die Fehlerkurve steigt insbesondere im hinteren Drittel 
deutlich an. Dies ist unabhängig von den eingesetzten Methoden, sodass 
auch mit Informationsfusionsverfahren kein Fortschritt erzielt werden 
kann. Eine erste Vermutung könnte sein, dass es einen Zusammenhang 
zwischen der durch die perspektivische Transformation des aufgenomme- 
nen Bildes in diesem Bereich geringere Auflösung und den Fehlklassifi- 
kationen gibt. Bei der Angleichung der Auflösung über das gesamte Bild 
wird allerdings nicht im letzten Bereich interpoliert, sondern die Auflö- 
sung des vorderen Bereichs herabgesetzt. Die Auflösung liegt mit 1 mm 
pro Pixel auch deutlich über dem Spitzenabstand von 14,8mm. Obwohl 
zu erwarten wäre, dass die Fehldetektionen abnehmen, wenn die Scha- 
blonen gezielt aus diesem Bereich der Palette stammen, ist dies nicht der 
Fall. Der Grund dafür ist wahrscheinlich, dass insbesondere in diesem 
Teil der Palette die Seiten der Stege gut zu sehen sind und damit auch 
der Verschlackungszustand einen großen Einfluss hat. Ebenso führt die 
perspektivische Überlappung von aufeinander folgenden Stegen in diesem 
Bereich des Bildes zur Verwechselung von Stegsenken und Stegspitzen. 
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2.5 Stegdetektion durch schmalbandige 
Nahinfrarot-Aufnahmen 


Die im vorhergehenden Kapitel angesprochenen Probleme haben zu nicht 
unerheblichen Teilen mit der unkontrollierten Beleuchtungssituation zu 
tun. Die Beleuchtung ist bei vielen Anwendungen der automatischen Sicht- 
prüfung entscheidend für Robustheit und Geschwindigkeit der Auswer- 
tung [30]. Der Ansatz einer NIR-Dunkelfeldbeleuchtung soll die stören- 
den Fremdeinflüsse bei Einzelbildern reduzieren. Dadurch sollten auch 
Aussagen über den Zustand einzelner Spitzen möglich werden, während 
gleichzeitig die Auswertungsgeschwindigkeit steigt. Da nicht die gesamte 
Form der Stege, sondern nur die Position der Spitzen von Interesse für die 
Messaufgabe ist, ist ein Bild der Palette denkbar, das im Hintergrund voll- 
ständig schwarz ist und nur die Spitzen als weiße Punkte in geraden Linien 
zeigt. Dies ist das Ziel der vorgeschlagenen NIR-Dunkelfeldbeleuchtung. 
Die Auswertung eines solchen Bildes ist offensichtlich mit einfachen Algo- 
rithmen möglich, die im Folgenden vorgestellt werden. 

Die in diesem Unterkapitel dargestellten Inhalte sind in ähnlicher Form 
bereits in [88] veröffentlicht worden. Um eine bessere Lesbarkeit zu errei- 
chen, wird auf wörtliche Zitate dieser Publikation verzichtet. 


2.5.1 Stand der Technik 


Die Dunkelfeldbeleuchtung ist in der automatischen Sichtprüfung ein oft 
benutztes Inspektionsverfahren [1, 7]. In einem üblichen Aufbau wird die 
Beleuchtung in einem sehr flachen Winkel bezüglich der zu überprüfenden 
Oberfläche platziert, über der senkrecht die Kamera steht. Das meiste Licht 
wird nicht zur Kamera geleitet, sondern in einem flachen Winkel von 
der Beleuchtung weg, wodurch ein dunkles Bild entsteht. Defekte in der 
Oberfläche streuen das Licht, das so in die Kamera geleitet wird. Defekte 
erscheinen also heller als ihre Umgebung. 

In [1] wird diese Beleuchtungsart gewählt, um Defekte in Stoffgewe- 
ben vor der Produktion feststellen zu können. In einer anderen Studie 
wurde die Dunkelfeldbeleuchtung gewählt, um ein dünnes Gewebe zu 
erfassen und daraus die besten Schneidwege des Werkzeugs zu bestimmen 
[4]. Ein ebenfalls flächiges Ausgangsmaterial, nämlich ein Separator für 
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Lithium-Ionen-Batterien, wurde in [43] betrachtet. Hier wurde neben der 
Dunkelfeldaufnahme auch eine Hellfeldaufnahme durchgeführt und die 
gewonnene Information hinterher zur Defektdetektion zusammengeführt. 

Andere Studien benutzen die Dunkelfeldbeleuchtung auch im nicht 
sichtbaren Spektrum. In [23] werden Bilder von Composit-Materialien 
mit Röntgen-Strahlung aufgenommen. Die Proben wurden vorher starker 
mechanischer Belastung unterzogen, um die Kollision eines Bauteils mit 
einem Gegenstand zu simulieren. Die entstandenen Bilder können genutzt 
werden, um Defekte im Composit-Material, die durch die mechanische 
Belastung ausgelöst wurden, zu detektieren. 

Die Unterdrückung von Fremdlichteinflüssen bei automatischen Sicht- 
prüfungsaufgaben ist ein breit erforschtes Gebiet. In [57] wurde zur Un- 
terdrückung von Fremdlichteinflüssen im sichtbaren Spektrum bei der 
Gesichtserkennung eine schmalbandige Aufnahme und eine Beleuchtung 
im nahen Infrarot benutzt. Der Vorteil einer NIR-Beleuchtung ist, dass die 
Unterdrückung von Fremdlichteinfliissen durch Sonnenlicht oder künstli- 
che Beleuchtung auch ohne Einhausung sichergestellt werden kann. 

Für die Assistenzsysteme von Automobilen sind NIR-Beleuchtungen 
ebenfalls betrachtet worden. So kann durch eine aktive NIR-Beleuchtung 
mit bestimmten Wellenlängen eine recht genaue Sichtweitenschätzung 
abgeben werden, ohne dabei durch Sonnenlicht gestört zu sein oder nachts 
durch sichtbare Beleuchtung die Verkehrsteilnehmer zu stören [12]. 

Arbeiten über den Einsatz von NIR-Beleuchtung zur Maschinenüber- 
wachung oder anderen Aufgaben in der Blechfertigung sind dem Autor 
nicht bekannt. 


2.5.2 Konstruktion einer schmalbandigen 
NIR-Dunkelfeldbeleuchtung 


Der Aufbau einer schmalbandigen Dunkelfeldbeleuchtung an einer LFM 
soll vor allem zwei Probleme lösen. Erstens soll mittels der schmalbandigen 
NIR-Komponente Fremdlicht sehr gut unterdrückt werden, sodass eine 
Dunkelfeldbeleuchtung ohne Einhausung ermöglicht wird. Zweitens soll 
durch die Dunkelfeldbeleuchtung die Analyse der Bilder sehr viel einfacher 
und dadurch auch schneller möglich werden. Die angestrebte Aufnahme ist 
zum Großteil schwarz und zeigt die Spitzen eines Steges als helle Punkte. 
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Die Umsetzung einer schmalbandigen NIR-Beleuchtung und -Aufnahme 
ist durch eine entsprechende Beleuchtung und Filter vor dem Kamera- 
objektiv möglich. Um eine möglichst gute Fremdlichunterdrückung zu 
erzeugen, sind die Wellenlängen 850 nm und 940 nm besonders geeignet, 
da sie auf der Erdoberfläche Bandlücken im Spektrum des Sonnenlichts 
sind. Dies hängt mit der Absorption durch Wassermoleküle in der At- 
mosphäre zusammen [32]. Deshalb sind diese Wellenlängen auch in der 
Sicherheitstechnik verbreitet und entsprechende Filter verfügbar. 

Die Konstruktion einer Dunkelfeldbeleuchtung für die Palette einer LFM 
unterliegt mehreren Anforderungen. Die Kamera muss an der Rückwand 
der Maschine angebracht werden, da die drei anderen Seiten sowie der 
Raum über der Palette für Be- und Entladeprozesse der Maschine zur 
Verfügung stehen müssen. Die Lichtquelle muss zum einen oberhalb der 
Stegspitzen positioniert werden, damit reflektiertes Licht nach oben zur 
Kamera geleitet wird. Zum anderen darf die Lichtquelle nicht in einem zu 
steilen Winkel von oben auf die Palette gerichtet sein, dasonst Gegenstände, 
die sich unter der Palette befinden, beleuchtet werden könnten. 

Außerdem darf die Lichtquelle nicht die der Kamera zugewandte Seite 
der Stege beleuchten. Sonst entstehen Bilder, die Bildern aus dem sichtbaren 
Spektrum sehr ähnlich sind und entsprechend auch die in Kapitel 2.4 
angesprochenen Probleme bei der Auswertung haben. 

Die Beleuchtung muss also seitlich, in einem flachen Winkel von oben 
und hinten erfolgen (vgl. Abbildung 2.24). Tests haben gezeigt, dass es 
zu teilweise starken Blendenflecken kommt, wenn Licht direkt von der 
Beleuchtung in die Kameraoptik gelangen kann. Dies ist vor allem dann der 
Fall, wenn sich die Beleuchtung auf der der Kamera gegenüberliegenden 
Seite der Palette befindet. Die Beleuchtung muss dann mit entsprechenden 
Abdeckungen versehen werden. 


2.5.3 Auswertungsmethoden für schmalbandige 
NIR-Aufnahmen 

Durch die Position der Kamera sind die Aufnahmen perspektivisch ver- 

zerrt. Stege und Stegspitzen, die weiter vorne auf der Palette sind und 


damit im Bild weiter unten erscheinen, sind größer dargestellt als Ste- 
ge und Stegspitzen weiter hinten. Die nötige Entzerrung kann mit einer 
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Homografie durchgeführt werden. Dazu müssen manuell vier Punkte im 
Bild bestimmt werden, die die Ecken eines Rechtecks sind. Diese Punkte 
sind durch die Stegspitzen sehr einfach manuell zu finden. Da alle Bilder 
aus derselben Perspektive aufgenommen werden, muss dies nur einmal 
geschehen und kann danach auf alle Bilder übertragen werden. Beim ge- 
richteten Bild sind sowohl die Abstände zwischen zwei Stegen als auch 
die Abstände zwischen zwei Spitzen bei verschiedenen Stegen nahezu 
konstant. 


Abbildung 2.21 Ein Beispiel einer gerichteten NIR-Dunkelfeldaufnahme von Auflagestegen 
einer LFM. 


Ob ein Steg eingelegt worden ist, lässt sich leicht mittels Parallelprojek- 
tion bestimmen. Die Grauwerte werden dabei reihenweise (von links nach 
rechts in Abbildung 2.21) addiert und es entsteht ein diskreter Verlauf der 
Helligkeitswerte im Bild. Stege stellen in diesem Signal ein lokales Maxi- 
mum dar. Die Maxima lassen sich mit einfachsten Methoden erkennen. So 
reicht es, Maxima durch den Vergleich mit den zwei direkten Nachbarn 
zu identifizieren. Aus dieser recht umfänglichen Liste werden dann alle 
Maxima entfernt, deren Differenz zwischen ihrem Wert und dem Wert 
des nächsten lokalen Minimums einen definierten Schwellenwert nicht 
erreicht. Zusätzlich werden so lange kleinere Maxima entfernt, bis alle Ma- 
xima eine gewisse Mindestdistanz zueinander haben. Die verbleibenden 
Maxima werden als Stege gewertet. 


55 


2 Vermessung der Auflagestege von Laserflachbettmaschinen 


Aufgrund von Punktreflexionen (siehe Abbildung 2.22) kann es mit die- 
sem Verfahren zu Fehldetektionen kommen. Um ihre Anzahl zu senken, 
kann außerdem für jede als Steg klassifizierte Reihe geprüft werden, wie 
stark die Grauwerte entlang der Reihe fluktuieren. Dazu wird im Abstand 
von wenigen Pixeln die Differenz der Grauwerte gebildet. Für die Imple- 
mentierung wurde ein Abstand von 10 Pixeln gewählt, wobei leicht andere 
Werte zu den selben Ergebnissen führen. Wenn der Durchschnitt dieser 
Differenzen über einem definierten Schwellenwert liegt, handelt es sich um 
einen Steg. Wenn der Durchschnitt unter dem Schwellenwert liegt, handelt 
es sich um eine Fehlklassifikation, die aus der Liste der detektierten Stege 
gelöscht wird. In den in Kapitel 2.5.2 beschriebenen Aufnahmen dieser 
Arbeit ergaben sich zwischen den korrekten und fehlerhaften Detektionen 
zwei klare Cluster, weshalb ein Schwellenwert einfach bestimmt werden 
konnte. Als Beispiel ist das aus der Parallelprojektion resultierende Signal 
für das in Abbildung 2.21 gezeigte Beispiel in Abbildung 2.23 dargestellt, 
wobei die nach der Filterung von Punktreflexionen erkannten Stege rot 
markiert sind. 


Abbildung 2.22 Ein Beispiel von zwei der manchmal auftretenden Punktreflexionen, zwi- 
schen drei Stegen. 


Aus dem extrahierten Steg können auch die einzelnen Spitzen detektiert 
werden, indem wieder eine Parallelprojektion durchgeführt wird. Dazu 
werden nur jeweils die 11 Zeilen des Bildes berücksichtigt, bei denen 
zuvor die mittlere Zeile als Steg erkannt wurde. Die Parallelprojektion 
wird diesmal in Richtung der Spalten durchgeführt, wodurch wie zuvor 
bei den Stegen auch ein diskretes Signal entsteht. Jede Stegspitze sollte in 
diesem Signal ein lokales Maximum darstellen. Wie bei der Stegerkennung 
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Abbildung 2.23 Das aus der Parallelprojektion resultierende Signal für das Bild aus Abbil- 
dung 2.21. Die nach der Filterung von Punktreflexionen erkannten Stege sind rot markiert. 


zuvor wird auch hier bei der Detektion von Maxima erstens ein minimaler 
Abstand zwischen zwei Maxima und zweitens eine minimale Differenz 
zum nächsten lokalen Minimum gefordert. 


2.5.4 Umsetzung 


Der Aufbau wurde wie in Kapitel 2.5.2 beschrieben an einer LFM der Serie 
TRUMPF TruLaser 3000 realisiert (siehe Abbildung 2.24). Die verwendete 
Kamera war eine FLIR Blackfly S mit Computar M0824-MPW2 Objektiv mit 
8 mm Brennweite. Die eingesetzten Schmalbandfilter haben eine maximale 
Durchlässigkeit von über 90 % bei 850 nm + 5 пт bzw. 940 nm + 5nm. Als 
Beleuchtung wurde eine breitbandige Rotlicht-Infrarot-Lampe für dem 
medizinischen Gebrauch genutzt. Die Lampen waren dabei so weit hinter 
und über den Stegen positioniert, dass der Einfallswinkel von hinten zwi- 
schen 15° und 20° und bei circa 15° von oben lag. Die Stegsteckplätze sind 
bei dieser Maschine 33,5 mm von einander entfernt. Der Abstand zwischen 
zwei Spitzen beträgt 14,8 mm. 

Die Kamera wurde mit dem verbreiteten Schachbrettmusterverfahren 
intrinsisch kalibriert [82]. Vor Beginn der Aufnahmen wurde durch Aufnah- 
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Abbildung 2.24 Der Versuchsaufbau am Palettenwechsler einer LFM. Die Beleuchtung ist 
ein breitbandiger Strahler, weshalb eine deutliche Rotfärbung zu sehen ist. 


men mit dem Schmalbandfilter, aber ohne NIR-Beleuchtung sichergestellt, 
dass NIR-Licht nur in unwesentlichem Maße in den Räumen vorhanden 
ist. Dies hätte bedeutet, dass die Fremdlichunterdrückung nur zum Teil 
funktioniert und hätte entsprechend die Ergebnisse beeinflussen können. 
Bei Aufnahmen mit der NIR-Beleuchtung können die Stege klar erkannt 
werden (siehe Abbildung 2.21). 

Um verschiedene Bedingungen abbilden zu können, wurden Stege mit 
verschiedenen Materialien, nämlich Baustahl, Edelstahl und Kupfer, ge- 
testet. Außerdem waren manche Stege neu, während andere nur gealtert 
waren, was einen Einfluss auf die Oberfläche der Stege hat, oder starke 
Verschlackung aufwiesen. Bei letzteren waren auch geschmolzene Spitzen 
vorhanden und einzelne Spitzen wurden absichtlich verbogen, um das 
vorgeschlagene Messsystem auch in der Detektion einzelner Spitzen testen 
zu können. Insgesamt wurden 20 unterschiedliche Bilder mit jeweils 12 
Stegsteckpositionen im Sichtfeld der Kamera aufgenommen. 
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2.5.5 Ergebnisse und Diskussion 


Die Ergebnisse der Auswertungen wurden manuell überprüft. Von den 
insgesamt 240 aufgenommenen Stegsteckpositionen waren 150 belegt und 
90 nicht belegt. Da ein leichtes Klassenungleichgewicht vorliegt, soll neben 
der Accuracy (siehe Formel 2.3) auch die Balanced Accuracy (siehe Formel 
2.4) betrachtet werden. 


RP I RN 
2 RP+EN’ 2° FP+RN 


Balanced Accuracy = (2.4) 

Ohne Korrektur von Punktreflexionen werden alle 150 eingelegten Stege 
erkannt, zusätzlich aber 31 falsch positiv klassifiziert. Die Accuracy beträgt 
somit 87,1%, die Balanced Accuracy 91,4%. Bei Korrektur der Punktrefle- 
xionen gibt es keine Fehlklassifikationen. 


Abbildung 2.25 Zwei Ausschnitte mit teilweise korrekt (grün eingekreist) und teilweise 
falsch erkannten Spitzen. Links: Ein Beispiel einer falsch positiv erkannten Spitze (rot ein- 
gekreist). Rechts: Ein Beispiel für zwei falsch negativ (orange eingekreist) und eine falsch 
positiv erkannte Spitzen (rot eingekreist). 


Wenn auf Basis dieser korrigierten Stegerkennung außerdem jede ein- 
zelne Stegspitze erkannt werden soll, kommt es, wie zu erwarten war, zu 
leicht schlechteren Ergebnissen. Da die Stege nicht symmetrisch sind, sind 
von einem Steg entweder 45 oder 46 Spitzen im Sichtfeld der Kamera. 
Insgesamt wurden 6870 mögliche Stegspitzen aufgenommen, davon waren 
6616 tatsächlich vorhanden. Die erreichte Accuracy beträgt 98,7% und 
die Balanced Accuracy 93,1%. In diesem Fall ist die Balanced Accuracy 
aussagekräftiger, da das Klassenungleichgewicht recht stark ist. 
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Es treten ungefähr doppelt so viele falsch positive Erkennungen auf wie 
falsch negative. Die falsch positiven Erkennungen kommen sehr oft da- 
durch zustande, dass eine halb abgeschmolzene Spitze trotzdem sehr hell 
erscheint, teilweise auch, weil der Hintergrund sehr dunkel ist. Ein Beispiel 
dafür findet sich in Abbildung 2.25 links. Man könnte diese Fehldetektion 
mit den bestehenden Algorithmen lösen, wenn man die Höhe des zur 
Detektion der Spitzen benutzten Teilbildes verringern würde. Dadurch 
wird aber das aus der Parallelprojektion entstehende Signal verrauschter. 

Falsch negative Erkennungen treten vor allem dann auf, wenn die Umge- 
bung der Spitze verrauscht ist. So kam es teilweise durch einen der beiden 
Querträger der Palette (siehe auch Abbildung 2.15) zu Fehlerkennungen, 
da dieser dann beleuchtet wird, wenn mehrere Stege fehlen. Ein Beispiel 
ist in Abbildung 2.25 rechts zu sehen. Da einzelne Spitzen einen Mindest- 
abstand haben müssen, werden die beiden Spitzen rechts und links des 
Querträgers nicht erkannt, weil dieser aufgrund der Helligkeit als Spitze 
gewertet wird. 
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In vielen Bereichen der Produktion und Logistik entstehen Optimierungs- 
probleme, bei denen mehrere kleinere Gegenstände in einen größeren 
Gegenstand gepackt werden müssen. Beispiele sind die Beladung eines 
Lastkraftwagens mit Paketen unterschiedlicher Größe oder die Aufteilung 
eines großen Ausgangsmaterials in kleinere Stücke durch gerade Schnitte 
in der Papier- und Glasindustrie. Auch bei der Blechbearbeitung auf einer 
LFM müssen die zu schneidenden Teile möglichst effizient angeordnet 
werden. 

Während bei diesen Problemen immer versucht wird, möglichst wenig 
Platz ungenutzt zu lassen, müssen alle kleineren Gegenstände in den grö- 
ßeren Gegenstand passen und dürfen sich nicht überlappen. Es handelt 
sich also um Optimierungsprobleme mit geometrischen Bedingungen. 
In der englischsprachigen Literatur werden diese Probleme oftmals als 
Cutting- and Packing Problems zusammengefasst, sollen im Folgenden aber 
als Schachtelung bezeichnet werden. Während die Probleme natürlich i. A. 
als n-dimensionale Probleme definiert werden können, soll in dieser Arbeit 
nur der 2D-Fall betrachtet werden, der in der Blechbearbeitung auftritt. Die 
kleineren Gegenstände werden deshalb als Teile bezeichnet, der größere 
Gegenstand als Tafel oder Blech. 


3.1 Stand der Technik 


3.1.1 Schachtelung 


Nach der in der Literatur gängigen Klassifizierung sind folgende Kriterien 
relevant, um Schachtelungsprobleme voneinander zu unterscheiden [81]: 
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= Output-Maximierung oder Input-Minimierung: Von Output-Maxi- 
mierung wird gesprochen, wenn ein Problem mehr Teile hat als mög- 
licherweise auf die Tafel passen. Es soll möglichst viel Fläche genutzt 
werden, wozu nicht nur die geometrische Anordnung der Teile be- 
stimmt werden muss, sondern auch ausgewählt werden muss, wel- 
che Teile überhaupt berücksichtigt werden. Bei Input-Minimierungs- 
Problemen muss diese Auswahl nicht getroffen werden, denn die 
gegebene Teilemenge muss vollständig platziert werden. Dazu soll 
möglichst wenig Rohmaterial verwendet werden. 


Beschränkungen des Rohmaterials: Das Rohmaterial kann entweder 
in beiden Dimensionen eine definierte Ausdehnung haben oder wird 
in einer Dimension als unendlich lang angenommen. Dies ist eine 
Annäherung für Rohstoffe, die auf einer Rolle angeliefert werden. 


Anzahl der Teiletypen: Wie viele geometrisch verschiedene Teilety- 
pen verarbeitet werden sollen, hat einen Einfluss auf die Komplexität 
des Problems. Wenn im Extremfall nur ein Teiletyp verwendet wird, 
sinkt der Aufwand für geometrische Berechnungen erheblich. Ent- 
sprechend werden Cutting Stock Problems mit wenigen Teiletypen 
und Bin Packing Problems mit vielen Teiletypen unterschieden. 


In der Blechfertigung soll immer eine gegebene Teilemenge hergestellt 
werden, es handelt sich also immer um Input-Minimierungs-Probleme. 
Dabei kann es sowohl vorkommen, dass von jedem Teiletyp nur ein Teil 
zu platzieren ist, als auch, dass nur ein einziger Teiletyp benutzt wird. 

Bleche werden oft als rechteckige Tafeln hergestellt, sind also in beiden 
Dimensionen beschränkt. Bei dieser Art Problem muss aber nicht nur das 
Layout festgelegt werden, sondern auch entschieden werden, wie viele Ble- 
che benutzt werden müssen und welches Teil auf welchem Blech platziert 
wird. Diese Aspekte sollen für diese Arbeit aber keine Rolle spielen, weil 
sie sich durch die später vorgestellten Nebenbedingungen, die spezifisch 
für eine LFM sind, nicht ändern. Zur Vereinfachung der Darstellung soll 
deshalb angenommen werden, dass die Breite des Blechs definiert ist, die 
Länge hingegen unbeschränkt. Es wird die kürzeste Länge gesucht, in die 
eine gegebene Teilemenge geschachtelt werden kann. Diese Art Problem 
wird in der Literatur oft Strip Packing genannt [81]. 
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én бй 


Abbildung 3.1 Verschiedene Bewertungskriterien. Links: Länge. Mitte: Konvexe Hülle. 
Rechts: Projektion. Eigene Darstellung nach [27]. 


Das Optimierungsproblem wird außerdem ungleich schwieriger, wenn 
jedes zu platzierende Teil in einem frei wählbaren Winkel platziert wer- 
den kann. Eine Zwischenstufe sind wenige, für jeden Teiletyp definierte 
Winkellagen. Oftmals sind dies in der Blechfertigung Drehungen um 90° 
von einer definierten Ausgangslage. 

Eine wichtige Nebenbedingung für das Laserschneiden ist der minimale 
Abstand zwischen zwei Teilen. Dieser ist aus thermischen Gründen not- 
wendig. Um die Darstellung der Probleme nicht unnötig komplizierter zu 
machen, wird im folgenden angenommen, dass der Minimalabstand zwi- 
schen einem Teil und dem Tafelrand und der Minimalabstand zwischen 
zwei Teilen identisch sind. 

Eine nicht unerhebliche Rolle bei der Bewertung eines Schachtelungs- 
algorithmus fällt dem Bewertungskriterium zu, mit dem die Güte einer 
Schachtelung bestimmt wird. Welcher Teil des ungenutzten Materials dem 
Verschnitt oder dem wiederverwendbaren Material zugeordnet werden 
soll, ist keine triviale Frage und hängt maßgeblich vom Anwendungsfall 
ab. Dazu zählt die einfache Länge, die vor allem im Strip Packing als Krite- 
rium genutzt wird (vgl. [16, 77]). Andere Möglichkeiten sind die konvexe 
Hülle oder eine Projektion [27], wie in Abbildung 3.1 dargestellt. Weitere 
mögliche Bewertungskriterien sind z. B. in [27] zusammengestellt. 


3.1.1.1 Modellierung 


Schachtelungsprobleme können auf unterschiedliche Art formalisiert wer- 
den. In allen Darstellungen wird dabei versucht, eine möglichst effiziente 
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Anordnung der Teile zu finden und dabei drei Nebenbedingungen einzu- 
halten: 


= Alle zu platzierenden Teile werden platziert. 
= Für jedes mögliche Teilepaar gilt, dass sie sich nicht überlappen. 
= Alle Teile liegen vollständig auf dem Ausgangsmaterial. 


Ein großer Unterschied zwischen verschiedenen Modellierungsarten 
ergibt sich daraus, dass das Blech und die Teile entweder diskret [77] oder 
kontinuierlich [2] dargestellt werden können. Zur Darstellung des Modells 
in Kapitel 3.2 soll eine diskrete Formulierung verwendet werden, da diese 
leichter zu formalisieren und zu implementieren ist. Deshalb soll der Fokus 
dieses Kapitels auch auf diskreten Modellen liegen. 

Das erste diskrete Modell, das eine exakte Lösung des Schachtelungs- 
problems mit Teilen, die eine beliebige Kontur haben, ermöglicht hat, 
wurde in [77] veröffentlicht. Hier wird der größere Gegenstand mit Punk- 
ten diskretisiert, die auf einem rechtwinkeligen Raster entlang der beiden 
Dimensionen liegen. Daher stammt auch der Name Dotted Board Model. 
Die Teiletypen werden jeweils mit einem Referenzpunkt an der Außenkon- 
tur versehen. Dieser kann zu Beginn willkürlich gewählt werden, ändert 
sich aber nicht. Um sicherzustellen, dass die Teile vollständig innerhalb 
des größeren Gegenstands liegen und sich nicht überlappen, werden zwei 
wichtige geometrische Konzepte benutzt, nämlich das Inner Fit Polygon 
(ТЕР) und das № Fit Polygon (МЕР) [5, 6]. Die Idee bei der Konstruktion die- 
ser beiden Polygone ist ähnlich. Das IFP beschreibt die Fläche des größeren 
Gegenstands, in der sich der Referenzpunkt eines Teiletyps befinden muss, 
damit dieses Teil vollständig im größeren Gegenstand liegt (siehe Abbil- 
dung 3.2). Die Berechnung, ob ein Teil nach seiner Platzierung vollständig 
im größeren Gegenstand liegt, lässt sich mit dem IFP auf die Berechnung, 
ob der Referenzpunkt im IFP liegt, vereinfachen. 

Das NFP ermöglicht diese vereinfachte Berechnung bei der Berechnung 
einer möglichen Überlappung von zwei Teiletypen. Dazu wird ein Teile- 
typ fest platziert und es wird das Polygon errechnet, das entsteht wenn 
der zweite Teiletyp in seiner fixen Rotationslage ohne Überlappung mög- 
lichst dicht am ersten Teiletyp um diesen herum geführt wird. Das NFP 
beschreibt also das Polygon, in dem sich der Referenzpunkt des zweiten 
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Abbildung 3.2 Links: Das NFP (rot) eines an zwei verschiedenen Positionen dargestellten 
Dreiecks mit Referenzpunkt (gelb) oben, bei einem feststehenden Quadrat. Rechts: Das IFP 
(rot) eines Quadrats mit Referenzpunkt (gelb) oben links bei einer rechteckigen Tafel. 


Teiletyps nicht befinden darf, damit sich die Teiletypen nicht überlappen 
(siehe Abbildung 3.2). 

Das Modell benutzt dann binäre Variablen, um die Lage eines Teiletyps 
auf dem Blech zu codieren. Dabei gibt es für jeden Punkt auf dem Blech so 
viele binäre Variablen wie Teiletypen. Die binäre Variable an einem Punkt 
ist für einen Teiletyp dann 1, wenn ein platziertes Teil dieses Teiletyps dort 
seinen Referenzpunkt hat. Sonst ist sie 0. Die Nebenbedingungen können 
nun über diese binären Variablen definiert werden. So müssen mindestens 
so viele binäre Variablen eines Teiletyps 1 sein, wie Teile von diesem Teile- 
typ platziert werden sollen. Außerdem müssen die binären Variablen, die 
1 sind, im IFP dieses Teiletyps liegen. Für jedes platzierte Teil muss gelten, 
dass der Referenzpunkt nicht in dem NFP eines anderen platzierten Teils 
liegt. Durch diese beiden Einschränkungen wird sichergestellt, dass die 
Teile vollständig auf der Tafel liegen und sich nicht überlappen. 

Eine weitere Möglichkeit das Schachtelungsproblem während einer Be- 
rechnung darzustellen, ist die Rasterung des Blechs und der Teile in kleine 
Elemente, die den Pixeln eines Bildes sehr ähnlich sind [6]. Dabei kann 
entweder eine binäre oder eine ganzzahlige Darstellung sinnvoll genutzt 
werden. Bei der binären Darstellung werden ausschließlich 0 und 1 zur 
Codierung des Blechs genutzt. Oftmals wird durch 0 dabei als nicht von 
einem Teil belegte Fläche definiert und 1 entsprechend als belegte Fläche. 
Ein Nachteil dieser Darstellung ist, dass sich die Grenzen zwischen zwei 
Teilen nicht mehr eindeutig feststellen lassen, wenn sich diese berühren 
(siehe Abbildung 3.3). Abhilfe schafft hier eine ganzzahlige Darstellung, 
bei der also auch Zahlen außer 0 und 1 zulässig sind. So kann z.B. 1 be- 
nutzt werden, um die Außenkontur eines Teils zu codieren, während 3 das 
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Abbildung 3.3 Ein Beispiel für die binäre Darstellung einer Schachtelung durch Rasterung. 
Links: Das Layout mit drei Teilen auf der Tafel. Mitte: Eine binäre Darstellung derselben 
Schachtelung. Rechts: Eine Darstellung derselben Schachtelung mit ganzen Zahlen, um den 
Rand eines Teils gesondert zu codieren. 


innere des Teils darstellt (siehe Abbildung 3.3, rechts). Die Teile lassen sich 
nun eindeutig rekonstruieren. 

Eine kontinuierliche, also geometrisch exakte Formulierung, kann durch 
trigonometrische Berechnungen erreicht werden. Es kann aus der Beschrei- 
bung der Teile als Polygone bestimmt werden, ob diese bei einer gegebenen 
Platzierung überlappen oder nur unvollständig auf dem Blech liegen [16]. 
Der große Vorteil dieser Modellierung ist, dass man sich nicht für ein Ras- 
terungsmaß entscheiden muss und damit keinen Quantisierungsfehler 
einführt, der je nach Geometrie der Teile und des Blechs einen mehr oder 
minder großen Einfluss haben kann [16]. 


3.1.2 Exakte Lösung von Schachtelungsproblemen 


Erste Veröffentlichungen, die exakte Lösungen für ein Schachtelungspro- 
blem betrachtet haben, sind oftmals von starken Einschränkungen ausge- 
gangen, wie 2. В. nur einem einzigen Teiletyp [58, 67] oder ausschließlich 
rechteckigen Teilen [62]. Durch diese Einschränkungen werden die geo- 
metrischen Nebenbedingungen deutlich einfacher darstellbar. 

Das erste veröffentlichte Modell, das Teile mit beliebiger Form zur exak- 
ten Lösung optimieren konnte, war das im vorherigen Abschnitt darge- 
stellte Dotted Board Model. Die Modellierung erfolgt hier bewusst so, dass 
ein gemischt ganzzahliges Problem (MIP, für eng. Mixed Integer Problem) 
entsteht. Bei einem gemischt ganzzahligen Problem muss mindestens eine 
Variable ganzzahlig sein, während dies für mindestens eine weitere Varia- 
ble nicht gilt. Für diese Probleme lassen sich durch iteratives Annähern 
exakte Lösungen finden, z. В. durch das Branch and Bound Verfahren [51]. 
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Hierbei wird der Suchraum des Optimierungsproblems in jeder Iteration 
in kleinere Suchräume unterteilt (engl. branching genannt). Dieser Schritt 
alleine würde dazu führen, dass alle möglichen Lösungen im Sinne eines 
brute force Ansatzes berechnet werden. Stattdessen wird aber der noch zu 
durchsuchende Lösungsraum immer weiter eingeschränkt, sodass ganze 
Teilräume nicht durchsucht werden müssen, weil die optimale Lösung 
nicht in diesem Teilraum liegen kann (engl. bounding genannt). Dafür wer- 
den oft obere und untere Schranken benutzt. Beim Strip-Packing kann 
eine untere Schranke aus der Teileflache hergeleitet werden. Die Länge des 
Strips kann nicht kleiner sein, als eine Schachtelung ohne Verschnitt wäre, 
denn das ginge nur mit Überlappung der Teile. Eine obere Schranke für 
das Strip-Packing lässt sich z. В. aus der Summe der längeren Dimensi- 
on aller Teile berechnen. Dann würden alle Teile nacheinander und eben 
nicht nebeneinander liegen, was die schlechteste zulässige Schachtelung 
darstellt, die nicht offensichtlich verbessert werden kann. 

Auch andere Modellierungsansätze, ob als Rasterung oder exakte Geo- 
metrie, nutzen die Formulierung als MIP, um eine exakte Lösung zu er- 
reichen [2, 15, 16]. Die dabei eingesetzten exakten Algorithmen wurden 
auf Problemgrößen mit bis zu 16 Teilen angewandt und benötigen dafür 
eine Rechenzeit zwischen ungefähr 2,5 und 5,5 Stunden. Probleme mit 6 
Teilen brauchen im Vergleich deutlich weniger, nämlich zwischen 5 und 
30 Sekunden. 

Die mit der Teileanzahl stark ansteigenden Rechenzeiten sind nicht ver- 
wunderlich, da Schachtelungsprobleme NP-vollständig sind [26]. Die Zeit, 
die benötigt wird, um eine exakte Lösung zu berechnen, steigt also im 
schlimmsten Fall vermutlich stärker als polynomiell mit der Größe der Ein- 
gangsvariablen, hier also vor allem der Anzahl der Teile. Deswegen wurde 
in der Vergangenheit neben der Beschreibung exakter Lösungsalgorith- 
men vor allem versucht, schnelle Heuristiken mit guter Lösungsqualität zu 
finden. Heuristiken sind approximierende Verfahren, die nicht unbedingt 
die optimale Lösung finden. 


3.1.3 Heuristische Lösung von Schachtelungsproblemen 


Für die in der Industrie auftretenden Schachtelungsprobleme mit bis zu 
mehreren Tausend Teilen sind exakte Lösungsverfahren zu langsam. Des- 
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wegen wird oftmals auf verschiedene Heuristiken zurückgegriffen und 
Ergebnisse von Hand kontrolliert und ggf. korrigiert. 

Die erste publizierte heuristische Lösung für das Schachtelungsproblem 
ist die sogenannte Bottom-Left-Heuristik [3]. Dabei werden Teile in einer 
gegebenen Reihenfolge nacheinander platziert. Es wird eine Haupt- und 
eine Nebenrichtung definiert, oftmals ist „nach unten” (engl. bottom) die 
Hauptrichtung und „nach links” (engl. left) die Nebenrichtung, woher 
der Name des Verfahrens stammt. Die Position eines Teiles wird durch 
Verschiebe-Operationen bestimmt. Das Teil wird zunächst in der Ecke 
platziert, die der Haupt- und Nebenrichtung gegenüber liegt. Von dort 
aus wird das Teil in jedem Schritt so weit wie möglich in Hauptrichtung 
verschoben. Wenn dies nicht möglich ist, wird das Teil so weit wie möglich 
in Nebenrichtung verschoben. Wenn kein Verschieben des Teils möglich 
ist, wird mit dem nächsten Teil fortgefahren. 

Die Bottom-Left-Heuristik stellt eine geometrische Anordnungsvorschrift 
für eine gegebene Reihenfolge der Teile dar. Zur Optimierung der Schach- 
telung muss nur die Reihenfolge der Teile optimiert werden. Dies ist mit 
Meta-Heuristischen Ansätzen möglich, wie z.B. einem evolutionären Al- 
gorithmus [25]. 

Evolutionäre Algorithmen gehören zu den stochastischen Optimierungs- 
verfahren. Der Name stammt aus der schemenhaften Anlehnung an die 
Grundprinzipien der Evolutionslehre. Es werden in vielen Publikationen 
deshalb auch die Begriffe der Evolutionslehre verwendet, um Verfahren 
der Optimierungstheorie zu beschreiben [74]. Um Unklarheiten zu vermei- 
den, wird im Folgenden so weit wie möglich darauf verzichtet. Bei einem 
evolutionären Algorithmus gibt es eine Menge an möglichen Lösungen, 
die zu einem Zeitpunkt untersucht wird. Die Anzahl der Lösungen pro 
Iteration ist vorab definiert. Die initialen Lösungen werden oftmals zufällig 
generiert. Je nach Problemgröße und Vorwissen können auch mit anderen 
Verfahren initiale Lösungen generiert werden. Bei der Schachtelung von 
rechteckigen Teilen bietet es sich beispielsweise an, die Teile nach der Größe 
zu sortieren [48]. Die Lösungen werden dann anhand einer Bewertungs- 
funktion bewertet und können so in eine eindeutige Reihenfolge gebracht 
werden. 

Der Übergang von einer Iteration zur nächsten passiert in zwei Stufen 
[22]. Erstens werden bestehende Lösungen durch ein Selektionsverfahren 
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ausgewählt, damit aus diesen Lösungen neue Lösungen generiert werden 
können. Als Selektionsverfahren kann eine zur Bewertung proportionale 
Selektion gewählt werden. Hierbei ist die Wahrscheinlichkeit, dass eine 
Lösung gewählt wird, von der absoluten Bewertung dieser Lösung im 
Verhältnis zur Summe aller absoluten Bewertungen abhängig. Der große 
Nachteil daran ist, dass sich gute Lösungen sehr schnell vollständig durch- 
setzen und dadurch nur ein kleiner Teil des potentiellen Suchraums auf 
gute Lösungen geprüft wird. Ein anderes Selektionsverfahren ist die Tour- 
nament-Selektion. Hierbei wird zuerst eine Teilgruppe aus der Menge der 
Lösungen zufällig ausgewählt. Diese Gruppe an Lösungen wird nach ihrer 
jeweiligen Bewertung sortiert und die beste Lösung ausgewählt. Da alle 
Lösungen mit der gleichen Wahrscheinlichkeit in die Teilgruppe kommen, 
setzen sich einzelne Lösungen nicht ganz so schnell durch [22]. 

Neue Lösungen werden dann aus der oder den gewählten bestehenden 
Lösungen durch verschiedene Veränderungen generiert. Die drei verbrei- 
tetsten sind Elitismus, Mutation und Kreuzung [22]. Elitismus bedeutet 
dabei, dass ein vorher spezifizierter Anteil der besten Lösungen in die 
nächste Iteration übernommen wird. Dies führt dazu, dass einmal gefun- 
dene gute Lösungen nicht wieder verworfen werden. Bei der Mutation wird 
eine Lösung in sich selbst stochastisch verändert. Für die Reihenfolge von 
Teilen bei der Schachtelung bedeutet das, dass eine bestimmte Anzahl an 
Teilen, die in der Reihenfolge benachbart sind, aus dieser Stelle der Reihen- 
folge ausgeschnitten und an einer zufälligen anderen Stelle der Reihenfolge 
eingefügt werden. Bei der Kreuzung werden zwei Lösungen kombiniert. 
Dafür wird die erste Lösung bis zu einer gewissen Stelle übernommen und 
der Rest der neuen Lösung stammt aus der zweiten originalen Lösung. Bei 
den Reihenfolgen von Teilen für die Schachtelung wird also eine Lösung 
bis zu einer zufälligen Stelle übernommen und danach die Teile einer zwei- 
ten Lösung eingefügt. Hierbei kann es zu zwei Fehlern kommen, nämlich, 
dass Teile nicht mehr vorkommen oder doppelt vorhanden sind. Die neu 
entstandene Lösung muss um diesen Effekt durch Einfügen und Löschen 
von Teilen bereinigt werden. Dabei muss definiert werden, welches Teil 
einer möglichen Dopplung gelöscht werden soll. Ebenso braucht es eine 
Regel, wie nicht mehr vorhandene Teile eingefügt werden sollen. [25]. Der 
evolutionäre Algorithmus läuft dann, bis ein definiertes Abbruchkriterium 
erreicht wird. Am häufigsten ist dies eine Anzahl an Iterationen bzw. eine 
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Zeitschranke, stattdessen oder zusätzlich kann aber auch eine bestimmte 
Lösungsqualität angegeben werden. 

Bei den publizierten Ergebnissen wird schnell klar, dass der hier skiz- 
zierte Algorithmus bestehend aus einem evolutionären Optimierer und 
einer Bottom-Left-Anordnungsvorschrift gut geeignet ist, um rechteckige 
Teile zu schachteln [48]. Bei beliebigen Konturen hingegen können schwer- 
wiegende Fehler auftreten, insbesondere wenn ein Bereich der Tafel nicht 
mehr erreicht werden kann, da ein Teil den Zugang zu diesem Bereich blo- 
ckiert. Abhilfe kann hier eine Erweiterung der Bottom-Left-Regel schaffen, 
die nicht erreichbare Bereiche der Tafel nach jeder Platzierung registriert 
und Teile in dem kleinsten Bereich der Tafel platziert, in den dieses Teil 
passt [11]. Die entstandenen Lücken werden also zuerst gefüllt, weshalb 
das Verfahren als Bottom-Left-Fill bekannt ist. 

Andere Meta-Heuristiken, wie z. B. Simulated Annealing [36] oder Particle- 
Swarm Algorithmen [71] sind ebenfalls auf das Schachtelungsproblem an- 
gewandt worden. Ebenso können verschiedene Anordnungsvorschriften, 
wie z.B. Methoden, die auf das NFP zurückgreifen, benutzt werden [35]. 
Diese weiteren Ansätze sollen hier aber nicht genauer betrachtet werden, 
da sie für die weitere Abhandlung nicht relevant sind. 


3.1.4 Schachtelung mit Nebenbedingungen aus dem 
Produktionsprozess 


In der wissenschaftlichen Literatur existiert keine Abhandlung über die 
beschriebene Problematik der Veränderung der Schachtelung zur Stei- 
gerung der Prozesssicherheit oder Verringerung des Verschleißes beim 
Laserschneiden. Im Patent [40] hingegen wird genau diese Problematik 
behandelt. Die dort präsentierte, recht einfache Lösung besteht darin, den 
Abstand zwischen zwei Teilen von vorne herein so groß zu wählen, dass 
ein Rand entsteht, in dem im Nachhinein jedes Teil bestmöglich verschoben 
werden kann. Dieser Ansatz optimiert also nicht den Verschnitt, sondern 
ausschließlich die Prozesssicherheit und nimmt einen deutlich gestiege- 
nen Verschnitt in Kauf, unabhängig davon, ob eine Verschiebung für die 
konkreten Teile überhaupt notwendig gewesen wäre. Durch den hohen 
Anteil der Materialkosten an den Teilekosten ist dieses Verfahren in der 
Praxis unwirtschaftlich. 
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Anders als bei der Schachtelung gibt es zur Planung der Verfahrwege 
des Schneidkopfs einer LFM mit Berücksichtigung des Verkippverhaltens 
von geschnittenen Teilen wissenschaftliche Literatur [17]. Bei Teilen, die 
verkippen könnten, wird dann entweder der Schneidkopf hochgezogen 
oder das Teil vollständig umfahren. Es bleibt allerdings offen, wie sich 
das Verkippen der Teile genau berechnen lässt und wie die Auflagestege 
konkret dazu beitragen. Die Schachtelung wird als gegebener Input der 
Planung der Verfahrwege angesehen und wird deshalb nicht verändert. 

Es wurde schon Anfang der 1990er Jahre erwähnt, dass eine möglichst 
enge Schachtelung der Teile nicht zwangsweise zu guten Ergebnissen in 
industriellen Prozessen führt [80]. Die weiteren Faktoren, die dabei berück- 
sichtigt wurden, waren eine mögliche Überproduktion bei regelmäßigen 
Aufträgen, um weniger Verschnitt zu erzeugen, und die Lagerkosten für 
die überschüssigen Teile. 

Eine weitere bei vielen Fertigungsverfahren zu beachtende Nebenbe- 
dingung ist die Anzahl der Rüstvorgänge, die zwischen zwei verschie- 
denen Produktionsaufträgen nötig werden und die Produktionszeiten 
entscheidend verlängern können. Es werden also möglichst viele gleiche 
Schachtelungen gesucht. In [63] wird dies grundsätzlich betrachtet, für die 
Anwendungsfälle einer Stanzmaschine auch in [38] und die Papierindustrie 
in [29, 47]. Auf Grund der Prozessunsicherheiten bei der LFM wird in der 
Praxis manchmal auch nach möglichst vielen gleichen Schachtelungen ge- 
sucht. Diese werden dann von Hand im Hinblick auf die Prozesssicherheit 
optimiert. 

Eine andere Nebenbedingung aus dem Produktionsprozess, deren Ein- 
fluss auf die Schachtelung untersucht wurde, ist die Walzrichtung des 
Blechs. In [38] wird ein Verfahren vorgestellt, bei dem jedes Teil eines 
Teiletyps in derselben, vordefinierten Rotation relativ zur Walzrichtung 
ausgerichtet wird. Die Rotation der Teile ist also vorgegeben und kann 
nicht durch den Schachtelungsalgorithmus verändert werden. 

Die hier erwähnten Arbeiten sind sich in der Hinsicht ähnlich, dass 
sie nie den Schachtelungsalgorithmus selbst verändern, sondern auf be- 
stehende Algorithmen zurückgreifen. Die im nächsten Kapitel vorgestell- 
ten Verfahren unterscheiden sich in diesem wichtigen Punkt vom Stand 
der Technik. 
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3.2 Schachtelung unter Berücksichtigung der 
Auflagepunkte 


Um die in Kapitel 1.1 angesprochenen Prozesssicherheitsprobleme zu 
lösen, wird nun versucht, die Auswirkungen der Auflagestege auf den 
Produktionsprozess schon während der Schachtelung zu berücksichtigen. 
Dazu wird das Problem zunächst formal modelliert. Anschließend werden 
Lösungen durch zwei Heuristiken betrachtet. Diese beachten zusätzlich 
zum offensichtlichen Ziel einer möglichst hohen Schachtelungsdichte wei- 
tere Nebenbedingungen, nämlich die Stabilität bzw. das Verkippen von 
Teilen in Abhängigkeit vom Freischnittpunkt und die Stegbeschädigung. 

Teile der folgenden Ausführungen, insbesondere zur formalen Problem- 
beschreibung und zur Lösung durch eine Konstruktionsheuristik, sind 
bereits in [87] veröffentlicht worden. 


3.2.1 Formale Problembeschreibung und exakte Lösung 


Es lässt sich in Anlehnung an das Dotted Board Model ein Modell zur Schach- 
telung unter Berücksichtigung der Auflagepunkte definieren. Der wich- 
tigste Unterschied besteht dabei in der Wahl der Zielfunktion, da jetzt 
nicht mehr nur die Packungsdichte (2,4) optimiert wird. Es sollen zwei 
zusätzliche Faktoren berücksichtigt werden, nämlich wie viele Teile kippen 
(z,) und wie stark die Auflagestege beschädigt werden (z,). Um den durch 
die verschiedenen Zielvariablen ausgelösten Zielkonflikt entscheiden zu 
können, sind Gewichte wpa, их und wy für jede Variable nötig, sodass eine 
additive Zielfunktion entsteht (Gleichung 3.1). Eine weitere Möglichkeit 
wäre hier die hierarchische Optimierung, die für die verschiedenen Zielva- 
riablen eine klare Hierarchie festlegt und weniger wichtige Zielvariablen 
erst dann optimiert, wenn die wichtigeren Variablen nicht mehr verbessert 
werden können. Diese Möglichkeit erscheint aber nicht sinnvoll für die vor- 
liegende Anwendung, da der beste Kompromiss unter den Zielvariablen 
gefunden werden soll. 

Für diese Modellierung soll das Blech als Matrix dargestellt werden, 
wobei die Matrixelemente analog zum Pixel eines Bildes zu sehen sind. 
Jedes Element der Matrix lässt sich also über die diskreten Koordinaten 
(и, v) eindeutig adressieren. Die verwendeten Teiletypen t; sind ebenfalls 
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als Matrix ihres einhüllenden Rechtecks dargestellt, wobei $ die jeweili- 
ge Rotation des Teils angibt. Die Form des Teils wird dabei durch binäre 
Einträge dargestellt. Die Teiletypen werden durch einen Referenzpunkt 
1 = (u,v) auf dem Blech platziert. Dazu wird t; auf alle entsprechenden 
(и, v) addiert. Wenn dadurch ein Matrixelement (и, v) größer als 1 werden 
sollte, stellt dies eine unzulässige Überlappung von Teilen dar. Der Refe- 
renzpunkt eines Teiletyps darf nur so gewählt werden, dass das Teil auch 
vollständig auf dem Blech liegt. Zur Überprüfung wird das IFP herange- 
zogen, wobei p(t;) die Menge aller Elemente im ТЕР für einen Teiletyp in 
einer Rotationslage ist. Die zu platzierende Anzahl an Teilen pro Teiletyp 
ist durch q, gegeben. 

Die Auflagestege werden durch eine Matrix Ф dargestellt, deren Einträge 
nur die Werte 0, 1 oder 2 annehmen können. Die Einträge stehen respektive 
für leeren Raum unter dem Blech, einen Steg aber keine Spitze oder eine 
Stegspitze. Die Außenkonturen der Teile und damit die Schneidkonturen 
sind durch die Matrix © gegeben, die für eine Schneidkontur 1 ist und für 
alle anderen Einträge 0. Ein Element dieser Matrizen wird durch einen 
entsprechenden Subskript notiert, z.B. Ө, „. Das Verkippen eines Teils 
ist von seiner Rotation, seinem Referenzpunkt und den Auflagestegen 
abhängig. Dies wird als ®,, (W) notiert, wobei die Einträge 0, 1 und 2 hier 
respektive für ein stabiles oder fallendes, ein unsicheres und ein kippendes 
Teil stehen. 

Die notwendigen Nebenbedingungen des Schachtelungsproblems sind 
im Grundsatz schon angesprochen worden. Es muss sichergestellt werden, 
dass alle Teile platziert werden. Dies wird durch Gleichung 3.5 sicherge- 
stellt. Außerdem dürfen sich Teile nicht überlappen, was durch Gleichung 
3.6 gegeben ist. Damit alle Teile auch vollständig auf dem Blech liegen, 
muss Gleichung 3.7 erfüllt sein. 

In dieser Formulierung des Problems muss außerdem sichergestellt 
werden, dass die Zielvariablen den korrekten Wert zugewiesen bekommen. 
Dabei wird nur eine untere Schranke angegeben. Für 24 bedeutet dies, 
dass die größte Spalte v ausgewählt wird, in der ein Eintrag größer als 1 
ist (Gleichung 3.2). Bei der Stegbeschädigung wird z, als Summe über alle 
Produkte der Außenkonturen und Auflagestege definiert (Gleichung 3.3). 
Das Verkippen z, wird ebenfalls als Summe über die Stabilität aller Teile 
dargestellt (Gleichung 3.4). 
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Das Optimierungsproblem als Ganzes ist in den Gleichungen 3.1 bis 3.10 
dargestellt. Die letzten drei Gleichungen definieren dabei ausschließlich 
den Wertebereich der Variablen. 


min 2ра ` ра + Zk Wk + Zp + Wp (3.1) 
u.d.N. max{v| Wo) аа (3.2) 
Е (u,v) < 2b (3.3) 
a, 2, Yat (8.4) 
> Le >q Vt (3.5) 
5 0) (u,v) < 1 У (u, v) (3.6) 
Т 
u(t) € plti) Ve (3.7) 
2р4, Žk % € N (3.8) 
P uv)? Viu) € 10, 1, 2} (3.9) 
On) € 10, 1} (3.10) 


Die in Kapitel 3.1.2 eingeführte Vorgehensweise zur exakten Lösung 
mittels der Branch and Bound Methode ist auf das oben aufgestellte MIP 
mit den Nebenbedingungen immer noch anwendbar. Eine interessante 
Frage ist, ob das Problem grundsätzlich durch die neuen Nebenbedin- 
gungen schneller oder langsamer zu berechnen ist. Dies ist zum Einen 
davon abhängig, ob die zusätzlichen Nebenbedingungen und unter Um- 
ständen notwendigen Zielfunktionsvariablen den Suchraum des Problems 
vergrößert oder verkleinert haben. Zum anderen haben auch die Struktur 
des Suchraums und entspchend angepasste Suchstrategien einen großen 
Einfluss. Für manche Probleme aus der industriellen Anwendung führen 
Nebenbedingungen zu einem kleineren Suchraum. Dies ist z. B. der Fall, 
wenn beim Schachteln die Rotationsmöglichkeiten von Teilen auf eine be- 
stimmte Rotation fixiert werden. Gleiche Reihenfolgen von Teilen verlieren 
dann komplett die Möglichkeit von verschiedenen Rotationen und stellen 
immer die gleiche Lösung dar. Der Suchraum ist deshalb erheblich kleiner. 
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Die Erweiterung um die Zielvariablen z, und 2р hat den Suchraum offen- 
sichtlich nicht verkleinert, da keine der möglichen Lösungen des einfachen 
Schachtelungsproblems entfallen. Der Suchraum ist durch die zusätzli- 
chen Zielvariablen aber auch nicht zwangsweise größer geworden, da die 
Lösungsmöglichkeiten immer noch dieselben Schachtelungen sind, auch 
wenn sie jetzt anders bewertet werden. Die Rechenzeit nimmt allerdings 
durch die zusätzlichen Rechnungen in der Bewertung auf jeden Fall zu. 

Exakte Lösungen sind schon ohne die Berücksichtigung von Nebenbe- 
dingungen aus dem Produktionsprozess für die industrielle Anwendung 
uninteressant, da bei Teilen mit beliebiger Kontur und Stückzahlen größer 
als ungefähr 20 die Berechnungen zu lange dauern. Deshalb werden Heu- 
ristiken benutzt, die schnell Lösungen mit relativ hoher Qualität ausgeben 
sollen. Im Folgenden werden zwei Heuristiken vorgestellt, eine Konstrukti- 
onsheuristik und eine Verbesserungsheuristik. Eine Konstruktionsheuristik 
errechnet eine zulässige Lösung direkt aus den Eingaben. Eine Verbesse- 
rungsheuristik benutzt eine bestehende zulässige Lösung und verbessert 
diese durch eine Veränderung der Variablen. 


3.2.2 Lösung durch eine Konstruktionsheuristik 


Eine Konstruktionsheuristik zur Lösung des vorgestellten Problems muss 
nicht nur eine gute Schachtelung der Teile erstellen können, sondern vor 
allem auch die prozessrelevanten Nebenbedingungen berücksichtigen. 
Wie schon in der Modellierung gezeigt, sollen das Verkippen von Teilen 
und die Stegbeschädigung zusätzlich in die Optimierung einfließen. 

Die Optimierung der Schachtelung wird durch einen evolutionären Al- 
gorithmus durchgeführt. Eine Lösung besteht dabei aus einer Permutation 
der Teile. Die Teile sind in der gegebenen Reihenfolge auf der Tafel zu plat- 
zieren, wobei die Bottom-Left-Platzierungsregel benutzt wird. Die Fitness 
wird nach Gleichung 3.1 berechnet. Zur Darstellung des Problems werden 
die Teile und das Blech als diskrete Elemente aufgefasst, genau wie in der 
Modellierung in Kapitel 3.2.1. Für jedes Teil, das platziert werden soll, wird 
dann das Verkippverhalten und die Stegbeschädigung berechnet. Dafür 
müssen die Position der Stege und die Position des Blechs relativ zu den 
Stegen bekannt sein. In der Praxis kann eine der in Kapitel 2 vorgestell- 
ten Methoden benutzt werden, um diese Informationen zu erlangen. Die 
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Verkippberechnung hängt von der Geometrie des Teils, insbesondere vom 
Schwerpunkt, den Auflagepunkten, dem Freischnittpunkt, der Gaskraft 
und dem Gewicht des Teils ab. Vereinfacht gesagt ist ein Teil dann stabil, 
wenn es ohne Gasdruck auf den Auflagestegen stabil liegen kann und das 
über die Gaskraft wirkende Drehmoment beim Freischneiden des Teils 
nicht größer ist als das durch die Schwerkraft eingebrachte Drehmoment. 
Der genaue Berechnungsablauf wird in Anhang A dargestellt. Ähnlich wie 
in der Modellierung in Kapitel 3.2.1 wird die Stegbeschädigung aus der 
Schnittlänge über den Auflagestegen berechnet. Dabei ist aus Versuchen 
bekannt, dass die Beschädigung schlimmer ist, wenn über einer Stegspitze 
geschnitten wird, im Vergleich zu einer Flanke oder Senke. Entsprechend 
gibt es einen Gewichtungsfaktor, der dafür sorgt, dass Stegspitzen im 
Vergleich zum sonstigen Steg doppelt gewichtet werden. 

Wenn das Schachtelungsverfahren ansonsten unverändert beendet wer- 
den würde, würde das bereits helfen, das zu Grunde liegende Problem zu 
lösen. Die grundsätzliche Anordnung der Teile auf dem Blech erfolgt zwar 
ausschließlich nach dem Bottom-Left-Prinzip, aber bei zwei in der Schach- 
telungsdichte gleichen oder zumindest ähnlich guten Schachtelungen wird 
jetzt diejenige gewählt, die gegeben die drei Gewichtungsfaktoren das 
beste Ergebnis liefert. Das hier vorgestellte Verfahren geht aber noch einen 
Schritt weiter. 


3.2.3 Lokales Verschiebverfahren und Berechnungen 
3.2.3.1 Lokales Verschiebverfahren 


Wenn ein Teil seine vom Bottom-Left-Prinzip festgelegte Position erreicht 
hat, wird geprüft, ob durch eine Verschiebung des Teils eine Verbesserung 
des Verkippverhaltens oder der Stegbeschädigung erreicht werden könnte. 
Damit dies nicht der Fall ist, muss das Teil stabil gelagert sein und es darf 
nur wenige Millimeter Schnitt über den Auflagestegen geben. Hier kei- 
nen Schnitt über den Auflagestegen zu fordern, würde nur dazu führen, 
dass die meisten Teile in fast jeder Lage eine Verbesserungsmöglichkeit 
aufweisen. Abhängig von der Größe des Teils und dem Gewichtungsfak- 
tor für die Stegbeschädigung wird deshalb ein geringer Schnitt über den 
Auflagestegen zugelassen. 
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Wenn eine Verbesserung möglich wäre, wird zuerst geprüft, in welchem 
Raum der Referenzpunkt unter Berücksichtigung der unterschiedlichen 
Gewichtungsfaktoren verschoben werden sollte. Dadurch soll verhindert 
werden, dass eine Verschiebung geprüft wird, bei der definitiv keine Ver- 
besserung erzielt werden kann. Dies ist der Fall, wenn die Schachtelungs- 
dichte so weit sinkt, dass die anderen beiden Faktoren definitv keine aus- 
reichend große Verbesserung erreichen können. Entscheidend ist an dieser 
Stelle die Definition des Verschnitts, der durch die Verschiebung entsteht, 
denn den tatsächlich am Ende der Schachtelung entstehenden Verschnitt 
kann man zu diesem Zeitpunkt in der Berechnung nicht direkt vorhersa- 
gen. Die Schachtelung vollständig weiter zu berechnen, bis dieser Effekt 
bekannt ist, lässt die Berechnungskomplexität sehr stark steigen, denn für 
jedes zu platzierende Teil müssten Hunderte Schachtelungen berechnet 
werden. Ein solches Vorgehen erscheint also nicht sinnvoll. 

Stattdessen soll der zusätzliche Verschnitt z, lokal am Teil abgeschätzt 
werden, ohne mögliche Abhängigkeiten mit den danach platzierten Teilen 
zu berücksichtigen. Dazu wird die Verschiebung in der jeweiligen Di- 
mension д,,д, und die Länge des einhüllenden Rechtecks b,, b, in dieser 
Dimension genutzt. Bei der Verschiebung in beiden Dimensionen muss 
außerdem die Fläche, die sich aus der Multiplikation der beiden Verschie- 
bungen ergibt, berücksichtigt werden. Zusammenfassend lässt sich also 
schreiben: 

Z5 = Ôx ` ba + y ` by T Ôn * Oy: 


Ein Beispiel für die Fläche z, ist in Abbildung 3.4 gezeigt. 

Der mögliche Suchraum besteht aus einem rechtwinkeligen Dreieck, in 
dem das Teil entgegen der Bottom-Left-Richtungen verschoben wird. Die 
maximale Verschiebung in mm kann in der jeweiligen Dimension durch: 


Umax 


max = (Tk ` Wk F Op Wp) - 7 


ра 
bestimmt werden, wobei о bzw. съ die maximale Verbesserung bei der 
Stabilität bzw. der Stegbeschädigung (nach dem Optimierungsmodell in 
Kapitel 3.2.1 zwischen 0 und 2 definiert), | die Länge des Teils in der jeweils 
anderen Dimension in mm und umax die Fläche einer Spalte der Blechmatrix 
in mm? ist. Damit die Verschiebung möglichst gering ausfällt, werden die 
möglichen Platzierungen anhand der Euklidischen Distanz sortiert und 
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LS 


Abbildung 3.4 Ein Beispiel für den bei der Verschiebung eines Teils veranschlagten Ver- 
schnitt. Links: Die Position des schwarzen Teils bei Anwendung der Bottom-Left-Regel. 
Rechts: Das schwarze Teil in einer verschobenen Position, der Verschnitt ist blau markiert. 


nacheinander geprüft. Sollte es dazu kommen, dass bei einer Position keine 
weitere Verbesserung bei der Stabilität des Teils und der Stegbeschädigung 
möglich ist, wird diese Position gewählt und das Verfahren bricht ab. Jeder 
weitere Punkt würde nur zu einer größeren, ungewollten Verschiebung 
führen. Sonst werden alle Punkte im Suchraum bewertet und die beste 
gefundene Position gewählt. Dabei muss jeweils auch die Überlappung 
mit anderen Teilen geprüft werden, da die Verschiebungen insbesondere 
in y-Richtung zu Überlappungen mit anderen Teilen führen können. 

Für die Stabilität eines Teils ist es sehr wichtig, an welcher Stelle der Kon- 
tur das Teil freigeschnitten wird. Zum Zeitpunkt des Freischnitts liegt der 
volle Gasdruck an und das Teil kann zum ersten Mal verkippen. Wenn der 
Freischnittpunkt ungeschickt liegt, ist der Hebelarm der Gaskraft lang und 
ein Verkippen des Teils wahrscheinlicher. Insbesondere kleinere Teile aus 
dünnen Materialien, die mit einem hohen Gasdruck geschnitten werden, 
sind gefährdet zu verkippen. 

Bei einer klassischen Programmierung der LFM können die Freischnitt- 
punkte nicht so gewählt werden, dass die Teile möglichst wenig verkippen, 
da die Position und der Zustand der Auflagestege nicht bekannt sind. Dies 
wird bei den Berechnungen durch zufällige Freischnittpunkte modelliert. 
Aus allen Punkten der Außenkontur eines Teiles wird dabei ein Punkt 
zufällig ausgewählt, wobei alle Punkte gleich wahrscheinlich sind. 

Im Gegensatz dazu stehen die bewusst unter Berücksichtigung des Ver- 
kippens und der Stegbeschädigung gewählten Freischnittpunkte, die im 
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Folgenden als beste Freischnittpunkte bezeichnet werden. Sie werden da- 
durch berechnet, dass alle Punkte auf der Kontur geprüft werden. Sobald 
ein Freischnittpunkt gefunden ist, der keine Stegbeschädigung verursacht 
und bei dem das Teil stabil liegen bleibt, wird dieser Punkt gewählt und 
das Verfahren bricht ab. Falls kein solcher Punkt gefunden werden kann, 
wird der beste bekannte Punkt gewählt. 

Einen Sonderfall stellen Teile dar, die so klein sind, dass es keine stabile 
Position geben kann. Sie können daran erkannt werden, dass die längere 
Seite des umgebenden Rechtecks des Teils kleiner ist als der Abstand 
zwischen zwei Stegen. Wenn das Teil in der ursprünglich durch das Bottom- 
Left-Verfahren bestimmten Position einen Steg unter sich hat, wird es als 
kippend angesehen. In diesem Fall muss geprüft werden, wie weit das 
Teil verschoben werden muss, damit das Teil durch die Stege fallen kann. 
Erst wenn kein Auflagesteg, also nicht nur keine Spitze, unter dem Teil 
ist, ist dies der Fall. Abhängig von der Distanz, die das Teil verschoben 
werden müsste, und den Gewichtungsfaktoren kann bestimmt werden, 
welche Option zur Platzierung des Teils die bessere ist. Die Definition des 
Verschnitts ist dabei gleich zum vorhergehenden Absatz über Teile, die 
stabil platziert werden sollen. 


3.2.3.2 Berechnungen 


Die Teilemenge für die Berechnungen besteht aus 20 Teiletypen, die teil- 
weise von Kunden stammen und leicht in der Geometrie verändert wurden 
und teilweise von TRUMPF zu Vorführungszwecken benutzt werden. Aus 
diesen 20 Teiletypen werden 100 Schachtelungen zufällig generiert. Zuerst 
wird dafür die Teileanzahl einer Schachtelung zufällig und gleichverteilt 
zwischen 40 und 62 Teilen festgelegt. Dann wird ein Vektor aus ganz- 
zahligen Zufallszahlen mit entsprechender Länge generiert, wobei die 
Zufallszahlen zwischen 1 und 20 gleichverteilt sind, was den 20 Teilety- 
pen entspricht. Die Teiletypen werden so oft in der Schachtelung benutzt, 
wie sie im Zufallsvektor vorkommen. Die Teiletypen als Bilder und die 
Zusammenstellung der 100 Schachtelungen in tabellarischer Form sind in 
Anhang B dargestellt. 

Es wird im Folgenden immer angenommen, dass der erste Steg direkt 
unter dem Rand des Blechs liegt. Der Abstand zwischen zwei Spitzen 
beträgt 15mm, der Abstand zwischen zwei Stegen 67mm. Das für die 
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Verkippberechnung angenommene Material ist 2 mm dicker Baustahl mit 
einer Dichte von 7856 kg/m, der angenommene Gasdruck beträgt 15 bar 
bei einem Düsendurchmesser von 2mm. 

Um Auswirkungen von verschiedenen Gewichtungsfaktoren zu testen, 
werden verschiedene Varianten gerechnet. In einer ersten Variante wird 
ausschließlich die Bottom-Left-Regel mit zufälligen Freischnittpunkten 
angewandt und es kommt zu keiner weiteren Verschiebung. Die zweite 
Variante unterscheidet sich nur darin, dass die besten Freischnittpunkte 
genutzt werden. 

Um die lokalen Verschiebungen von Teilen zu demonstrieren, werden 
dann mehrere Varianten mit verschiedenen Gewichtungsfaktoren berech- 
net. Bei der ersten Variante hat die Packungsdichte weiterhin die höchste 
Priorität mit wpa = 5, ик = 1, шь = 1. Neben einer Variante, bei der die 
Gewichtungsfaktoren gleichberechtigt sind, gibt es zwei Varianten, bei 
denen die Stabilität der Teile wichtiger als die anderen beiden Faktoren ist 
ра = 1, Wk = 5, иь = 1 bzw. wpa = 1, wk = 10, w, = 1. In der letzten 
Variante wird die Stegbeschädigung als wichtigstes Kriterium gesehen mit 
Фра = 1, Wk = 1, Wp = 5. 

Bei der Optimierung mit dem evolutionären Algorithmus ergibt sich 
das Problem, dass die Konvergenzzeiten nicht nur für unterschiedliche 
Schachtelungen variieren, sondern auch bei mehreren Berechnungsdurch- 
läufen derselben Schachtelung. Der Vergleich zwischen den verschiedenen 
Parametrierungen kann durch diesen Effekt gestört werden. Deshalb wird 
zum Vergleich der verschiedenen Parametrierungen als Startreihenfolge 
eine Vorsortierung verwendet, wobei die Teile zuerst nach der längsten 
Seite aufsteigend und nachrangig nach der kürzeren Seite aufsteigend 
sortiert werden. Die Teile werden dann immer mit der längeren Seite quer 
zu den Auflagestegen platziert. 

Trotzdem soll zusätzlich gezeigt werden, welche Ergebnisse der vor- 
gestellte Optimierungsansatz erzielen kann. Dafür wurde der Elitismus- 
Faktor so gewählt, das immer nur die beste bekannte Lösung in die nächs- 
te Generation übernommen wird. Die Kreuzungs- und Mutationswahr- 
scheinlichkeiten, die sowohl auf die Reihenfolge der Teile als auch ihre 
Winkellage wirken, wurden ebenfalls vor der Berechnung fixiert. Als mög- 
liche Winkellagen war für jedes Teil 0° und 90° zulässig. In jeder Generation 
gab es 20 Lösungen, wovon jeweils 5 an der Tournament-Selektion teilneh- 
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men. Als Abbruchkriterium wurden 50 Iterationen gewählt, da die Angabe 
einer Lösungsqualität schwierig ist und eine Zeitschranke zu unterschied- 
lich vielen Optimierungsschritten führen könnte. 


3.2.4 Ergebnisse und Diskussion der 
Konstruktionsheuristik 


Die ersten drei der 100 Schachtelungen unter Anwendung der Bottom- 
Left-Methode ohne Verschiebungen mit zufälligem Freischnittpunkt sind 
in Abbildung 3.5 dargestellt. Es ist sofort ersichtlich, dass die Schachtelung 
mit dem Bottom-Left-Ansatz bei der benutzten Teilemenge zu Resultaten 
mit geringer Packungsdichte führt, wenn ein gerader Trennschnitt als 
Bewertungskriterium herangezogen wird. Dies liegt daran, dass die Teile 
beliebige Polygone sind und nicht nur Rechtecke oder Teile, deren konvexe 
Hülle ein Rechteck stark annähert. Die Teile haben dadurch eine größere 
Chance, an kleinen Ausbuchtungen der schon platzierten Teile hängen zu 
bleiben, weil sie nicht zurück nach rechts verschoben werden können, um 
dann nach unten verschoben zu werden. Deshalb werden die Teile eher in 
der oberen Hälfte des Blechs platziert, was natürlich zu einer schlechteren 
Packungsdichte führt. Der Effekt kann in den Abbildungen 3.5 und 3.6 
beobachtet werden. 

Veränderungen bei den Gewichtungen der Zielvariablen haben auch 
einen deutlichen Einfluss auf die Ergebnisse. Während bei dem Verfahren 
ohne Verschiebung durchschnittlich nur 15,8 Teile stabil sind, steigt diese 
Zahl deutlich, wenn Verschiebungen möglich sind, auf durchschnittlich 
bis zu 42,7 Teile. Die Anzahl der verkippenden Teile nimmt entsprechend 
ab. Die Anzahl der fallenden Teile nimmt ebenfalls leicht zu, wenn Ver- 
schiebungen möglich sind. Dies ist aber fast unabhängig von den exakten 
Gewichtungen. Die Anzahl der Teile, bei denen das Ergebnis der Verkipp- 
berechnung unsicher ist, nimmt bei geringerer Gewichtung der Fläche 
leicht ab, sinkt aber nicht unter im Durchschnitt 2 Teile. Die Länge des 
benötigten Blechs steigt wie erwartet mit geringerer Gewichtung der Pa- 
ckungsdichte an. Allerdings handelt es sich dabei im Durchschnitt um 
Unterschiede im einstelligen Millimeterbereich. Nur die beiden Gewich- 
tungen, die die Stabilität der Teile sehr wichtig nehmen, benötigen deutlich 
längere Bleche. Die Ergebnisse der verschiedenen Parametrierungen sind 
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Abbildung 3.5 Beispiele von Schachtelungen ohne Verschiebung mit der Konstruktionsheu- 
ristik. Grüne Teile sind stabil, rote Teile kippen, orangene Teile sind unsicher und blaue Teile 
fallen. Die x-Achse verläuft horizontal, die y-Achse vertikal. 
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Abbildung 3.6 Vergleich der Ergebnisse der Konstruktionsheuristik bei verschiedenen Се- 
wichtungsfaktoren und zufälligem Freischnittpunkt. Oben links: ира = 1, wk = 1, wp = 1. 
Oben rechts: wpaq = 1, шу = 5, шь = 1. Unten links: wpa = 1, wg = 1, wp = 5. Unten rechts: 
ра = 1, wg = 10, wp = 
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Abbildung 3.7 Vergleich der Ergebnisse bei zufälligem Freischnittpunkt bei verschieden 
parametrierten Gewichtungsfaktoren. BL bezeichnet dabei Bottom-Left Schachtelung ohne 
Verschiebung, ansonsten sind die Gewichtungen in der Reihenfolge wpaq, Wk, Wp angegeben. 


in Abbildung 3.7 zu sehen. Ein konkretes Beispiel einer Schachtelung mit 
unterschiedlichen Parametrierungen ist in Abbildung 3.6 dargestellt. Ein 
Grund, warum es selbst bei der Gewichtung wpa = 1,u, = 10,0 = 1 
zu kippenden Teilen kommt, ist, dass die Teile nicht notwendigerweise in 
y-Richtung verschoben werden können, obwohl der lokale Suchraum groß 
genug dafür wäre. Zum Teil liegt das daran, dass es sonst zu Überlappun- 
gen mit bereits platzierten Teilen kommen würde, oder aber daran, dass 
der Tafelrand eine weitere Verschiebung verhindert. 

Die Wahl des besten Freischnittpunkts führt zu deutlich mehr stabilen 
Teilen. Ohne Verschiebung und damit bei gleicher Packungsdichte sind 
nun 31,6 Teile stabil und damit fast doppelt so viele wie bei einem zufälligen 
Freischnittpunkt. Der Unterschied bei den stabilen Teilen zwischen dem 
zufälligen und optimalen Freischnittpunkt sinkt deutlich ab, wenn die 
Packungsdichte eher gering gewichtet wird. Der deutlichste Unterschied 
besteht bei den Teilen, deren Stabilität unsicher ist. Der Wert sinkt für alle 
Gewichtungen auf ungefähr 0,1 Teile ab. Wenig überraschend bleibt die 
Zahl der fallenden Teile im Grunde gleich. Die Zahl der kippenden Teile ist 
für Gewichtungen mit Fokus auf der Stabilität der Teile im Grunde genau so 
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Abbildung 3.8 Vergleich der Ergebnisse bei bestem Freischnittpunkt und verschieden ра- 
rametrierten Gewichtungsfaktoren. BL bezeichnet dabei Bottom-Left Schachtelung ohne 
Verschiebung, ansonsten sind die Gewichtungen in der Reihenfolge wpa, Wk, шь angegeben. 


groß wie bei einem zufälligen Freischnittpunkt. Dies spricht dafür, dass die 
Teile nicht durch den Gasdruck kippen, sondern schon ohne Betrachtung 
des Gasdrucks nicht stabil auf den Auflagestegen liegen. Entsprechend 
hat ein anderer Freischnittpunkt keine Auswirkung. Die Ergebnisse der 
verschiedenen Parametrierungen mit optimalem Freischnittpunkt sind in 
Abbildung 3.8 zusammengefasst. 

Die Stegbeschädigung ist ausschließlich von der gewählten Gewichtung 
abhängig und nicht davon, ob ein optimaler oder zufälliger Freischnitt- 
punkt gewählt wird. Die durchschnittliche Schnittlänge über dem Steg 
liegt zwischen 3880 mm und 3920 mm bei den Gewichtungen, die keinen 
großen Wert auf die Stegbeschädigung legen. 

Bei den Schachtelungen mit der Gewichtung wpa = 1, wp = 1,w, = 5 
verbessern sich diese Werte deutlich. Es kommt nur noch zu 2491 mm 
Schnittlänge über dem Steg. 

Die Optimierung mit dem evolutionären Algorithmus ist in der Praxis 
durch hohe Rechenzeiten erschwert. Die Ergebnisse einer Schachtelung 
liegen bei den erwähnten Optimierungsparametern nach ungefähr 4 bis 7,5 
Stunden vor. Deshalb wurden nur Ergebnisse für die ersten zehn Schach- 
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telungen berechnet. Das Problem bei der Optimierung mit dem evolutio- 
nären Algorithmus besteht darin, dass manchmal sehr gute Ergebnisse 
gefunden werden und manchmal nicht, was am recht großen Einfluss des 
Zufalls auf das Verfahren liegt. So liegt die durchschnittliche Länge des 
Blechs bei der Gewichtung wpa = 1, wk = 1, w, = 1 und der oben vorge- 
stellten Vorgehensweise mit einer Vorsortierung bei 1695 mm mit einer 
Standardabweichung von ungefähr 242 mm. Dieser Durchschnitt sinkt auf 
1482 mm, wenn die Startwerte des evolutionären Algorithmus rein zufäl- 
lig sind. Wenn die Vorsortierung als ein Startwert genommen wird und 
alle anderen Startwerte zufällig sind, lässt sich eine weitere Verbesserung 
auf 1388 mm beobachten. Gleichzeitig steigt jedoch die Standardabwei- 
chung auf ungefähr 442 mm bzw. 460 mm. Im direkten Vergleich von einer 
Schachtelung in den verschiedenen Verfahren ergeben sich natürlich noch 
größere Abweichungen. 


3.2.5 Lösung durch eine Verbesserungsheuristik 


Die Verbesserungsheuristik versucht von einer zulässigen Lösung aus die 
Schachtelung so zu verändern, dass nicht nur die Packungsdichte opti- 
miert wird. Die Zielfunktion bleibt grundsätzlich als additive Funktion 
wie im vorhergehenden Abschnitt bestehen. Allerdings wird hier neben 
der Packungsdichte nur die Stabilität der Teile in der Zielfunktion be- 
trachtet, da die Stegbeschädigung in der Praxis eine untergeordnete Rolle 
spielt. Das Ziel ist nun, einzelne Teile so zu verschieben, dass die gesamte 
Schachtelung im Sinne der Zielfunktion verbessert wird. 

Bei der Verbesserungsheuristik wird als erstes durch einen beliebigen 
Algorithmus eine Basislösung generiert. Der alleinige Fokus liegt dabei auf 
einer möglichst hohen Packungsdichte. Für alle Berechnungen in dieser 
Arbeit gilt, dass die Basislösung mit der Software-Bibliothek des Unter- 
nehmens Optalog SAS erstellt wurde. Grundsätzlich gilt natürlich, dass 
eine möglichst gute Basislösung erreicht werden sollte. 

Die Teile werden dann in eine Reihenfolge gebracht und nacheinander 
ggf. verschoben oder erneut platziert. Die Verschiebung findet dabei aus- 
schließlich in x-Richtung statt, da in y-Richtung kein weiterer Platz zur 
Verfügung steht. In der Basislösung wird schließlich die gesamte Breite 
des Blechs benutzt. Vor der endgültigen Platzierung werden die Teile ggf. 
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zuerst verschoben. Grundsätzlich sollen die Teile, die sich näher am linken 
Rand befinden, also bei x = 0, zuerst eingefügt werden. Jedes andere Ver- 
fahren widerspricht dem Ziel des Strip Packings. Später eingefügte Teile 
liegen dann natürlicherweise weiter rechts und müssen ggf. zusätzlich 
weiter nach rechts verschoben werden. 

Die Verschiebung hat dabei zwei Ziele. Zum einen muss sichergestellt 
sein, dass sich Teile nie überlappen und der Mindestabstand eingehal- 
ten wird. Zum anderen soll natürlich eine Position gefunden werden, bei 
der das Teil stabil durch die Auflagestege gelagert ist. Wie auch im vor- 
hergehenden Abschnitt ist dabei die Position des Freischnittpunktes ein 
bedeutender Parameter. 

Die beiden sinnvollen Methoden, um die Teile einer Schachtelung in 
eine Reihenfolge zu bringen, sind die jeweils aufsteigenden Ordnungen 
nach der kleinsten x-Koordinate eines Teils oder der größten x-Koordinate 
eines Teils. Je nach dem, wie die Teile ineinander greifen, kann dies große 
Auswirkungen auf die Verschiebungen der Teile und damit auf das finale 
Layout haben. Ein Beispiel ist in Abbildung 3.9 gezeigt. Von den oberen 
beiden Teilen wird das obere, längere Teil zuerst eingefügt, wenn die Teile 
nach der kleinsten x-Koordinate aufsteigend eingefügt werden. Wenn die 
Teile nach der größten x-Koordinate aufsteigend eingefügt werden, wird 
das untere, kleinere Teil zuerst eingefügt. Die unteren beiden Teile werden 
unabhängig von der Sortierung immer gleich eingefügt, nämlich zuerst 
das linke und dann das rechte. 


Abbildung 3.9 Vier Teile, an denen die verschiedenen Sortierungen erklärt werden können. 
Von den oberen beiden Teilen wird das obere Teil zuerst eingefügt, wenn nach der kleinsten x- 
Koordinate aufsteigend sortiert wird und das untere Teil, wenn nach der größten x-Koordinate 
aufsteigend sortiert wird. Bei den unteren beiden Teilen wird in beiden Fällen zuerst das 
linke und dann das rechte eingefügt. 


87 


3 Schachtelung für Laserschneidmaschinen 


Eine nachrangige Sortierung anhand der y-Koordinaten ist natürlich 
möglich, hat aber bei Vorversuchen zu keinen Unterschieden in den Er- 
gebnissen geführt. Dies liegt wahrscheinlich daran, dass die geometrische 
Abhängigkeit von Teilen in der y-Richtung zufällig ist. Die gewählte Sor- 
tierung stammt aus der Beschreibung der Basislösung und ist zufällig. 

Das Einfügen eines Teils besteht dann darin, zuerst die Position aus 
der Basislösung zu prüfen. Liegt an dieser Position eine Überlappung mit 
anderen Teilen vor, weil diese bereits verschoben worden sind, wird das 
Teil so lange verschoben, bis keine Überlappung mehr vorliegt. Sonst wird 
geprüft, ob das Teil stabil liegt, verkippt oder durch die Stege fällt. Falls das 
Teil verkippen sollte oder das Ergebnis der Verkippberechnung ungewiss 
ist, wird das Teil um 1 mm verschoben und wieder eine Verkippberechnung 
durchgeführt. Dieses Verfahren wird so lange wiederholt, bis entweder 
eine stabilere Position gefunden werden kann oder eine vorher definierte 
maximale Verschiebung erreicht ist. Sollte für eine Platzierung des Teils 
ohne Überlappungen mit anderen Teilen eine größere Verschiebung als 
die maximale Verschiebung notwendig sein, wird diese selbstverständlich 
durchgeführt, da das oberste Ziel ist, tatsächlich alle Teile zu platzieren. 


3.2.6 Berechnungen und Diskussion der 
Verbesserungsheuristik 


3.2.6.1 Berechnungen und Ergebnisse 


Die 20 Teiletypen und 100 Schachtelungen sind die gleichen wie in Kapitel 
3.2.3.2. Falls nicht weiter vermerkt, werden auch alle sonstigen Parameter 
beibehalten, um möglichst vergleichbare Ergebnisse zu erhalten. Da nur 
noch in eine Richtung verschoben wird, ist eine Verschiebung, die größer 
als der Abstand zwischen zwei Stegen ist, nicht sinnvoll. Die maximale 
Verschiebung ist deswegen auf 67mm beschränkt. Wird keine Position 
gefunden, an der das Teil stabil gelagert werden kann, verkippt das Teil an 
dieser y-Position immer. 

Um die Ergebnisse mit denen der Konstruktionsheuristik vergleichen zu 
können, werden hier wieder zuerst die Berechnungen mit einem zufälligen 
Freischnittpunkt und danach mit dem besten Freischnittpunkt für die 
Basislösung durchgeführt. Bei Anwendung der Verschiebeheuristik von 
Teilen ist es nicht sinnvoll, einen zufälligen Freischnittpunkt zuzulassen, 
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Abbildung 3.10 Beispiele von Schachtelungen der Basislösung der Verbesserungsheuristik 


da Teile dann weiter verschoben werden könnten, als notwendig wäre. 
Insgesamt ergeben sich also vier Berechnungsvarianten, die in Abbildung 
3.11 an einem Beispiel gezeigt sind. 

Die ersten drei der 100 Schachtelungen als Basislösung ohne Verschie- 
bungen und mit zufälligem Freischnittpunkt sind in Abbildung 3.10 dar- 
gestellt. Die durchschnittlich benötigte Länge des Blechs über alle 100 
Schachtelungen beträgt 1052,1mm. Dabei sind durchschnittlich 11,3 Teile 
stabil, 33,5 Teile kippen, 4 Teile fallen und bei 0,9 Teilen ist das Ergebnis 
der Verkippberechnung nicht eindeutig. Wie schon bei der Konstruktions- 
heuristik verbessert sich dieses Ergebnis, wenn der beste Freischnittpunkt 
gesucht wird. Hier verkippen nur noch durchschnittlich 23 Teile und 22,6 
Teile sind stabil. Die Anzahl der fallenden Teile bleibt natürlich konstant, 
während die Anzahl der ungewissen Teile auf 0 sinkt. 

Die Verschiebung bei Sortierung der Teile nach aufsteigender kleinster 
x-Koordinate führt zu einer deutlich höheren Anzahl an durchschnittlich 
stabilen Teilen, nämlich 33,7. Es verkippen noch 5,1 Teile, 10,8 Teile fal- 
len und bei keinem Teil ist die Verkippberechnung ungewiss. Durch die 
Verschiebungen steigt die durchschnittliche Blechlänge auf 1337 mm. 

Bei der Sortierung der Teile nach aufsteigend größter x-Koordinate be- 
trägt die durchschnittliche Blechlänge 1391,3 mm. Durchschnittlich sind 
33,5 Teile stabil, 4,8 Teile verkippen, 11,3 Teile fallen und bei 0 Teilen ist 
das Ergebnis der Verkippberechnung ungewiss. 
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Abbildung 3.11 Schachtelung 1 als Beispiel für die Verbesserungsheuristik. Oben links: 
Basislösung mit zufälligem Freischnittpunkt. Oben rechts: Basislösung mit besten Freischnitt- 
punkt. Unten links: Verschiebung sortiert nach aufsteigend kleinster x-Koordinate. Unten 
rechts: Verschiebung sortiert nach aufsteigend größter x-Koordinate. 


3.2.6.2 Diskussion 


Es lässt sich durch den Vergleich der Abbildungen 3.5 und 3.10 erkennen, 
dass die Basislösungen für die Verbesserungsheuristik deutlich bessere 
Schachtelungen sind, als die Konstruktionsheuristik erstellen kann. Dies 
zeigt auch der Vergleich der durchschnittlichen Blechlängen mit 1052,1mm 
und 1952 mm. 

Wie schon bei der Konstruktionsheuristik führt die Auswahl des bes- 
ten Freischnittpunkts zu einer deutlichen Verbesserung der Ergebnisse. 
Die Anzahl an stabilen Teilen wird verdoppelt. Da keine Verschiebung 
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Abbildung 3.12 Schachtelung 94 als Beispiel für einen deutlichen Unterschied der Blech- 
länge in Abhängigkeit von der Sortierung. Oben: Aufsteigende Sortierung nach größter 
x-Koordinate. Unten: Aufsteigende Sortierung nach kleinster x-Koordinate. 


stattfindet, bleibt gleichzeitig die Anzahl der fallenden Teile und die Blech- 
länge gleich. Die Verschiebung von Teilen führt zu einer deutlich erhöhten 
Anzahl von stabilen Teilen, nämlich von durchschnittlich 22,6 Teilen auf 
33,5 bzw. 33,7 Teile. Dies ist zwar auf den ersten Blick schlechter als die 
Zahlen bei der Konstruktionsheuristik, aber der Tatsache geschuldet, dass 
bei der Verbesserungsheuristik die Teile in zwei Drehrichtungen platziert 
werden können. Wenn die Summe aus stabilen und fallenden Teilen be- 
trachtet wird, steht die Verbesserungsheuristik der Konstruktionsheuristik 
kaum nach. Bei der Verbesserungsheuristik sind dies 44,5 bzw. 44,8 Teile 
und bei der Konstruktionsheuristik im besten Fall mit der Parametrierung 
wpa = 1, Wk = 10, шь = 145,1 Teile. Die Verbesserungsheuristik bietet aber 
den Vorteil von deutlich kompakteren Schachtelungen. 

Der Einfluss der Drehrichtung der Teile kann eindeutig beim Vergleich 
der Ergebnisse der beiden Heuristiken gesehen werden. Während die 
Konstruktionsheuristik die Teile immer mit der längeren Dimension in 
x-Richtung legt, können die Teile bei der Verbesserungsheursitik um 90° 
gedreht werden. Entsprechend fallen hier statt ungefähr 2,9 Teilen im 
Durchschnitt ungefähr 11 Teile im Durchschnitt durch die Auflagestege 
und haben keine stabile Position. Dass der Einfluss der Drehrichtung so 
drastisch ausfällt, hat vor allem mit der zugrunde liegenden Teilemenge zu 
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tun, die recht viele Teile enthält, die insgesamt eher klein sind (z. В. Teil 5, 
Teil 7 oder Teil 8) oder eine relativ lange und eine kurze Dimension haben 
(z.B. Teil 9, Teil 15 und Teil 19) und deswegen in einer Rotationslage fallen 
können. Der Einfluss der Rotationslage würde geringer ausfallen, wenn es 
weniger solche Teile geben würde. 

Während die Unterschiede zwischen den Sortierregeln bei der durch- 
schnittlichen Stabilität der Teile sehr gering ausfallen, gibt es einen syste- 
matischen Unterschied bei der Blechlänge. Die Sortierung nach kleinster 
x-Koordinate ist im Durchschnitt 54,3 mm kürzer bei einer Standardabwei- 
chung von 146,2 mm. Die maximale Abweichung bei einzelnen Schach- 
telungen kann bis zu 725mm betragen. Diese Unterschiede lassen sich 
dadurch erklären, dass bei der Sortierung nach größter x-Koordinate grö- 
ßere Teile später eingefügt werden und dann unter Umständen nicht aus- 
reichend Platz vorhanden ist, da die kleineren Teile um das Teil herum 
verschoben worden sein können. Das größere Teil muss dann recht weit 
verschoben werden, um ohne Überlappung platziert werden zu können. 
Ein Beispiel dafür ist die Schachtelung 94, die in Abbildung 3.12 gezeigt 
ist. 

Auffällig ist, dass insbesondere die zuletzt platzierten Teile unter Um- 
ständen eine besonders ungünstige Position einnehmen können. Dies liegt 
daran, dass die Teile nie in y-Richtung verschoben werden. Sie bleiben also 
an ihrer x-Position aus der Basislösung, was aber gegeben die verschobene 
Position aller anderen Teile eine wenig kompakte Schachtelung ergeben 
kann. Vier Beispiele dafür sind in Abbildung 3.13 dargestellt. Bei den lin- 
ken beiden Schachtelungen könnten die Teile weiter unten besser platziert 
werden. Bei den rechten beiden Schachtelungen müssten schon mehrere 
Teile neu platziert werden, um zu einer kompakteren Schachtelung zu 
gelangen. 
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Abbildung 3.13 Vier Schachtelungen, bei denen die am weitesten rechts platzierten Teile 
auch eine andere Position einnehmen hätten können, die zu einer kompakteren Schachtelung 
geführt hätte. 
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4.1 Zusammenfassung 


Die Prozesssicherheit von LFMn beim Schneiden von Blechen kann erhöht 
werden, indem die Position und der Zustand der Auflagestege detektiert 
und die geplante Bearbeitung so angepasst wird, dass mögliche Probleme, 
wie z.B. eine Kollision des Schneidkopfs mit einem verkippten Teil oder 
die unnötig starke Beschädigung der Auflagestege, vermieden werden. In 
dieser Arbeit wurden dazu drei verschiedene Arten der Auflagestegerken- 
nung vorgestellt, weiter entwickelt und neue Verfahren zur Schachtelung 
unter Berücksichtigung der Auflagestege präsentiert. 

Die Stegdetektion mittels Laser-Triangulation ist ein im Grundsatz lange 
bekanntes Verfahren. Um die Auflagestege ohne störende Reflexionen und 
möglichst genau vermessen zu können, wurden sowohl die Laserlinie als 
auch die Kamera nicht direkt über dem zu messenden Steg platziert. So 
konnten Abschattungen verhindert werden und die wenigen Artefakte 
durch Reflexionen lassen sich fast vollständig entfernen. Die Messmetho- 
de wurde validiert und konnte Genauigkeiten im Submillimeterbereich 
erreichen. 

Ein im Aufbau deutlich einfacheres Verfahren ist die Stegdetektion auf 
einzelnen Bildern im sichtbaren Spektrum. Hier wurden verschiedene 
aus der Bildverarbeitung bekannte Methoden auf ihre Leistungsfähigkeit 
bei dieser Anwendung überprüft. Kanten- und Eckendetektoren sowie 
Texturanalysemethoden konnten dabei keine guten Ergebnisse erzielen. 
Die Filterung der Frequenzen, die dem Stegabstand entsprechen, konnte 
die besten Ergebnisse erzielen. Die meisten Fehler wurden dabei in dem 
Drittel der Palette gemacht, bei dem die Auflagestege bedingt durch die 
Kameraposition von der Seite gesehen werden und Störeinflüsse z.B. durch 
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die Verschmutzung der Stege oder durch Überlappung von Spitzen und 
Senken benachbarter Stege stärker ausgeprägt sind. 

Um dieses Problem besser bewältigen zu können, wurde eine aktive Be- 
leuchtung im Nah-Infrarot Spektrum getestet. Durch eine Position, bei der 
die Stegspitzen von oben und hinten angeleuchtet werden, entsteht ein fast 
vollständig schwarzes Bild, bei dem die Stegspitzen als helle Punkte dar- 
gestellt werden. Die Vorderseiten der Stege können also keine Störeinflüsse 
mehr ins Bild einbringen, wordurch die Genauigkeit der Stegerkennung 
entsprechend steigt. Zudem konnten nun gute Ergebnisse bei der Erken- 
nung einzelner Spitzen erzielt werden. 

Die mit den Bildverarbeitungsmethoden gewonnenen Informationen 
können dann zur Optimierung des NC-Programms eingesetzt werden. In 
dieser Arbeit wurde als Ansatz für diese Optimierung vor allem die Schach- 
telung unter Berücksichtigung verschiedener Freischnittpunkte betrachtet. 
In der zuerst vorgestellten Konstruktionsheuristik wird die Bottom-Left Re- 
gel benutzt, um die Teile nacheinander auf dem Blech zu platzieren. Nach 
der ursprünglichen Platzierung wird die Stabilität und die Stegbeschä- 
digung geprüft und das Teil ggf. entgegen der Bottom-Left Richtungen 
verschoben, um eine bessere Position zu finden. Die Optimierung der 
Freischnittposition bringt alleine schon eine bedeutende Verbesserung. 
Zusammen mit dem Verschiebeverfahren konnte eine sehr große Verbesse- 
rung erreicht werden, bei der nur noch vereinzelt Teile verkippen werden. 

Bei der zuletzt betrachteten Verbesserungsheuristik werden vorhande- 
ne Schachtelungsalgorithmen benutzt, um Schachtelungen mit möglichst 
wenig Verschnitt als Basislösung zu generieren. Die Teile werden dann auf- 
steigend sortiert, entweder nach der kleinsten oder größten x-Koordinate 
und in dieser Reihenfolge platziert. Dabei wird ein Teil so lange in zx- 
Richtung verschoben, bis es ohne Überlappung platziert werden kann. 
Zusätzlich wird die Stabilität des Teils geprüft und im negativen Fall wird 
versucht, durch Verschieben entlang der «-Achse eine stabilere Position zu 
finden. Die Verbesserungsheuristik liefert ähnlich gute Ergebnisse wie die 
Konstruktionsheuristik beim Verhindern von verkippten Teilen, kann aber 
deutlich kompaktere Schachtelungen generieren, da schon die Basislösun- 
gen sehr viel kompakter sind. 
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4.2 Ausblick 


Für die in dieser Arbeit vorgestellten Methoden gibt es verschiedene Mög- 
lichkeiten zur Weiterentwicklung. 

Zur Verbesserung der Stegvermessung mittels Laser-Triangulation könn- 
ten zukünftige Arbeiten einen ähnlichen Aufbau in einer Laborumgebung 
wiederholen. Dies würde eine deutlich exaktere Bestimmung aller Mess- 
fehler und damit einer bessere Validierung ermöglichen. 

Eine weitere interessante Frage ist, die Stegbeschädigung und Verschla- 
ckung bei produktiv eingesetzten Maschinen zu vermessen. Eine solche 
Messung wird das Verständnis über die Beschädigungsprozesse und ih- 
re Einflussfaktoren verbessern. Darauf aufbauend können Methoden zur 
Verminderung der Beschädigung vorgeschlagen und getestet werden. Des 
weiteren könnte das Veständnis über die Verschlackung von Stegen auch 
neue Erkenntnisse zum Laserschneidprozess ermöglichen. 

In zukünftigen Arbeiten zur Stegdetektion auf einzelnen Bildern aus 
dem sichtbaren Spektrum sollte untersucht werden, ob mit weitergehen- 
den Verfahren der Bildverarbeitung eine exakte Detektion der Stegspitzen 
möglich ist. Denkbar ist beispielsweise der Einsatz künstlicher neurona- 
ler Netze, insbesondere von CNNs, die in vielen anderen Anwendungen 
der Bildverarbeitung außerordentlich gute Ergebnisse erzielen. Dafür ist 
allerdings eine größere Datenbasis notwendig. Es ist auch interessant, 
künstliche neuronale Netze bei der Auswertung von Bildern aus dem NIR- 
Spektrum einzusetzen. Hier stellt neben einer reinen Lokalisierung auch 
die Zustandsmessung eine realistische Option dar. 

Nicht zuletzt können weitere Einsatzgebiete einer NIR-Beleuchtung an 
einer LFM erforscht werden, wenn man die Betrachtungsgrenze anders als 
in dieser Arbeit nicht direkt nach dem Schneidprozess zieht. Beispiele sind 
die Vermessung der Blechlage außerhalb des Maschinenkörpers oder die 
Detektion verkippter Teile nach dem Schneidvorgang. Letzteres ist wichtig, 
um Kollisionen von automatisierten Absortieranlagen zu verhindern. 

Auch bei der Programmierung von LFMn gibt es Ansätze für weitere For- 
schungsarbeiten. Bei der Konstruktionsheuristik hat sich gezeigt, dass der 
in dieser Arbeit verfolgte Ansatz nur in wenigen Fällen gute Schachtelun- 
gen erstellen kann. Eine ähnlich tiefe Einbettung in einen Schachtelungs- 
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algorithmus, der das NFP als Ansatz benutzt, wäre für einen objektiven 
Vergleich mit Verbesserungsheuristiken notwendig. 

Die in dieser Arbeit vorgestellte Verbesserungsheuristik könnte noch 
kompaktere Ergebnisse erzeugen, wenn das Problem von einzelnen rela- 
tiv zum Schluss platzierten Teilen, die weit aus der restlichen Schachte- 
lung herausragen, adressiert wird. Dazu könnten diese Teile an anderen 
y-Positionen oder auch in anderen Rotationen platziert werden, möglicher- 
weise sogar durch den Schachtelungsalgorithmus, mit dem die Basislösung 
erzeugt wurde. Dazu würden einige Teile neu platziert werden, wobei die 
meisten Teile in ihrer von der Verbesserungsheuristik bestimmten Lage 
verbleiben würden. Ein weiterer Ansatz, der bei der Verbesserungsheuris- 
tik untersucht werden sollte, ist die Verschiebung der Teile in y-Richtung. 
Dafür müsste das Blech auch in dieser Dimension nicht voll belegt sein. 
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A Verkippberechnung bei 
geschnittenen Teilen 


Das hier vorgestellte Vorgehen zur Verkippberechnung von Teilen ist schon 
länger bei TRUMPF Werkzeugmaschinen GmbH + Co. KG bekannt und 
wurde schon in Abschlussarbeiten dokumentiert (siehe z.B. [70]). Der 
beschriebene Ablauf der Verkippberechnung ist nochmal schematisch in 
Abbildung A.1 dargestellt. 

Wenn Teile auf einer LFM geschnitten wurden, können sie abhängig 
von ihrer Lagerung durch die Stegspitzen und die von oben durch den 
Gasdruck in das Teil eingeleitete Kraft (im Folgenden ”Gaskraft” genannt) 
entweder durch das Restgitter nach unten fallen, im Restgitter verkippen 
oder stabil liegen bleiben. 

Um a priori zu berechnen, welchen Zustand ein Teil einnehmen wird, 
müssen vier Dinge bekannt sein: 


= Geometrie des Teils und Position des Teils relativ zu den Stegen, 
= Position des Schwerpunkts und Größe der Schwerkraft, 
= Freischnittpunkt und Größe der Gaskraft. 


Im ersten Schritt der Verkippberechnung werden die Auflagestellen des 
Teils berechnet. Dazu wird geprüft, welche Stegspitzen in der gegebenen 
Position des Teils relativ zu den Stegen das Teil stützen. Entweder das Teil 
hat keine Auflagepunkte. Dann soll im Folgenden angenommen werden, 
dass das Teil nach unten durch das Restgitter fallen wird. Dies ist in der 
Realität nicht zwangsweise der Fall, da z.B. ein Teil auch im Restgitter 
verhaken kann. Wenn das Teil genau einen oder mehrere Auflagepunkte, 
die aber alle zum gleichen Steg gehören, hat, wird das Teil verkippen. 
Entweder wird das Teil unter Einfluss der Gaskraft verkippen oder nach 
dem Schneiden ohne Einfluss der Gaskraft. Es ist theoretisch möglich, dass 
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der Schwerpunkt des Teils und der Angriffsort der Gaskraft genau über 
dem Auflagesteg sind und das Teil somit stabil liegen bleibt. Bei diesem 
Fall handelt es sich offensichtlich um ein labiles Gleichgewicht, das nicht 
lange Bestand haben würde. Deshalb werden Teile, deren Auflagepunkte 
zu einem Steg gehören, als verkippend klassifiziert. Wenn ein Teil mehr 
als zwei Auflagepunkte hat, die zu verschiedenen Auflagestegen gehören, 
muss zur Verkippberechnung erst geprüft werden, ob das Teil ohne Einfluss 
der Gaskraft stabil ist. Dazu muss sich der Schwerpunkt des Teils innerhalb 
der konvexen Hülle der Auflagepunkte befinden. Wenn das nicht der Fall 
ist, wird das Teil als verkippend klassifiziert. 

Wenn sich der Schwerpunkt des Teils in der konvexen Hülle der Aufla- 
gepunkte befindet, muss als letzter Schritt geprüft werden, wie sich das 
Teil unter Einfluss der Gaskraft verhält. Falls sich der Angriffspunkt der 
Gaskraft ebenfalls in der konvexen Hülle der Auflagepunkte befindet, ist 
das Teil stabil. Wenn nicht, muss die Drehachse des Teils bestimmt werden 
und die Drehmomente, die durch die Schwerkraft und Gaskraft ausgelöst 
werden, müssen berechnet werden. Die Drehachse des Teils kann bestimmt 
werden, indem der Schnittpunkt der Geraden zwischen dem Schwerpunkt 
und dem Freischnittpunkt und der konvexen Hülle der Auflagepunkte 
berechnet wird. Die Seite der konvexen Hülle, auf der dieser Schnittpunkt 
liegt, ist auch die Drehachse des Teils. Falls der Schnittpunkt auf einer Ecke 
der konvexen Hülle liegt, werden beide angrenzenden Achsen jeweils ge- 
prüft und das Ergebnis mit der höheren Kippwahrscheinlichkeit gewertet. 
Das Drehmoment, das durch die Schwerkraft ausgelöst wird, lässt sich aus 
der kürzesten Distanz dp, zwischen dem Schwerpunkt des Teils und der 
Drehachse, sowie der Schwerkraft des Teils Fs bestimmen. Der Hebelarm 
der Gaskraft ist die kürzeste Distanz dp „ zwischen dem Freischnittpunkt 
und der Drehachse. Die Gaskraft Ес wird durch den Durchmesser der 
Schneiddüse dp und den Gasdruck рс abgeschätzt: 

Foz: E "Ра 

Das Gesamtdrehmoment des Teils ergibt sich aus der Summe der beiden 

Drehmomente: 

Maes = Fs: dr; — Ес ‘dp. 
Dabei wird die Richtung des durch die Schwerkraft ausgelösten Drehmo- 
ments als positive Richtung definiert. Da der Angriffspunkt der Schwer- 
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kraft in der konvexen Hülle der Auflagepunkte liegt und der Angriffspunkt 
der Gaskraft außerhalb, ist das durch die Gaskraft ausgelöste Drehmoment 
negativ. 

Ist M,., positiv, überwiegt das Drehmoment der Schwerkraft und das 
Teil bleibt stabil liegen. Ist M,., negativ, überwiegt hingegen das Drehmo- 
ment der Gaskraft und das Teil verkippt. Um die Ungewissheit, die durch 
verschiedene Abschätzungen in diesem Verfahren eingeführt wird, etwas 
abzufangen, wird die Stabilität des Teils als unsicher klassifiziert, wenn 
М,» nahe an null ist. Dies wird durch eine festzulegende Schranke M, 
getan. Der Kippzustand eines Teils gilt als unsicher, wenn |М,.;| < My gilt. 
In allen Auswertungen dieser Arbeit wurde M, = 10 Nmm gewählt, was in 
vorhergehenden Arbeiten experimentell als angemessener Wert ermittelt 
wurde [70]. 
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Abbildung A.1 
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Um die erwähnten Berechnungen durchzuführen, wurden die auf den 
nächsten beiden Seiten dargestellten 20 Teile benutzt. In der darauf folgen- 
den Tabelle ist dargestellt, wie oft ein Teil in einem der 100 Schachtelungen 


vorkommt. 


Teil 6 
Teil 2 wi 
Teil 5 
oO 
Teil 14 
Teil 4 Teil 3 | T. 
Teil 11 


Teil 9 
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